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The behaviour of a Pitot tube in transverse shear 


By P. O. A. L. DAVIES 


Aeronautical Engineering Department, University of Southampton 


(Received 30 August 1957) 


SUMMARY 

‘This paper presents an experimental investigation of the 
displacement of the total-head profile in a shear flow from its true 
position when it is measured with a flat-ended round Pitot tube. 
The experiments showed that at least two forms of displacement 
may exist depending on the flow conditions. 

In a wake the displacement of the effective centre of pressure 
towards the higher velocity was found to obey the law deduced 
theoretically by Hall (1956), provided that changes in flow 
direction were not large. In the presence of large eddies the 
behaviour of the Pitot is further modified by its characteristics 
in yaw. ‘This introduces a total-head error reducing the observed 
total head. 

In a turbulent boundary layer the experiments showed that 
the displacements were an order of magnitude smaller than in 
a wake and could be regarded as negligible for a round Pitot, even 
when it was touching the wall. 

An empirical method for correcting the observations made 
with a round Pitot which is of the same width as the wake is given 
in an appendix. 


1. INTRODUCTION 

In the experimental studies of boundary layers and wakes, the velocity 
profile is generally established by a Pitot traverse. The displacement of 
the effective centre of pressure recorded by a total-head Pitot tube in a 
wake has been studied experimentally by Young & Maas (1936) and 
theoretically by Hall (1956). A more exact theoretical study has just been 
published by Lighthill (1957) and his results substantiate Hall’s simplified 
theory. 

Young & Maas studied the wake close behind a symmetrical aerofoil, 
and found a displacement towards the region of higher velocity that could 
be expressed as a constant fraction of the Pitot diameter and was independent 
of the velocity gradient. However, as Young & Maas pointed out, this 
implies a sharp discontinuity in the displacement correction as the velocity 
gradient approaches zero, which is not physically plausible. ‘Their observa- 
tions indicate a velocity correction at the centre of the wake, but they did 
not discuss this in their paper. 
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Halls’ theoretical investigation was based on the consideration of the 
vorticity field alone and neglected viscosity. He pointed out that in front 
of a three-dimensional object in uniform shear the vortex tube must be 
stretched, giving a local increase in vorticity. ‘This gives a larger displace- 
ment of the stagnation streamline than pure two-dimensional flow alone. 
He then calculated the displacement of the stagnation streamline for a 
sphere in uniform shear, and showed it was a function of the velocity 
gradient as well as the mean velocity and the diameter of the Pitot. This 
avoids the discontinuity described in the previous paragraph. 

‘This paper first describes an experimental investigation of the validity 
of Hall’s assumptions, to resolve the discrepancy with the results of Young 
& Maas. Further experiments are described which were undertaken to 
establish the form of corrections which apply to other examples of shear 
flow. ‘lhe experiments indicated that the different conditions existing in 
other forms of shear flow introduce different corrections. Examples of 
shear flow where conditions may differ significantly include (a) flow in a 
wake after large eddies have decayed; (6) flow in a wake just downstream 
of the trailing edge; (c) flow in a turbulent boundary layer; (d) flow in 
a laminar boundary layer; (e) flow where compressibility is important; 
(f) pulsating flow. 

‘These experiments were confined to the first three of the above examples. 
Hall’s work is confined to the first shear flow example, which is discussed 
in §2 of this paper. ‘lhe results of Young & Maas may be classed with 
example (6), which is discussed in §3. Flow of the type in example (c) is 
discussed in § 4. 

Other experimental work on displacement corrections has been published. 
Livesey (1956) has studied the corrections for cylindrical Pitots, and also 
gives some results for flat-ended Pitots of the conventional pattern. All 
his results were obtained in a turbulent boundary layer after it had been 
thickened by a wire cloth screen. Mawson & Lilley (1956) include results 
for high speed flow, and also include some observations for low speed flow 
which may be classed with examples (a) or (6). They were, however, 
primarily concerned with corrections for total-head tubes larger than the 
wake. Macmillan (1957) gives careful measurements which may be classed 
with example (c) both for a pipe and a flat plate. The effect of flow 
pulsations on total-head tube readings has been investigated by several 
authors, although very few observations have been recorded. 


2. 'TOTAL HEAD MEASUREMENTS IN A WAKE 


Hall (1956) has shown that, when the vorticity field alone is considered, 
the displacement 6 of the stagnation streamline from its position at infinity 


for a sphere in uniform shear can be expressed by 

5/D = 0-6200K — 0-5786K? + O(K*), (1) 
where D is the sphere diameter, and K represents the maximum difference 
between the velocities at the edges of the sphere divided by twice the mean 
velocity. (Lighthill’s (1957) asymptotically exact theory for flows with 
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small shear (i.e. small K) gives 6/D = 0-45K—CK®, where Cis probably 
greater than 1:35. In his paper, K is expressed as Aa/u.) Hall suggested 
that a conventional flat-ended round Pitot tube may be represented by a 
sphere of slightly larger diameter. 

It still remains to test Hall’s two assumptions, namely that the viscosity 
may be neglected, and that a round Pitot of slightly smaller diameter may 
replace the sphere. ‘lo this end, traverses were made through a wake by 
a series of flat-ended round Pitots of graded size but constant shape, the 
ratio of bore to outer diameter being 0-6. Other experiments were made 
to investigate the effect of changing this ratio. 

The traverses were carried out behind two models in a N.P.L. wind 
tunnel of 3 ft. x 2 ft. working section, the maximum available wind speed 
being 70 ft./sec. One of these was a two-dimensional symmetrical aerofoil 
of 10 in. chord and 5°, thick. ‘lo permit the use of larger tubes an axi- 
symmetric double ogive (cigar-shaped), of 40 in. chord and 7 in. maximum 
diameter, was also employed. ‘This form was chosen for the larger model 
to minimize the effects of tunnel interference. It was considered that at 
radii large compared with the Pitot diameter the flow could be considered 
for practical purposes two-dimensional. For checking Hall’s assumptions, 
the traverses were carried out several wake widths downstream to allow 
for the decay of the large eddies produced by oscillation of the rear 
stagnation point, particularly in two-dimensional flow. Most of the 
measurements were made at a fixed wind speed, though check traverses 
made at other speeds gave similar results. 

The undisturbed velocity head py was measured with a_ standard 
round-nosed N.P.L. Pitot-static combination, of 0-5 in. outside diameter, 
placed well clear of the wake and just downstream of the trailing edge. 
The reading p of the total-head tube under test was compared with that 
from the standard total-head tube tapping, the difference (p)—p) being 
observed. Care was taken to ensure this difference was exactly zero when 
the tube under test was outside the wake. ‘This provided a check on the 
uniformity of the flow in the tunnel, and of the absence of leaks. The 
tunnel was an open N.P.L. type, and it was found that any leaks into the 
tunnel disturbed the flow there considerably. 

The measurements were made in regions of appreciable turbulence, 
and therefore the pressure differences observed within the wake were 
fluctuating. ‘l'o investigate a Pitot’s behaviour it was essential to establish 
that the mean observed pressure was unaffected by changes in the dynamic 
behaviour of the manometer system due to changes of Pitot bore. It would 
have been ideal if the manometer system always recorded the true mean 
Pitot pressure, but it was not possible to establish that it did so in turbulent 
flow. Experiments did show, however, that variation of the manometer 
system response by throttling the leads produced no change in the recorded 
mean pressure in turbulent flow. ‘This indicated that for the purpose of 
comparison of Pitot performance, variation of the dynamic behaviour of 
the manometer could be neglected. 
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Small temperature variations also limit the minimum Pitot bore which 
may be used satisfactorily with a given manometer system in unsteady 
conditions. It only requires a temperature change of 0-01°C to produce 
pressure variations of 0-015 in. of water in a constant volume system at 
atmospheric pressure. When the Pitot bore was 0-018 in., the manometer 
system would respond to a small disturbance of approximately } second 
period. When the bore was reduced to 0-005 in., however, the period 
rose to 40 seconds, which was considered too long. 
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Figure 1. Wake velocity profiles observed by geometrically similar Pitots behind 
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the aerofoil model. 


‘The manometer was a simple vertical U-tube made of 3 mm precision 
bore glass tubing. A small quantity of detergent was added to the manometer 
water to minimize the effect of surface tension variations. ‘The level was 
read with a travelling microscope fitted with a graticule to facilitate 
observation of fluctuating readings. Careful calibration showed that 
readings consistent to 0-001 in. of water were possible, though the precision 
could be increased by using a higher power microscope. When the amplitude 
of the fluctuations exceeded 0-010 in. of water, spot readings were taken 
every five seconds for half a minute and averaged to give the mean pressure. 

‘Typical results are shown in figures 1 to 3, where the dimensionless 
wake velocity (p/py)!? is plotted against transverse displacement y. From 
graphs of this type it was not possible to obtain the displacement correction 
and the corresponding value of A until the undisturbed profile was known. 
This was obtained by extrapolation from the profile measured with the 
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2. Wake velocity profiles observed by geometrically similar Pitots near the 
trailing edge of the axi-symmetric model. 
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smallest tube. When K was small, the corrections for the smallest Pitot 
were negligible (less than 0-001 in. which was the limit of accuracy for 
measurements of y) compared with those for the largest tube, which were 
of the order of 0-05 in. Once the correction for low K was established, 
it was used to obtain the undisturbed velocity profile by extrapolation from 
that measured by the smallest Pitot. ‘This prejudiced the accuracy of the 
results at high A values (K 0-2, say), but it was also difficult to obtain 
constant velocity gradients of sufficient width for measurements 1n this 
region. 

‘The displacements observed have been plotted against the corresponding 
values of K in figure 4. This only includes those observations where the 
velocity gradient across the face of the tube was practically uniform. It 
can be seen from figure 4, particularly where K < 0-2 and the accuracy 
is reasonably good, that the results show reasonable agreement with Hall’s 
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Figure 4. Pitot displacement corrections in two-dimensional uniform shear. 


theoretical expression. ‘The initial slope indicates a larger correction for 
a Pitot than for a sphere of the same diameter as Hall suggests. Lighthill 
(1957) has shown that the coefficient of — K? in Hall’s expression (1) is 
much too small, which is in agreement with the results for flat-ended Pitots 
in figure 4. If Lighthill’s expression is multiplied by two it fits these results. 

Hall’s theory applies to the displacement of the stagnation streamline 
in uniform shear flow. It might seem reasonable to suppose that an 
asymmetrical pressure distribution exists over the face of the tube corre- 
sponding to the incident flow. It is desirable therefore to check whether 
variations of the ratio of bore to diameter will alter the observed pressure. 
To this end a set of Pitots of constant diameter but varying bore were 
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tested, the results being shown in the upper part of figure 5. It can be seen 
that no significant difference was observed even though the ratio of bore to 
diameter was varied from 0-075 to 0-800. ‘The observed pressure was, 
however, found to vary when the measurements were repeated close 
behind the trailing edge, as is demonstrated later in this paper. 

In many instances (e.g. in supersonic flow) it is not possible to employ 
a Pitot small compared with the width of the wake. Results for traverses 
made with large Pitots, the largest spanning the wake, are included in 
figures 1 and 2. ‘The observed displacements are less than those obtained 
by applying Hall’s theory using the slope of the undisturbed profile at the 
centre of the Pitot. ‘his is to be expected since Hall’s theory applies to 
uniform shear and here both the sign and magnitude of the vorticity changes 


continuously across the face of the tube. 
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Figure 5. Effect of changing Pitot bore on observed wake velocity profile. ‘The 
lower curve gives the corrected wake velocity profile. 


It is, however, possible to obtain an approximate correction from a 
traverse witha large tube. ‘The results in the lower part of figure | show that 
the corrections are zero at the centre of the wake. If the velocity in this 
position is taken as the minimum velocity, and if the maximum velocity 
is taken as unity, A being determined from these values, the corrections 
obtained are 6 = 0-031 (AK = 0-058) compared with 0-033 observed with 
the large tube at the mean velocity. ‘This correction is a little smaller than 
that observed, but not seriously so. It is well known that the shape of the 
velocity distribution in the wake obeys a universal law (Goldstein 1938, 
p. 581). If the velocity at the centre of the wake is known, and the width 
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of the wake is also known at the mean wake velocity, it is possible to 
establish the whole profile. ‘lhe details of the process are given in the 
appendix to this paper, and the corrected profiles so obtained for the big 
Pitot are plotted in the lower part of figure 5. It can be seen that this 
empirical process can yield satisfactory results under favourable conditions. 
3. WAKE TOTAL-HEAD MEASUREMENTS CLOSE TO THE TRAILING EDGE 

The results obtained when traverses were made close to the trailing edge 
differed from those already described. It was found that the profiles did 
not coincide at the centre of the wake, where the displacement correction 
should be zero as in figure 1. ‘This implies a total-head correction also 
exists, as was indicated by the results obtained by Young & Maas (1936). 
The results obtained just behind the symmetrical aerofoil were all similar 
to those plotted in figure 6. It can be seen here that the displacement is 
proportional to the Pitot diameter, although it differs on each side of the 
wake. 
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LT 1 j 1 1 1 1 
Young & Maas obtained the correction by plotting the observed wake 
width at a given velocity versus the Pitot diameter. “They obtained points 


i 


lying approximately on straight lines of constant slope, the correction being 


equal to half this slope. Using their technique a displacement correction 
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equal to 0-24D was obtained for the results in figure 6. However, traverses 
3/16 in. closer to the trailing edge yielded a correction equal to 0-28D, 
while traverses an equal distance downstream gave a correction of 0:16D 
when reduced by the same method. 

‘The matter can be carried further. ‘The momentum thickness calculated 
for each profile measured behind the aerofoil has been plotted against 
distance downstream from the trailing edge in figure 7. It can be seen that 
there is a tendency for the observed momentum thickness to rise sharply 
as the trailing edge is approached. If the correction were simply one of 
displacement towards the higher velocity one would expect the curve to 
drop as the trailing edge is approached, not rise. ‘The effect could be 
explained by the existence of a velocity-head correction close to the trailing 
edge. A reason for its existence is not hard to find. 
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observed by geometrically similar Pitots in the aerofoil wake. 


It is well known that a Pitot tube is very sensitive to yaw of more than 
a tew degrees (‘Todd 1949). Not only is the mean direction of the flow 
inclined inwards near the trailing edges, but the flow direction oscillates 
rapidly about the trailing edge. ‘This will produce large rapid variations 
in incident flow direction, and due to its directional characteristics the 
Pitot will record a pressure lower than the mean total head. Further 
downstream the fluctuations in the incident flow direction will be much 
smaller, and the Pitot will then record the mean total head. 
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It can be seen from figure 6 that this error in the total-head reading 
increases as the size of the Pitot is reduced. One would expect a Pitot close 
to the trailing edge would modify the flow there. ‘Che undisturbed flow is 
two-dimensional, but when a Pitot is introduced, the flow in its neighbour- 
hood must be considered three-dimensional. In addition, the presence of 
a rigid body wiil introduce extra shear stresses if the flow direction departs 
from that of the body. ‘This would modify the fluctuations in flow direction 
in the vicinity of the Pitot. ‘lhe stabilizing of the flow would increase with 
the size of tube, with a corresponding reduction in the velocity-head error. 
[he observations indicate that the reduction in the velocity-head error is 
roughly proportional to the diameter of the tube. 
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Figure 8. The effect of changing the bore of a yawed Pitot on its total head response. 


‘This explanation of the results requires further confirmation. ‘lo this 
end, use may be made of the behaviour of yawed Pitot tubes. A series of 
Pitots of constant diameter but varying bore were tested in an empty tunnel 
at various angles of yaw. ‘The results obtained are plotted in figure 8, which 
shows observed total head plotted against yaw angle. It can be seen that 
changes in bore produced substantial changes in total-head error. ‘This 
seems reasonable as the forward stagnation point must be close to the 
leading corner of the Pitot and a pronounced pressure gradient must exist 
across the tapped face. Changes in bore will thus vary the observed pressure 
at any given yaw angle as the results indicate. 

It has been shown earlier in figure 5 (upper part) that a series of Pitots 
of constant outside diameter but varying bore have the same displacement 
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correction in a stream of uniform vorticity. Furthermore, it seems reasonable 
to assume that their influence on the flow around them must be the same. 
‘The results of traverses made by a series of such Pitots are plotted in figure 9. 
In this case any difference between the profiles should be due to changes 
of total-head error of the form indicated in figure 8. ‘The results indicate 
that substantial total-head errors do exist particularly at the centre and edges 


of the wake where the largest fluctuations are known to occur. 
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Figure 9. Wake velocity profiles displaying total head errors due to yaw. 


A further experiment was performed to ensure the results in figure 9 
were independent of the manometer system response. ‘The 0-098 in. bore 
Pitot was fitted with a drilled plug to reduce its bore to 0-018 in. ‘T'raverses 
with the plug flush with the Pitot face reproduced the results obtained with 
the normal 0-018 in. Pitot, as would be expected. When the plug was 
pushed back, however, the behaviour approached that of an unplugged 
Pitot. ‘The results are plotted in figure 10 and agree with those in figure 9, 
When the plug was displaced more than 0-050 in. the Pitot behaved as it 
it was unplugged. ‘lhe dynamic response of the manometer system must 
have remained unchanged during these tests, demonstrating that the 
results in figure 9 were independent of the response. 

It should .be possible to obtain a profile corrected for velocity-head 
error corresponding to the results in figure 9, which are mean values. ‘lhe 
results in figure 8 show how the correction varies with incident flow angle, 
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but it would be necessary to determine the distribution of incident flow 
angles with time before a mean correction could be calculated. An empirical 
correction can be estimated if it is assumed that the velocity-head error 
approaches zero as the Pitot wall thickness does so. This is true for small 


angles as shown in figure 8. 
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Figure 10. Wake velocity profiles observed with a manometer system of constant 
response obtained by fitting a sliding plug in the Pitot bore. 


If this is done for the results in figure 9, the profile is very nearly that 
shown recorded with the largest bore Pitot, there being small corrections 
at the centre and edges of the wake. The profile so obtained must now be 
corrected for displacement error, and as the Pitot is large compared with 
the wake width the method outlined in the appendix is employed. 

he profile finally corrected gives a momentum thickness of 0-025 in. 
(instead of the uncorrected value of 0-015 in.), which compares favourably 
with the values 0-025 in. and 0-0245 in. observed further down the wake using 
the smaller Pitot of 0-031 in. diameter. 

‘The existence of a velocity-head error accounts for the discrepancy 
between Hali’s results and those of Young & Maas. ‘The presence of 
velocity-head corrections at the centre and edges of the wake, where the 
velocity gradient is small, explains why Young & Maas obtained large 


displacement corrections in this region. 
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4. ‘T'RAVERSES IN A TURBULENT BOUNDARY LAYER ALONG A FLAT PLATE 

In boundary layer flow it seems reasonable to expect that the results 
described above will be considerably modified by the presence of the wall. 
This will tend to reduce the displacement of the stagnation streamline as 
it will modify the flow around the Pitot. For the same reason one would 
expect that the random changes in flow direction will be much smaller. 
To illustrate this one may consider the extreme case of a Pitot in contact 
with a wall. 











spgcchaic gpcinnaisahegnates _ . oe 
|-OOF —— 
e ot 
-. 
| ill 
0:90! —* 
| al 
| ra 
: ge 
O-80- se 
> rh 
it - 
U j y 
fe) | : 
= a 
w | a PITOT DIAMETER 
O70F 2 © 003! IN 
8 @ 0060 JIN. 
I 4 ® 0090 IN 
| 7 ® O120 IN. 
| ¢ © 0240 IN 
O-60r J 
| ® 0480 IN 
i? 
aid. a ea a re i 
oO. Oo! 02 03 04 O5 0-6 0:7 


DISTANCE FROM WALL (IN.) 


Figure 11. Comparison of the measurement of a turbulent boundary layer velocity 
profile by a series of geometrically similar Pitots. 


‘Traverses made through the boundary layer near the trailing edge ot 
a flat plate of 3 ft. chord at zero incidence confirmed the above suggestions. 
The results are plotted in figure 11, where it can be seen that no appreciable 
displacement or velocity-head error was observed, although the diameter 
of the largest Pitot was 16 times that of the smallest. ‘This result was 
unexpected as it is reasonable to suppose some displacement would exist. 
A small displacement, 0-02 of the Pitot diameter, was recorded by the largest 
Pitot tube. 

There is some evidence (Preston 1954) that corrections exist close to 
the wall for flat Pitots in a boundary layer. ‘The results in figure 11 indicate 
that the correction for a round Pitot is negligible in this region, even when the 
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Pitot is touching the wall. Working with round Pitots, Macmillan (1957) 
observed displacements of the order of 0-15D in a pipe, and of the same 
order in the boundary layer on a flat plate. His results also implied a 
total-head correction close to the wall which was expressed as a negative 
displacement. ‘The absence of such corrections in the present results may 
be due to the higher response rate of the manometer system used. 

Livesey (1956) obtained corrections for a flat-ended Pitot in a boundary 
layer that had been thickened considerably by a wire cloth screen. It is 
possible that the flow conditions in this case were more nearly those 
obtaining in a wake than those in a boundary layer. Livesey also reported 
that the corrections for a conical Pitot with a sharp lip were negligible 
indicating, that his other corrections were probably of the total-head type 
and not displacement corrections as he suggests. Assuming this is so the 
results in figure + can be applied to his measurements with the sharp-lipped 
Pitots. For his largest (0-7 in. diameter) Pitot at the centre of the velocity 
profile, we get K = 0-2 giving a displacement of the order 0-03D which is 
less than the scatter of the observations. 

The results of both Livesey and Macmillan indicate that both types of 
corrections may be observed in a turbulent boundary layer. ‘The conditions 
under which these corrections will occur have not been established. 


5. CONCLUSIONS 

The results indicate that the displacement of the effective centre of 
pressure recorded by a total-head tube in a wake obeys a relationship of the 
form put forward by Hall (1956) and by Lighthill (1957). ‘This is expressed 
by equation (1), for a sphere in uniform shearing flow, and by the curve 
in figure 4 for a flat-ended round Pitot of conventional form. 

When the Pitot is large compared with the dimensions of the wake, the 
displacement is reduced and may be determined approximately by the 
process outlined in the appendix. 

In the region close behind a trailing edge, a velocity-head error also 
exists due to the characteristics of a yawed Pitot. Since this error is a 
function of the dynamic characteristics of the flow, it cannot be specified 
in terms of the Pitot geometry and the velocity gradient alone. It is possible 
to obtain a reasonably accurate corrected profile empirically, although the 
process is tedious. ‘The results suggest that, for accurate drag measurements, 
it is best to make wake traverses well downstream of the trailing edge, and 
an advantage to use a Pitot with a sharp lip. 

For a flat-ended round Pitot of conventional form and with bore 0-6 times 
the outer diameter, no corrections were observed in a turbulent boundary 
layer, even when the Pitot was touching the wall. 


‘The majority of the experimental work was carried out in the Mechanical 
Engineering Department of the University of Adelaide. 
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APPENDIX 


The estimation of displacement corrections for a Pitot large compared with the 
width of the wake 





Schlichting (1930) has established that the velocity profile in a wake 


obeys the relation 
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where the notation is illustrated in figure 12. 
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Figure 12. Illustration of notation used for displacement corrections to wake velocity 





profiles. 


The value of U,,, is given by the Pitot traverse regardless of the Pitot 

size, and y, can be found from (A 1) given the value of y corresponding to 
a known value of U. ‘The most convenient for this purpose is y,, corre- 

| sponding to the mean wake velocity U, = 3(U,+U,,), and (A1) shows 
that vo = 0-441y,. 

| ‘The value of y, must first be corrected for displacement. If the tube 
is small, the corrections can be determined from figure 4 using the 
appropriate value of AK. When the Pitot diameter is of the same order 
of magnitude as the wake width, it has been shown experimentally that 

Uy — U.,, 

2U, 
The corresponding value of 6/D may be read from figure 4. 
Once yy is known, the whole profile may be calculated from (A 1). 
‘This profile will be in error near the edge of the wake, but it can be seen 


K= 
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from the lower part of figure 5 that otherwise the agreement is good. It is 
not suggested that this empirical procedure may be applied to a case where 
the Pitot diameter is very much larger than the wake width. Such a 
condition is unlikely to arise in practice; but, if so, it is not difficult to 
establish a suitable correction experimentally. 
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SUMMARY 

This paper describes theoretical and experimental studies 
of the effects on shock tube flows of a monotonic convergence 
at the diaphragm section. Systematic flow equations are developed 
for tubes of uniform bore and tubes having either a monotonic 
convergence or a convergence—-divergence in the diaphragm 
section. Except across the shock front itself, isentropic processes 
and ideal-gas behaviour have been assumed. Simplified pro- 
cedures are presented for predicting the ideal-flow parameters 
over a wide range of operating conditions, as well as for comparing 
straight and convergent tubes. Such comparisons made by 
other investigators are found to be incomplete or in error. ‘The 
experiments described utilize a very simple device for altering 
the diaphragm section convergence and a multi-station measure- 
ment of shock velocity. ‘The expected effect of convergence is 
verified over a wide range of Mach numbers. Even at Mach 
numbers where the processes of shock formation can no longer 
be ignored, it is found that the relative performance between a 
uniform and convergent tube is preserved. 


INTRODUCTION 


The basic shock tube is one of uniform cross-section throughout, in 
which the rapid removal of a diaphragm separating two gases at different 
pressures leads to the generation of a shock wave in the region of lower 
pressure. The regions initially at high and low pressure are designated 
as driver and driven sections respectively. It has now been well established 
experimentally that if all other initial conditions are equal, one will realize 
a stronger shock wave in a tube having in the diaphragm region an area 
reduction from driver to driven section than in a tube of uniform bore. 
So-called convergent shock-tube flows have been considered by Bannister 
& Mucklow (1948), Lukasiewicz (1952), Resler, Lin & Kantrowitz (1952), 
Wallace & Mitchell (1953), Wallace & Nassif (1954), Hooker & White (1955), 
and Alpher & White (1957). ‘This last reference contains the theoretical 
work of this paper together with more extensive derivations and calculations. 

* A preliminary account of this work was given at the 1957 Annual Meeting of the 
American Physical Society in New York City. 


F.M. 2G 








458 R. A. Alpher and D. R. White 


In this paper we shall develop a general idealized theory of shock-tube 
flows, valid for uniform, convergent, or convergent-divergent tubes, 
whether operated in the subsonic or supersonic cold-flow regimes. 
Experiments will be described which show the difference between uniform 
and convergent tubes to be as predicted over a wide range of Mach numbers, 
even though with increasing shock strength an increasing deviation from 
absolute predictions of the theory is noted. 
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Figure 1. Schematic diagram of shock tube with convergent geometry at the dia- 
phragm section. The dashed lines in the diaphragm section denote a 
convergent—divergent geometry with minimum area at 3h’. The pressure 
distribution is, of course, only symbolic 


The model of the flow situation which has been employed is given in 
figure 1. After diaphragm removal, energy is extracted from the driver 
gas through an unsteady expansion from state 4 to state 3a. The transition 
section presents three possibilities. If it is uniform, te. if it connects 
equal areas, states 3a and 36 coalesce and denote the state at the diaphragm 
location. If the cold flow is supersonic, this state is the point in an unsteady 
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expansion at which the Mach number becomes unity. If the transition is 
montotonically convergent, the flow from 3a to 36 is steady nozzle flow, 
and the convergent flow at 3b may be either subsonic or sonic. Following 34 
one has either uniform flow in state 3 with M, | or an unsteady supersonic 
expansion to state 3 with M, l. Finally, if the transition section is 
convergent—divergent, i.e. if 3b’ represents a minimum area, one has either 
a subsonic or supersonic nozzle between states 3a and 36. In the supersonic 
case the exit flow at state 3b may be subjected to a further unsteady 
expansion to state 3. Boundary-layer effects, attenuation and the process 
of shock formation have not been considered. ‘lhe driver gas is treated 
isentropically, and departures from ideality are not explicitly considered 
for the shocked gas. As will be seen, it is a simple matter to modify this 
presentation to include the real changes in the driven gas states. 

Basic calculations involve determining, for initially given gases and 
thermodynamic states in the driver and driven sections (i.e. the pressures p, 
and p,, sound velocities a, and a,, and specific heat ratios y, and y,) the 
resultant shock strength and the thermodynamic state of both the expanded 
and shocked gases. 


‘THEORY 

Uniform shock tube 

The situation in a uniform tube is sufficiently well known that it need 
have only brief mention. One wants to relate the initial pressure ratio 
across the diaphragm, ~,/p, = 3, to the pressure ratio across the shock 
(i.e. the shock strength p,/p, = y) or alternatively to the shock Mach 
number M,. In the uniform tube there is a continuous unsteady expansion 
from state 4 to state 3. If M,; > 1 one has sonic flow at the diaphragm 
location. ‘The particle velocity imparted by an unsteady expansion from p, 
to ps3, where we note that p = kp’ and that u, = 0, can be written (see 
Liepman & Puckett (1947), p. 461) 


° (dp 12 dp 2a, ( p3\ ed 2" 
Us io) - = ae | — ( “al ; (1) 
J 0s \4 p p V4 Ps 


This velocity can be subsonic, sonic or supersonic in value. ‘The flow or 
particle velocity behind the normal shock is determined from the shock 
Hugoniot; or, if treating the shocked gas as ideal gives an adequate 
approximation, then the Rankine—Hugoniot relations give 


Us 21/2(y — 1) 2 M2 


| 

SG “i iwi <a ia (2) 
ay (YiL(1 = (v1 T y}j 2 wild 8 

We can combine (1) and (2) since at the interface between shocked and 

expanded gas we must have us, = u, and ps = p.. Solving for = yields, tor 

the uniform tube, 
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Tube with diaphragm section area change 

Relationships between = and y or M, for the tube with an area change 
in the diaphragm section are obtained as for the uniform tube. Referring 
again to figure 1, one proceeds as before to connect the shock front with the 
driving gas by matching pressure and flow velocity at the contact surface, 
except that now steady or unsteady processes, as the case may be, must be 
considered at the diaphragm section. Consider the general case of a 
convergent—divergent diaphragm section. When the shock-tube flow is 
established, the pressure ratio x may be expanded as follows (see figure 1): 


_, — Ps _ Pa Poo Pav Pay Bo Be " 

“Pi Psa Pav Pav Ps Pa Pr’ 
Here the sufhx c refers to the non-uniform tube, whereas suffix s will be 
used to refer to the uniform tube. ‘The various pressure ratios in (4) can be 
interpreted directly. ‘Thus p,/ps,, is that ratio required to accelerate the 
driver gas by an unsteady expansion from rest to Mach number M,,.. The 
ratio ps,,/P3, is that required to bring the driver gas by a steady expansion 
from M;,, to M;,,.. One can introduce the area ratio A,/A4,, into p5,/P3y by 
using the condition of conservation of mass. Expansion from 30’ to 3b also 
will be a steady process, supersonic or subsonic according as the flow at 36’ 
has become sonic or not. (By our restriction to isentropic channel flow 
for the driver gas, we exclude that range of z which leads to an over- 
expansion for a convergent—divergent section.) In either event the ratio 
P,,,/Ps, is that required to bring the gas by a steady expansion from Msg,, 
to M;, through an area ratio As,,/As,. ‘The ratio ps,/p3 is then that needed 
for an unsteady expansion from M,, to M,. At the contact surface, ps = po, 
and p,/p, is defined as the shock strength y. Equation (4) can therefore 
be written as 


i ae Ya ‘lay 2+(¥4— Mp Wr 2+(ya—1)Mz ory 
- 4 20" WL2 + (ya DMG L2 +44 - M0 J) ~ 
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Since this relation involves M;, M;,, and M3, as well as z, and y, we require 
additional relationships. One of these involves the overall area ratio in 
the transition section, which can be written as follows, regardless of whether 
M;,, is subsonic or sonic or whether 36’ is a section of minimum area or an 








A,_ My E +(y4- om DID 6) 
A, ~ Mon L2+(y4—1)My : 
The second relation we require connects M, with M3, M;,, and M,,. One 
expands , in the same way as 2, in (4); thus 
“uy... 8.8/9.3) (3 dsp Usa 4), (7) 
GQ, G3, 4 | a, 4G] \Gsp Gaq Gy Q 


Substituting for the sound velocity ratios, we obtain 
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where the quantity g, which we call an ‘ equivalence’ factor (see next section), 
was first used by Resler, Lin & Kantrowitz (1952). It is defined as 


an (P2+(ys—1)M5, ]!7[ 24 (v1— 1) Ms, lla (9 
* LZ +@s- DM] L2+0.- DMs, : 
Using g, one can rewrite (5) as 


y A ee 27'4/(Vs 1) y Us a a; 4 | 271 O's—-D 
2, ==) i+ SM, ae) fa oe iin . (10) 
 - 2 g &@ ££ ” 


Equations (6) to (10) enable the calculation of shock-tube flows with a 
convergent or convergent-divergent diaphragm section, in the case of 








subsonic, sonic or supersonic cold flow. ‘The results for a uniform tube are 
recovered by noting that with 4,/A4, = 1, one has M,, = M,, with g = 1. 

For subsonic cold flow (the assumption of isentropy excludes a transition 
to subsonic flow through a shock), the transition section is a subsonic nozzle 
with M,, = M3, ps, = ps and ay, = a3. Equations (6), (8) and (9) reduce 
to three simultaneous equations in g, 4/, and M;,,, which can be solved by 
iteration with an ab initio choice of a,, a, and M,. In this subsonic case, 
only the ratio 4,/A, enters, and the existence of a minimum section A3,, is 
irrelevant. 

For supersonic cold flow, M, > 1, one must have M,,, = 1 (or M;, = 1 
if the section is monotonically convergent so that A,, = A, is the minimum 
area). ‘The usual relationships between area ratio and Mach number in 


) 


supersonic nozzles apply, as used in obtaining (5), to give 47, and M3). 
From these M3, g, and z, follow. 


Computation of shock-tube parameters 

‘The procedure used in developing the shock-tube equations suggests 
a rather simple graphical means of solving the equations for a given tube. 
The reader will recall that the shocked gas and the driver gas are connected 
by matching pressure and velocity across the contact surface. ‘The graphical 
material can be divided into that pertaining to the driver and that to the 
driven gases, the one being independent of the other except as matched 
at the contact surface. ‘The advantages of this are made clear in the following. 

The shocked gas may be characterized by a plot of either shock strength y 
or shock Mach number versus w./a,, where w, is the flow velocity imparted by 
the shock. Such curves may be computed from (2) for ideal gases, but 
require a knowledge of the shock Hugoniot if the shock strength is 
sufficient to merit this (see, for example, Déring 1949, Alpher 1957, o1 
Lighthill 1957). 

The driver gas is brought into the calculation in the following way. 
A statement of A,/A,, y, and a, determines the equivalence factor g as 
a function of us/a, = uy/a,. Plots of g vs (wy/a,)(a,/a,) are given in figure 2 
for y, = 5/3 and in figure 3 for y, = 7/5 for both subsonic and supersonic 
cold flow. It may be noted that for V/, 1 these plots apply either to a 
monotonically convergent or a convergent-divergent transition section, 
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while for W@, > 1 they hold only for the monotonically convergent section, 
since in this latter case M3, (see (9)) is taken as unity. Figure 4 is a plot 
of g vs y, for various values of A,/A, in the case M, > 1. Again the plot 
holds for monotonic convergence, since M3, is taken as unity. Such a plot 
is of use with combustion drivers, in which helium or hydrogen plus the 
products of hydrogen-oxygen combustion give an intermediate y, value. 








EQUIVALENCE FACTOR 7 
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Figure 2. ‘The equivalence factor g vs (us/a,)(a,/a,) for a selection of area ratios, with 

v4 1-4. The connecting line at the right of the curves is the locus of g 

values for which the cold flow is sonic (V/, 1). For larger values of 

(3/a,)(a,/ay) corresponding to MW; 1, gis constant and equal to the values on 


this sonic line. 


‘Vhe usual calculation to be performed for a shock tube 1s the deter- 
mination of the diaphragm pressure ratio z required to produce a shock of 
given strength y. One is given y,, yy, @,, 4, 44/4, and y. From y the 
How velocity u,/a, is computed. ‘The appropriate equivalence factor is 
then read from figures 2, 3 or 4. Finally from figure 5 (which is a plot of 
equation (10)) one reads gz/y, whence z is determined. It may be noted 
that the quantity 2/y is the ratio of pressures through which the driver gas 


expands. When one has at hand a tube with fixed A,/4,, then curves 
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of 3/y vs uy/a, may be prepared which incorporate the appropriate 
equivalence factor. Such curves are particularly useful when various gases 
are to be used or when high temperature shock Hugoniots are to be taken 
into account. 
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Figure 3. The equivalence factor g vs (us/a,)(a,/a,) for a selection of area ratios with 
Vs 1:667. ‘The connecting line at the right of the curves has the same 
meaning as in figure 2. 

To illustrate the results of such calculations, we have prepared figure 6, 
which gives the diaphragm pressure ratio required to produce given shock 
Mach numbers for helium driving into ideal argon (electronic excitation 
and ionization are negligible up to M = 10 in argon) for several area 
ratios. An ideal y vs M, curve is included for reference. 

‘Gain’ factors 

Although the use of a convergent transition section may have such 


advantages as smoothing out flow irregularities associated with combustion 
driving, as well as locating the diaphragm where it will not choke the flow 





if it does not open completely, the main reason for using such sections is 
that they provide an easy way of getting a stronger shock, other things 
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being equal. It is therefore pertinent to consider various measures of the 
‘gain’ to be realized with a convergent tube as compared with a uniform 
tube. One such measure is the ratio G,, = 2,/z,, where z, and 2, are the 
respective values for a uniform and a convergent tube in which the same 
Mach number is achieved, all other things being equal. ‘The simplest way 
of obtaining G,,, other than computing it from (3) and (10), is to read off z 
for the given A,/A,, and for 4,/A, = 1 for a given M,, from a plot such as 
figure 6; the ratio in then computed between these z values. As may be 
readily shown from the analytical expression, the factor Gy, becomes 
infinite, as does z,, when M™, is increased without limit. 
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ADIABATIC INDEX OF DRIVER GAS y, 
igure 4. Limiting value of the equivalence factor g, corresponding to sonic or super- 
sonic cold flow, vs y, for various convergence area ratios in a monotonically 


converging tube. 


A second measure of gain, G., may be defined as the ratio of the shock 


Mach number achieved in a convergent section tube to that in a uniform 


tube for the same diaphragm pressure ratio, all other things being equal. 
Analytic expression of such a factor is not possible, but values of G. are 
readily obtained from plots such as figure 6 by reading the values of M, 
for the given A,/A, and for A,/A, = 1 at a given zs. It is not difficult to 
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show that for ideal gas behaviour G. approaches a limiting value g%~)" 
as increases without limit. ‘hese limiting values are implicit in table 1, 
which will be explained immediately below. Examination of equations (3), 
(8), and (10) indicates that the diaphragm pressure ratio and the cold flow 
Mach number both become infinite at a finite value of u,. Hence, under 
ideal conditions (as represented by equation (2)), there is a maximum 
shock strength attainable in a given tube. ‘lable 1 gives maximum values 
of the shock Mach number (i.e. the values for p,/p, > ©) for various 
values of 4,/A,, both for the temperature ratio 7,/7, = 1 and for 7,/T, = 2. 











































i) 
10° Sere 
EA 
+ sasa5 sees essee 
Soe i #5 3525955858 335255555: 
AH 
| : +H 1232: 
E +4 = rt ' 













































































































































































































































































































































































Figure 5. Curves for computing shock-tube performance. ‘The point at which the 
cold flow becomes sonic is indicated. 

‘Table 1 points out quite strikingly that the Mach number gain resulting 
from the use of a convergent transition section is not particularly large. 
‘The advantage of increasing sound velocity in the driver should be noted. 

Comparisons of performance in convergent and straight tubes have 
been made by Lukasiewicz (1952), Resler, Lin & Kantrowitz (1952), and 
Yoler (1954) (see Alpher & White 1957 for a more detailed discussion 
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Figure 6. Diaphragm pressure ratio vs shock Mach number, for helium driving into 
A plot of shock strength p,/p, vs shock 


argon 


Mach 
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| N,/N2 6:18 6°34 6°47 6°76 
| N,/A 7:48 7:69 7:84 8-20 | 
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H./A | 38-9 40-0 40-7 426 
Table 1. Maximum values of shock Mach number ™, for z Pilpy -> ©. 
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of their work). Lukasiewicz obtains limiting values of the equivalence 
factor g for A,/A, > «, but erroneously restricts them as holding only 
for equal cold flow Mach numbers M, in straight and convergent tubes. 
His work also suggests as a measure of gain the ratio of values of z in straight 
and convergent tubes required to achieve a given value of M;. ‘This does 
not appear to be a particularly useful measure since shocks of different 
strength will result. Resler, Lin & Kantrowitz discuss relative performance 
of the two situations from two viewpoints. First, one may regard the 
equivalence factor g as the ‘gain in effective pressure’ on driving a shock 
by using a steady flow (convergent) rather than an unsteady flow (uniform) 
transition section at the diaphragm. Second, they point out that a tube 
with a convergent section is ‘equivalent’ to a straight tube having an 
initial pressure ratio 


2%, 
(see equation (10)). Their discussion of these matters is quite brief and 
limited to the case of supersonic cold flow. Finally Yoler has used the 
equivalence factor g as though it were the gain factor G,,, which invalidates 
some of the data reduction and conclusions in his work. 


and an initial sound speed ratio (a,/a,)g°"")?" 


EXPERIMENTAL WORK 

Described below are some experiments with a shock tube in which the 
area ratio A,/A, across the diaphragm section was alternatively 1-0 or 1-51. 
‘hey show that the shock enhancement due to a convergent diaphragm 
section is satisfactorily predicted by the preceding analysis. Measured 
quantities were the initial pressures p, and p,, and the shock velocity as 
a function of distance from the diaphragm. ‘The shock tube had a driven 
section 3} in. square in cross-section and 40 ft. long. The driver was 
circular in cross-section, with a 4-5 in. inside diameter, and is 6-5 ft. long. 
Driver gases were helium and hydrogen, while driven gases were air, 
nitrogen, and argon. ‘lhe most extensive data were obtained with 
helium/air, and only these will be reported here since the experiments with 
other gas combinations merely corroborated these results. In all tests the 
driver pressure was slowly raised until the diaphragm ruptured. ‘The 
diaphragms were clamped at the driver end of the transition section, and 
were generally of 0-010 to 0-030 in. thickness stainless steel scored in a cross, 
one-third to a half of the thickness being removed by the scoring tool. 
Driver pressures were in the range of 100 to 800 Ib./in.* and were measured 
with a calibrated Bourdon gauge. Driven section pressures were measured 


with a manometer for p,; > 50 mm of mercury, a Wallace and ‘Tiernan 
absolute pressure indicator for 50 mm > p, ~ 2mm, or a McLeod gauge 
for py < 2mm. 


Since the processes of shock formation and attenuation cause the shock 
velocity to change with distance from the diaphragm, instrumentation was 
devised to enable measurement of the shock velocity as a function of position 
in the tube, and hence to determine the maximum shock strength wherever 
it occurred relative to the diaphragm. ‘The arrival of the shock at an 
instrument station along the tube was detected by a thin-film heat detector 
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of the sputtered platinum type (Rabinowicz, Jessey & Bartsch 1956). 
‘To permit a number of such signals from different stations to be coupled 
to a recording system, it was desirable to lock out signals subsequent to the 
first from each station, and also to make the duration of this first signal very 
short. ‘The heat detector output was amplified and fed to a 6D4 thyratron 
(the more commonly used 2D21 thyratron was found to be unsuitable 
due to its having firing delays of from 8 to 20 microseconds for weaker 
signals). ‘he outputs from six of these thyratrons were fed through 
individual decoupling diodes to a blocking-oscillator type of pulse generator. 
Hence, through two such circuits, brief pulses were generated corre- 
sponding to shock arrivals at up to twelve stations spaced along the tube 
at 32-6 in. intervals. A vertical raster pattern was formed on the screen 
of a cathode-ray oscilloscope by connecting to the y-input a crystal- 
controlled sawtooth voltage generator of 100 microseconds period, the 
v-deflection being actuated by the oscilloscope’s time base with a much 
longer period of about 10 milliseconds. In this manner about 500 cm of 
sweep at a writing speed of 20 microseconds/cm were obtained for a single 
sweep of the oscilloscope’s time base. ‘Time markers at 10 microsecond 
intervals were applied through intensity modulation of the oscilloscope. 
The amplified and shaped (i.e. shortened and freed from extraneous signals 
following the initial one) pulses originating from the heat gauges were 
applied directly to the x-plates of the oscilloscope, and so caused horizontal 
breaks in the vertical raster. ‘The oscilloscope was triggered by the output 
from a barium-titanate transducer, which was fixed in the wall of the shock 
tube near the diaphragm and was sensitive to mechanical vibrations of the 
tube as the diaphragm broke. ‘The shock velocity at different positions in 
the tube could readily be deduced from a photographic record of the raster, 
on which time intervals could be measured with a possible error of about 
a microsecond. 

The shock tube as constructed had an area ratio A,/A4, = 1-51. For 
these experiments an area ratio of unity was obtained by inserting down 
the centre of the driver a wooden dowel which terminated short of the 
diaphragm. ‘This resulted in quite dissimilar cross-sections in the driver 
and driven parts of the tube, but permitted use of the same diaphragms as 
before the modification, so that the mechanics of the rupturing diaphragm 
was not a variable in the comparison. 

In general a rapid acceleration of the shock to a maximum velocity was 
noted relatively near the diaphragm, followed by a gradual deceleration as 
the shock proceeded down the tube. Since the treatment in the present 
paper takes no account of shock formation and attenuation processes, the 
maximum shock velocity observed was judged to be that most suitable for 
comparison with the theoretical results. Since this observed value is actually 
the average over a significant interval, it is in general slightly smaller than 


the actual peak velocity. 
The data are shown in figure 7. Note that the experimental shocks in 
each case are somewhat weaker than the theoretical predictions, but that the 
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shock enhancement due to a convergent transition section is in excellent 
agreement with theory. 

Our experiments have been continued to the study of stronger shocks 
than those described here; and with the present shock tube the indications 
are that, for diaphragm pressure ratios above about 1000 with hydrogen 
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Figure 7. Comparison of ideal shock-tube theory with experiment for a 3} in. square 
shock tube to demonstrate the effect of a convergence at the diaphragm section. 


or helium as the driver gas, the performance of the tube can no longer be 
depicted by an ideal theory. Shocks stronger than predicted by ideal theory 
are observed, and the maximum shock velocity in an experiment can occur 
quite far from the diaphragm. ‘To explain such observations, one must 
consider the finite time for diaphragm rupture, mixing effects in the contact 
zone, shocked gas non-ideality, and wall attenuation effects. Nevertheless, 
even with these stronger shocks, the enhancement due to convergence 
appears to be predicted fairly accurately by ideal theory. 
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SUMMARY 

‘The boundary layer on a hot flat plate which is fixed at zero 
incidence in a slow stream carrying a progressive sound wave is 
investigated. Formulae are obtained for the skin friction and the 
heat transfer in the extreme cases when the frequency is very low 
and very high. In addition, different methods of simplifying the 
boundary layer equations in unsteady compressible flow are 
briefly compared. 

1. INTRODUCTION 

During the last few years there has been a revival of interest in the theory 
of unsteady laminar boundary-layer flows. ‘The solution for the initial 
growth of the boundary layer on a body started from rest has long been 
known, but the recent work has been concerned with the later stages of 
such a motion and with the boundary layers associated with fluctuating 
external flows. One of the first contributions to these new developments 
was made by Moore (1951), who considered the compressible boundary 
layer on a heat-insulated flat plate moving with variable velocity in a uniform 
medium. Moore’s work has since been supplemented by that of Ostrach 
(1955) who treats the same problem for an isothermal flat plate. Both 
these papers are relevant to the boundary layer on a missile whose velocity 
changes continually during flight. 

Another important contribution, dealing with unsteady incompressible 
boundary layers, has been made by Lighthill (1954). He investigates the 
boundary layer in plane flow over a fixed cylinder of arbitrary cross-section 
when the free stream is fluctuating with small amplitude about a steady 
mean value. For incompressible flow, this problem is mathematically the 
same as if the cylinder were moving with the same fluctuating velocity in 
a uniform medium at rest. For this reason, Moore’s and Ostrach’s results 
when interpreted for incompressible fluctuating flow link up with Lighthill's 
predictions concerning a flat plate. 

Lighthill mentioned the Rike tube among the possible fields of 
application of the theory of fluctuating boundary layers, and it was, in 
fact, the phenomenon of the Rijke tube that suggested the work described 
in this paper. Rayleigh (1894) describes how such a tube may be set up. 
He quotes the case of a tube 5 ft. long and 4? in. in diameter, open at both 
ends and held vertically. When a fine wire gauze stretching across the tube 
about 1 ft. from the bottom was made red hot by a flame, and then the 











472 C. R. Illingworth 


flame was suddenly removed, almost immediately a sound of considerable 
power was emitted which lasted for several seconds. In the tube there is a 
stream of air composed of two components, namely the fairly slow convection 
current up the tube combined with the longitudinal motion associated with 
standing sound waves. ‘lhe energy of the sound waves is derived from the 
varying transfer of heat from the gauze to the surrounding air, and this 
changes with the fluctuations in the speed of the air current past the gauze. 
In the Rijke tube, therefore, we have an example of a hot obstacle in a 
fluctuating air stream. It is true that the Reynolds number for the flow past 
the fine wires of the gauze described by Rayleigh is below the range of 
Reynolds numbers for which boundary-layer theory is valid, and therefore 
the theory which follows really applies to an obstacle of larger longitudinal 
dimensions than the gauze. 

The principal aim of the present paper is to investigate the effect of 
high wall temperature on the skin friction and heat transfer on a flat plate 
fixed in a stream carrying sound waves. In addition, we shall take some 
account of the effects of the fluctuations in density and temperature that 
accompany the sound waves in the external stream. (In Lighthill’s theory, 
for incompressible flow, the obstacle may be warm, but not hot, and the 
density and temperature are constant in the external flow.) Our investigation 
is closely related to, but not equivalent to, Moore’s work on a moving plate. 
For if, in Moore’s problem, we change from axes fixed in the medium to 
axes fixed in the plate, so as to give the appearance of a plate at rest in a 
moving medium, we thereby introduce an apparent pressure gradient 
—p dU/dt in the boundary layer, where p is the density and U(t) is the 
velocity of the plate. On the other hand, the pressure gradient for a fixed 
plate in a stream U(t) which fluctuates slightly and has a uniform mean 
density p, is —p,,dU/dt. ‘Vhe two pressure gradients are different unless 
p=p,. Thus it is only for incompressible fluids that Moore’s and 
Ostrach’s results provide the answers to our present problem. 


2. EQUATIONS OF MOTION 
We shall consider the boundary layer on a cylindrical obstacle fixed 

in an unsteady stream. An example of such a stream occurs when sound 
waves are propagated downstream in steady flow past the obstacle, and this 
particular flow when the obstacle is a flat plate will be investigated in some 
detail in §3. ‘The flow is supposed to be perpendicular to the axis of the 
fixed cylinder, and so only two space coordinates are involved, x measured 
along the surface of the cylinder from the leading edge, and y measured 
perpendicular to the surface. ‘he boundary-layer equations then are 

dp (pu)  A(pz) 





ideal gas oe (1) 
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where uw and wv are the velocity components in the x and y directions, 
p, I and w are the density, specific enthalpy and viscosity respectively, 
and P is the Prandtl number, supposed to be constant. ‘The suffix 1 refers 
to the external stream at the edge of the boundary layer and the external 
velocity u,, density p,, pressure p, and enthalpy /,, which are all functions 
of ¢t and x in general, obey the equations 


Op, O(p; u,) 


RE —— = ( 
ae (4) 
Ou, Ou,\ Oy (5) 
PI\ Oi 1 Ox Ox’ 
Bie me 
S144, 21)-(Ft +4, 2) = 0. 
a(S; ia =) (5 My =) ’ (6) 


The last equation expresses the fact that the entropy of each fluid elemen 
in the external stream is conserved, but in the example considered in $3 
we shall use the even stronger condition that the external stream is 
homentropic. ‘To these equations must be added the equation of state 


yo 
ti: Tegel 
p 3 p (7) 


valid for a gas with constant specific heats, and the equation 
pl = pl, (8) 
which, in conjunction with (7), expresses the fact that the pressure p does 
not change across the boundary layer at any station. 
It is convenient now to follow Moore’s method of analysis. This is an 
extension of the Howarth (1948) transformation to unsteady boundary 
layers. First, the y-coordinate is replaced by Y, where 


1 sy 
¥(t,x,9) = — | ~ plt, 38) ds (9) 


in which the suffix oo refers to some standard state of the fluid. Later we 
shall identify this state with the mean state of the free stream. Next, Moore 
introduces a function y satisfying 

Pb a 


u a = 730 10 
p oy 7 (10) 
which closely corresponds to the stream function in steady plane flow. 
(For example, mass flux is measured by p,#s.) Then, the equation of 
continuity (1) is satisfied if 


p, [op oY 
c= - 2 (245) (11) 


and, with ¢, x, Y as the independent variables, (2) and (3) are immediately 
transformed into 


) os o dw 2 \ dub m l Op, I l o O7ubs (12) 
uaa” marley awh * eee 
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‘These equations may be simplified further by making the special 
assumption that the viscosity varies as the temperature 7’ across each station 
in the boundary layer. According to this, 


be r= Kw T i om (14) 
where the viscosity at the wall ,, may be related to the wall temperature 77, 
by Sutherland’s law 
as T\3? T+ Ts 
fo Sat Whar 


where 7’, is a characteristic temperature of the fluid, 114° K for air. The 
formula (14) which was discussed by Chapman & Rubesin (1949), is a 
good approximation to the true viscosity-temperature relation (Sutherland’s 
law) near the wall, all along the boundary layer. Of course it is less accurate 
in the outer part of the boundary layer, where it would be better to use 
yu = p, 7/7,, but errors in the viscosity and conductivity there are mitigated 
by the smallness of the velocity and temperature derivatives. Equations (14) 
and (15) together imply that 


po Cp, I, (16) 

where T,.+T. 7 

C = y!2? 3 17 

xX xT, +7," ( ) 

in which xy = T,,/7.,.. It now follows, with the help of (7) and (8), that 

LP = L*Pos (18) 

where silk ae ses (19) 
Ps 


When pp is replaced by p:*p,, according to (18), the boundary layer equations 
take the simplified form 
2.2 AVY 1B oh 
oO oY ox ox oY] 3} Py ont, eee) 


3 


4 
c 





=, (20) 


‘These are simultaneous differential equations for the mass-flux function 7: 
and the enthalpy /, which are the two quantities we have to determine. 


5 


3. FLAT PLATE FIXED IN A LOW SPEED FLUCTUATING STREAM 

We have already mentioned the flow in the Rijke tube. In order to 
gain some information about the behaviour’ of the boundary layer on the 
obstacle in this kind of flow, we shall examine the boundary layer on a hot 
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tlat plate which is fixed at zero incidence in a low speed stream containing 
a progressive sound wave. We choose a flat plate as the obstacle rather 
than a circular cylinder, which would be more representative of the wires 
in the gauze, because the external flow past it is simpler. A flat plate at 
zero incidence does not, in a first approximation to the boundary-layer 
problem, affect the incident stream. With a circular cylinder, unless we 
restricted attention to wavelengths of sound much greater than the radius 
of the cylinder, we should have to include the effects of the scattering of 
the sound wave in the velocity external to the boundary layer. (Of course, 
in the Ryke tube, the sound waves are in fact much longer than the 
dimensions of the obstacle.) With a flat plate there is no scattering 
whatever the frequency of the sound waves. 

If the stream past the plate has a mean velocity U,, and if there is a 
sound wave of frequency w moving downstream in it, the resultant velocity 
may be written, in complex form, as 


; a: 
u,(t, x) = ¢ J ‘ exp( io > aa i) | (22) 
L M4 


in which « < 1. Here, M is the Mach number of the mean flow and s is 
the frequency parameter wx/U,. ‘The corresponding pressure and 


temperature are given by 


WV \ 
pit, x) = p I ty Me exp( ion a4 is) | (23) 
‘i ‘at ; uv .\ 
Tilt, x) l ; (y —1)Me exp( ion Wo is) | (24) 


‘These expressions are valid for all values of Mand w; the only approximation 
that they involve is that e? and higher powers of « should be negligible. 

We shall confine attention to the case when WV is small, partly because 
this is appropriate to the low speed stream in the Rijke tube and partly 
because the solution of the boundary layer equations, especially for the 
temperature, is thereby shortened. In the extreme case M = 0, which 
applies to an incompressible external stream, the speed of sound then being 
infinite, (22) becomes 


u,(t, x) = U,,(1+ «e™), (25) 
and the corresponding pressure, density, and temperature, are constants 
with the values p,, p, and 7’, respectively. ‘hus, by setting 1/ = 0, the 


progressive nature of the sound wave is obliterated, and the problem is 
reduced to that of the compressible boundary layer on a flat plate in an 
incompressible stream whose velocity is fluctuating in magnitude. ‘This 
problem for a fluid that is incompressible in the boundary layer as well as 
in the external stream has been discussed by Lighthill (1954) for a cylindrical 
obstacle of arbitrary cross-section, and by Gibellato (1955) for a flat piate. 
By taking .V/ = 0 in the following theory we shall be able to show how large 


temperature differences modify their results. 
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‘The simplification to the boundary-layer solution that is possible 
when M is small is brought about by neglecting the dissipation in the 
energy equation (21). ‘This is a crucial step that sets a limit to the magnitude 
of M. For the discarded dissipation is O(M?), while the terms retained 
are O(yx—1), where x = T,,/T.,, and we shall assume that x is O(1). For 
a red-hot wire, y would be in the region of 2. We are therefore committed 
to neglecting any term whose weight is O(M?) compared with those 
retained. 

It follows that the external stream may be expressed as 


u(t, x) = U,(1+ 2), (26) 
pilt,x) = p.(1+yMB), (27) 
T,(t,x) = Toll +(y—1)ME], (28) 


where E = eexp(iwt— Mis). ‘This stream, like the one from which it has 
been derived, is homentropic, and so p,dJ, = dp,. Consequently, with 
dissipation neglected, the energy equation (21) takes the simpler form 
0 ob 0 Opa\I Crunp, & I 
(a+ apa ~ Seav)i,~ P pear(7): 
In solving equations (20) and (29) it will be convenient to separate the 
cases of low and high frequency. ‘The frequency parameter s is equal 
to kx/M, where k is the wave number. In the Rijke tube the quantity kx, 
being of the order of the ratio of the diameter of the wire to the length 
of the tube, is small, but since we are considering / to be also small it 
is quite possible for s to be large. More precisely, we shall examine the 
cases of small and large values of s in turn. 


Case of low frequency 
When s is small, the external stream is given by 


u,(t, x) = U,,[1+(1— Mis)ee”], (30) 
pi(t, *) = poll +yMee'], (31) 
T(t, x) = Ta[1 +(y —1)Mee), (32) 


when terms involving M? are neglected. Care is required in evaluating 
terms like (—1/p,)ép,/ox in equation (20). ‘This is equal to 





Ou, ou, U2 a ‘ ne 
an “ae [(1 — M)is — M(is)*]ee'’” 


according to (26), and does not vanish as (31) would suggest. The 
appropriate form of the function ¢ in the solution for small s is 


¥ = (CU, 7,2) A F(n) +e" y (yf. (33) 
n= 
where 7 = [U,,/(Cv., x)]"?Y. It follows that 
Os 


Hy =u = Uo Pete ¥ Cyfcen | (34) 
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ays a ee 
— = = 4(CU,,v,|x)¥ x 


x | nF) = Fla) +e Zz (s)M(nfeen)— n+ Ifa} | (35) 


n 0 
where the dashes denote differentiation with respect to 7. ‘lhe temperature 
may be written in the form 


@ 


T/T, = x—(x- 1G) +«e"" & (is)"g,(n)- (36) 

The boundary conditions require that _ 

F(0) = F'(0) = 0, F’(o) = 1; 

G(0) = 0, G(o«) = 1; 

fo(0) = fo(0) = 9, fo(%) = 1; 

F,(0) = fi(0) = 9, fi(%) = —M; : (37) 

f,(0) = f,(0) = 9, f,(%0) =0 (nu 2 2); 

£(9) = —(y—-1)My, So ©) = 0; | 

g, (0) = 0, gloj=0 dS). 


When the expressions (34), (35) and (36) are substituted into (20) we 
obtain, from the term independent of «, Blasius’s equation 
F"+3FF’ =0. (38) 
The Blasius function F is tabulated for instance by Schlichting (1955), 
and the numerical values required in the present problem were taken 
from that source. Considering next the terms involving « (€? is of course 
neglected), we obtain from the terms independent of s 


f+ hf + LF 'fy = —yMF" = WyMFF’, (39) 
and, from the terms with the factor s, 
Dy fy = fo-(1— M)[x—(x- 1G], (40) 
where D, denotes the operator 
da? a2 d 
=. UY 2 eee ey gael 3 fr” 
dj +38 oR I = + 3F". 
Similarly, equation (29) yields the equations 
1 
— G'+ihFC' = ( 
Pp G"+3FC ), (41) 
1 ”" yf v yM vit \ fv ' ’ 
poo t FB, = (x | Aer +6 |- (x-1)G(fo-yMF), (42) 
and Dz £8, = Ax-1NGA +o (43) 
where D, denotes the operator 
le d 
eae, Seon ‘eG! 7 
P dy © af dy - 
The required solution of (41) is 
Eom © [F’(s)]” ds / | [F’(s)]? ds, (44) 
) / ) 
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and this tunction is tabulated in table | tor P = 0-72, which has been taken 
as a representative value of the Prandtl number for air. Now, equations (39) 
and (42) would still survive if w were zero, when the external velocity 
would be U (1+ €). In fact, fg and gy represent the perturbations in the 
steady flow solution when the parameters U, and v,, that occur in the 
steady solution are changed to U,(1+e) and v,(1+yMe) respectively. 
‘his last expression arises from the fact that the factor Cv,, p,/p. in (20) 
and (29), which is Cv, when the flow is steady, takes the value Cv,,(1 + yMe) 


in the quasi-steady case (w = 0). ‘Thus 


whence 


f (eZ mh ) Sy M(F yt ) (45) 
Similarly, 
¥ \ 1). Wy(1 (5) l any Vi =~ lox LX 1)G(y)} 
a 
( Ll) My(1-—G) »(X 1)(1 yM)nG’, (46) 
which the first term arises because gy and y—(y—1)G have to satisty 
different boundary conditions. 
he solution of (40) may be written as 
i / f 
1 
] hii Kh a | y (4 Fy) ,(1 [ if ry)ot 


iMat(x—-1+ivfiet+ is | (+7) 


D, fy = G—-(ynF’ + F’) 48) 
D, Ihe G—1l, (49) 
Dy fry = F?—1, (50) 


subject to tie boundary conditions /,,(V) =/,,(U) =f,,(0©) = 0 tor 
n= 1,2,3. hese three differential equations are solved by first finding the 
solution of the homogeneous differential equation 

Dy, fay = B, 


subject to fy)(0) = f,,.(0) = 0, fy,(0) = 1. ‘Then, to find f,, for example, 
the particular integral f,,, of (48) that satisfies f,,(0) = f;,(0) =f7,(0) = 0 


is next determined, and the required function is given by 


fa 7 he (51) 


The functions f,,, (vn = 0, 1, 2, 3) are tabulated in table 1. 
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Finally, the solution of (43) may be written as 

$1 = (x— VD (811 — X82) + Mx(y — gas + (x— Ita — ay nO" — yen t 

+ (x—1+dy)812+ 81 — ay )eua— 807/F'"(O))gis$], (52) 


where £44) £12) £13) £14 and g,, are the solutions of 


Ds gn = —33futrn)G', (53) 
Ds £12 = - 3fi2G’, (54) 
Ds £13 = G- i, (55) 
Ds 814 = —fiaG’; (56) 
Dg HPS, (57) 


subject to the boundary conditions g,,(0) = g,,(«©) = 0 for n = 1, 2, 3, 4, 5. 
The solutions of these are found, as for the corresponding f equations, by 


first determining the solution of 
Dz £19 = 0 (58) 


subject to g19(0) = 0, gi9(0) = 1, and then combining it in turn with the 
particular integrals of (53)—(57) for which g,,,(0) = g;,(0) = 0. The functions 
g,, (n = 0, 1, ..., 5) so determined are tabulated in table 2. 

The skin friction is 


sPis U 12 /ou ee) A hy one Bray) 
T(t, x) = Pure 2 ,,, U2 (= 2. 
ACE: v) Ps (= :) ( “) f f (; ; 4 p 


and if (7,,), denotes the skin friction for steady flow, viz. 








it follows that 





— = 1 + €e"[ 1-5 + (3-661, — 1-105)rs 4 
M'\0-7 — (3-661 + 1-228)2s}], (60) 


1 


for y = 1-4, provided that s is small enough for s? and higher powers to be 
negligible. 
Similarly, the heat transfer from the plate to the air is 


PwPwfl U. \"? fol 
Iu, v) = “Pp, (= -) (=) 


1 (CU nva\% p 7 ae 
pe: —~*] Pt a 1)G'(0) ~ ee” > dye. f (OD 


and if (q,,), 1s the heat transfer in steady flow, viz. 


1 (CU, v,\" 
(Qw)s = 0-296 P Pa (<-" ) (i, i, )s 





x 
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it follows that 


cm} = 1+ee | os + (0-840 —0-910)is + 
Vw s 


M 
+ = {0-7x— 1-1 - (0-840y?— 0-080 +0-198)is} | (62) 
y= 
for y = 1-4, again provided that s* and higher powers are negligible. 


Case of high frequency 

When the frequency parameter s is large, it is appropriate to expand 
in inverse powers of s. If « = (7s)-!* and B = (iw/Cv,,)"”Y, it follows that 
7» = «8, and we may write, instead of (33), 





joe €f, (S2)"|5 F(aB)+E 4 ah,(B) | (63) 
‘hen, - 
= = ¢ = U.| F F’(apB) + EY « ah’ (8) |, (64) 
and 


a) 2 
ae V(iwCv,,)'"| «2B F’(«aB)—aF(«#B)+E no2+2M)«"h,(B) |. (65 
- ) | Af 
Corresponding to (36) we write 


T/T, = x—(x—1)G(aB)+£ ye a"k,,(B). (66) 
For small «, 
F(aB) = }F"(0)a?B? + O(a5), 

and G(«B) = G'(0)aB + O(2*), 
provided f is not too large. With the help of these approximations, we 
obtain, by substituting in (20) and considering successive powers of «, the 
following differential equations 

hy —h, = —x(1-™), 

hi, —h, = (x-1)(1— M)G'(0)8 + MF'"(0)(Ay — Bh}), 

h} —h, = MF"(0)(h, — Bhi), 

hy —h, = \ F’(0)(4 82h), — Bh.) + MF'"(0)(h. — Bha). 
In the complementary function of each differential equation, the term 
involving e® must be suppressed because it has the wrong behaviour for 
large 8, and therefore there are only two constants of integration. ‘These 
are determined by the boundary conditions /,(0) = h)(0)=0. The 
appropriate solutions are 


hy = (1— M)x(B +e? — 1), 


h, = —\Mx—1)G'(0)8? + M[4(y — 1)G"(0) B? 4 
+ yF"(0) B+ teF(9 + 5B + B*) — 33], 
h, = —M(y-1)F'( argo 1 B34 e-B— 1), 


hs = +(1—M)xF ae B? + eP(13 +138 +58? + 28%) — 13]. 
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Similarly by substituting in (29) we obtain the differential equations 
P %, —-ky = 9, 
PR’ — ky = — MLF"(0)Bhy + (x—1)G'O)hl, 





P-'k, —k, = — M[F’(0)BR, +(x — 1)G’(0)Ay], 
Pk; —k, = 4F"(0)B?k, + 4(x—1)G'(0) Bhi, — M[F’(0)BRo + (x — 1)G’(O)hy]. 
| The appropriate solutions, for moderate f, are 
hy = —M(y—1)xe-#¥?, 
k, = My(x- 1)G'(0)| B + =i tet - 1], 
) 


a 
ro 
II 


1 
MOx— 1G P| B(eFY?-1)-1" |, 





2P Ee = 
ks = x(x—-1)G'(0) | —B+ (1— Pp e-BvF ati —?p B+ (1- = he |- 


\ 


{2 


= 1 ree wv ae ee —B il ey 
mM AxQc-D60){ -8+ q—ppeP? - (BAH 
aE 


| + pale th + tly xP OMB + PIER gP BME AY? |, 


the term e®¥? in the complementary functions having been discarded in 
each case, and the solutions having been made to satisfy the boundary 
conditions 

k,(0) = —M(y—-1)x, k,(0)=0 (n > 1). 


It now follows that for large frequencies (« small) 


bw Pw iw \12 ou 
Tull #) = Px ie) (), 0 
72 Cr. 12 Dy wy . = — 
= p,, U2, (7 :) a (O)+E 2 


x 




















grt m0) |, 


and so 





Tw = 


Gy, = 1+ BLS-O1 x(i5)!* — 0-8903(y — 1) + 0-3125 y(i8)* + 





+ M{—3-011y(is)!? + 1-140 + 0-5097 — 0:2956(y — 1)(és) 42 
~()-3125y(is)}], (67) 


provided the cube and higher powers of s-'? are neglected. 
Similarly the heat transfer from the plate is 





a MwPw ta “en él 
gAt, x) = Pp, (=) (FR). 
1 CU,, v,.\12 0 
= peo (— “2) Pi 7 | (x -1)G'(0)—E > a (0) | 





Xv 


Pu 
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and so 
GV | eae M dae sicidaa ae 
ee LE] (3945 y(7s)-! 4 f__ 1-148 y(is)2 — 0: 2 
o [ er + a rn 


+ 2:3410y — 1-8 + 0°3484(y — 1)?(is)-!? — (0-3945x? — 
—0-4226y)(is)1} | (68) 


for P = 0-72, with the cube and higher powers of s~!? neglected. 
If, in order to compare the formulae for small and large frequency, we 


write 
ha re , 
Z.* 1 + ce) 4 + 7B), 
and 
fre _ 1 4 eeilot-0C 4 4D), 
(Vw) 
where 4, B, C and D are real, then 
1 = 1:5+0-7M (small s) | 
2:129s12 — 0-8903(y — 1) + M[—2-129s!? + 1-1403x 4 - (69) 
+ 0-5097 — 0:2090(. — 1)s-!?] (large s), 
B = (3-661, — 1-105)s — M(3-661y — 0-272)s (small s) 
2:1295!2 — (:3125s-! + M[ — 2-129s1? + - (70) 


+ O-2090(y — 1)s-2* + 0-3125s—4] (large s), | 








0-7 — 1-1 
C=054+M- cs Stl (small s) 
x-1 
" (71) 
0-81 l / X o = ‘ 9 
i| = = we st2 4. 1-8 — 0-5410y + 0-2464(y — 1)s "2 | darge s), 
y 
V 
D = (0-840, — 0-910)s— > (0-840? — 0-080 +0°198)s__ (small s) 
0-8117y 
0-3945 ys"! 4 wl ——" sli2 — ()-2464(y — 1)s-1* 4 (72) 
= 
()-3945y? — 0-4225x 
elise! Epica, s ‘| (large s). 
—— | 


‘The angles tan-'(B/A) and tan-'(D/C) are the amounts by which the phases 
of the skin friction and the heat transfer respectively are in advance of the 
mainstream fluctuation. 

In the important special case M = 0, 


7 = 1+ ee'™[1-5 + 13-661 x ~ 1-105)s] iii 
1 + ce'™[2-129512 — 0-8903(y — 1) 4 


1(2-129s!'? — 0-3125s—')] (large s). 


When further the flow is completely incompressible, y = 1, as we have 
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already mentioned, and we then obtain 





— = 1+ee"(1-5 + 2-556) (small s) 


(Tw)s (74) 


1 + ce!"{2 129512 + i(2-129s!2 —0-3125s-1)] (large 5). 








) (We may notice here that Moore (1951) gives 2-555 as the coefficient of is 
in the formula for small s.) ‘The formulae given by Lighthill for the skin 
e friction in this case are 
— = l+ee'(1:54 31 7s) (small s) 
(Tw)e . (75) 


= 1+ee[2-062s!? + 7(2-062s!?)] (large s). 
The good agreement between (74) and (75), and particularly the fact that 
the coefficients of zs for low frequencies are practically equal, strongly 
supports the use of the Pohlhausen approximate method by Lighthill in 
this problem. Lighthill also showed that the whole range of values of the 
frequency might reasonably well be covered by using (75) only. For at 








Figure 1. Variation of 4, B, and tan-'(B/A) with frequency parameter s when 


VJ 0 and y 1. 


the place s = sy where the values of tan-'(/A) are equal in the two parts 


of (75), the values of (42+ B?)!2 agree, and so it is reasonable to assume 


that the two parts of (75) apply to the respective ranges s s) and 
§ Sy). ‘This point is illustrated in figure 1, where the quantities 4, 6 and 
tan-1(B/A) as given by (74), rather than (75), are plott d for various values 


of s. For if the term involving s ! in (74) 1s discarded, so as to bring the 


formula into line with (75), the B-curve for large s will coincide with the 


given A-curve for large s. ‘Then, A B for both large and small s at 


=~ - ~ 
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s = sy = 0-59 and here the two A-curves are close together. Figure 2, 
based on (73), shows the same information for y = 2 as figure 1 does for 
x = 1. However, for y = 2, the whole range of values of s does not appear 
to be as well represented as it was for y = 1 simply by the formulae for small 
and large s. ‘lo make a strict comparison with the previous discussion, 
we should retain only the terms involving s'? in the part of (73) referring 
to large s. Then, the A- and B-curves for large s would coincide with the 
A-curve for large s in figure 1, and we should have A = B for both small 
and large s at s = s, = 0:24. Here the values of A would be 1:5 for small s 
and 1-043 for large s (from the A-curve on figure 1). This involves a much 
larger discrepancy than the case y = 1, where the corresponding numbers 
are 1-5 and 1-64 respectively. In fact, figure 2 suggests that the formulae (73) 
will not be adequate to cover all values of the frequency for an arbitrary 
value of x. 
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Figure 2. Variation of A, B, and tan~\(B/A) with frequency parameter s when 
M 0 and x 2 


The quantity tan''(B/A) is of course the phase advance of the skin 
friction relative to the main stream. As (69) and (70) show, this has the 
asymptotic value }z for large frequencies whatever the values of y and M 
(small). For any specified small value of s, on the other hand, the phase 
advance is increased by increasing the value of y, that is, by increasing the 
temperature of the plate and therefore reducing the inertia of the gas near 
the surface. 

As regards the heat transfer from the wall, it can be seen from (71) and (72) 
that the phase advance for small frequency is small and negative in the 
incompressible case (y = 1) as Lighthill found. In fact, for the case M = 0, 
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D y = 1, Lighthill’s work gives 
: SY = I + ee'(0-5 —0-037s) for small s, 
r (Vwis 
l whereas the present calculations give the same formula with the coefficient 
; 0-03 replaced by 0-070. (Ostrach’s (1955) value is 0-069 to the present order 
g of accuracy.) ‘ihe Pohlhausen method has in fact predicted this coefficient 
C reasonably well, giving the correct sign and a small quantity of the correct 
order of magnitude. As y is increased, the phase advance increases (becoming 
5 positive at x = 1-08), as it does for the skin friction. Since, for large values 
L of s, C and D are 0 in the case M = 0) provided that terms smaller than s'/* 
are neglected, the heat transfer does not lend itself to a graphical 


| representation corresponding to figures 1 and 2. 


APPENDIX 





A note on transformations of the compressible boundary-layer equations for 
unsteady flow 

The analysis of §2 shows that the Howarth transformation as extended 
by Moore provides a very useful method of treating the boundary-layer 
equations for unsteady compressible flow. Now, in the theory of steady 
boundary layers, besides the Howarth transformation there are the von Mises 
and the Crocco transformations. All these transformations hinge upon the 
replacement of.the coordinate y perpendicular to the wall by a convenient 


alternative variable. ‘lhe Howarth transformation uses — | p dy instead 
Polo — 


of vy, as we have already seen; the von Mises transformation uses the 
stream function, which exists in steady plane flow, andthe Croccotransforma- 
tion uses the velocity component u. ‘These have all been described by 
Howarth (1953). ‘The Crocco transformation has been mainly employed 
in studying the steady boundary layer on a flat plate, whilst the other two 
transformations, besides their application to the flat plate, have both been 
used to show that a compressible boundary layer with a non-zero pressure 
gradient in certain circumstances may be given in terms of the solution of 
an associated incompressible boundary-layer flow. It is therefore of some 





interest to examine whether the von Mises and the Crocco transformations 
| can also be conveniently extended to the case of unsteady boundary-layer 
How. 


The von Mises equations 
In the von Mises transformation for unsteady flow, the new variables 
are t, x and &, where ¢& is the mass flux function defined by the equation 


pus abs 


Px OY 
In terms of the new variables 
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oy\-1 voy Op = pa Op oy a 
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ptt = Po (3) and p }, (5 a i) dis + pu ae (76) 
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It is easy to show that the momentum equation (2) becomes 
Ou ou ¢* 1 ou ou 1 op, 1 uo Cu = 
—+u~ | =a dbt+u - -->~ +> =(ppu = }, (77) 
oft Ob | 9 u* ot ON Pp, ON ip p=, Cus Cys, 


in which the first two terms on the ain side may be replaced by 


ra / ( Ou 
u S(u |" 2 dis 
ods Lhe 


and that the energy equation (3) becomes 


= dh) 
ol oly’ 1 ol 1 fo l 
a tt x | = dip tu = (4 u 2) 
Ol ow Jo UC Ox Ot 1, 


| 3 ol pur (3) 2 78 
Pp? u ou pu =) + PS om (/i ) 


Of course these von Mises equations have little practical value if they are 
in all circumstances more difficult to solve than the corresponding 
equations (12) and (13) that result from the Howarth transformation. It 
is scarcely possible to assert that they will be more difficult without a detailed 
study of some special cases, but it will perhaps be worth while to notice 
how unsteadiness of the flow increases the difficulties of solution in one 
or two simple cases of steady flow that have been investigated by means 

f the von Mises transformation. ‘To take an easy example, consider steady 
flow past a flat plate for which the momentum equation in the von Mises 


eu 1 o Ou - 
— —=~ | upu — }. (79 
OX o Ous ( f 3) 7 


This equation was solved by von Karman & Tsien (1938) under the 


form 1s 


assumptions that px = w,(7/T,)" with n = 0-76, P = 1, and that there is 
no heat transfer at the plate. ‘lhe last two conditions ensure that 
f+hwt=l 'U~ , where the suthx refers to the uniform main stream, 
and the first condition shows that wp = pn, p,(I/I,)” ', so that pp is a 


function of w only, given by 


[ (Uj2 : n—I 
[Lp Vea E (1 a) | : (dS0) 
! a UT 


‘Thus, (79) reduces to a second order differential equation for the single 
dependent variable vu. If the expression (80) is substituted in equation (12), 
a third order equation for zy results, and so there is something to be said 
for preferring the von Mises momentum equation in this problem. 

In the field of unsteady flow a relatively simple problem of the same type 
concerns a flat plate in a constant stream when the temperature of the plate 
is a function of time. ‘hen there are no pressure gradients and the 


von Mises momentum equation reduces to 


c Lou ,, Cu ‘ Cu - 
“are (oa ae ene ais — Swe pu = }. (SL) 
ow ,u* ct OX p~ow \ ( ub 


Even if we assume that P = 1, it is no longer possible to write 


T+ he? = 1,4+1U' 
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now that there is heat transfer at the plate, and so yp cannot be expressed 
in terms of u only, as in equation (80). 

However, if, instead of assuming that « < 7", we adopt the approxi- 
mation referred to in §2 and write 

up = Crs p 

where C, given by (17), is a known function of ¢ only, equation (81) is a 
differential equation involving only the one dependent variable u. ‘Uhus, 
by altering the assumed viscosity temperature law, we have kept to a problem 
of the Karman~T'sien type, namely, the solution of a second-order partial 
differential equation for u, as tar as the momentum equation is concerned. 
Of course it still remains to solve the energy equation (78), for information 
about the temperature (and density) in the boundary layer. 

Actually, if we use the rather formal viscosity-temperature law p oc T, 
so that yp = 4, p,, this problem can be simplified still further. For the 


velocity u then loses its explicit dependence on ¢, and the momentum 


ou oO Co u \ 
= ~ Vs Fe Ma = fe 
Ax ous ous ) 


Thus u(x, %) 1s simply the velocity (expressed in von Mises coordinates) 


equation 1s 


} 


for steady incompressible flow past a fat plate. ‘The temperature would 


then be given by solving the energy equation 
ol ol V ( { C 1\ , Ou - 
— +u~ ~u—lu—)+v,u(—)}. 
ct OX Pod ( a) ( i) 


The Crocco equations 


In the Crocco transformation the new variables are ¢, x and x, where z 


is the non-dimensional velocity distribution function given by 
x = u(t, x, v)/u,(t, x) ($2) 


When u, is a constant this is precisely the transformation that Crocco used 








to investigate the steady boundary lay on H plate Follow1 if Crocco, 
we eliminate pz between the transformed equations of continuity and 
momentum, and introduce the shearing stress 
Ou fcy\- 
T b= ijLAU — 
\ COs 
as a dependent variable in place of u. In this wav, (1) and (2) lead to the 
equation 
Ou, ( [Lp \ C (uP ly , < Lp \ tite Oo { L4p \ 
~ os ily = Uy oe bis { } 
i O3\ T GEN FT ON OZ \ T CENT, 
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Similarly, by transforming (2) and (3) and eliminating pv between them 


212 
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we obtain 


ols ou, od ou, , al <n 
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‘The boundary conditions that go with these equations concern the 
values of 7 and J/ at the wall, z = 0, and at the outer edge of the boundary 
layer x = 1. Equation (2) shows that ¢é7/dy = dp,/ex at the wall, and so 


Or Op, ) 
T os = by Wy = at 2 =U. (85) 


Since the shearing stress vanishes at the outer edge of the boundary laver 
7r=0 at z=1. (86) 
We shall suppose that either the wall temperature (enthalpy J,,,) or the heat 
flux, q,., from the wall is specified. It follows that / must satisfy either 
[=I,, or (rol/oz) = —Pu,q, at z=0. (87) 
Finally, the condition 
f=f, atz=1 (88) 
completes the list of four boundary conditions. 

Equations (83) and (84) are quite complicated, especially because of 
the large number of terms on the left-hand side in cach case. In this respect 
they are much worse than either the Howarth—Moore equations (12) and (13) 
or the von Mises equations (77) and (78). However, it would be feasible 
to use them in some circumstances. Some preliminary work showed that 
it would have been practicable to obtain the results of §3 by solving these 
equations, but the simpler calculations involved in starting from the 
Howarth-Moore equations turned out to be preferable. One advantage 
of the Crocco equations is that the momentum equation is of the second 
order, but one of the complications encountered in the work just mentioned 
arose from the frequently occurring term pp/t. It will be recalled that in 
s 
arising in pp/7, is another series with complicated coefficients. Of course 


$3, 7 was expressed as a power series, and the reciprocal of such a series, 


these difficulties would also be present if the Crocco equations were applied 
to a steady boundary layer with a non-uniform main stream. So far as the 
author is aware, even this has not been done. 

‘The Crocco equations are considerably shortened by assuming that the 
external stream is uniform. In fact for the simple problem, mentioned 
above, of a flat plate with a varying temperature and with up = p, p,, 


they would reduce to 
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The solution of (89) for 7 is known from the work of Crocco for the case 
of steady flow, and the problem therefore reduces to solving (90) for J. 

in conclusion, the small amount of evidence from the problem considered 
in this paper suggests that the Howarth—Moore transtormation is the simplest 
to apply of the three available transformations in unsteady compressible 


boundary-layer theory. 
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SUMMARY 

‘The linearized equations of motion show that in a viscous 
heat-conducting compressible medium three modes of fluctua- 
tions exist, each one of which is a familiar type of disturbance. 
‘he vorticity mode occurs in an incompressible turbulent flow, 
the entropy mode is familiar as temperature tiuctuations in low 
speed turbulent heat transter problems, and the sound mode is 
the subject of conventional acoustics. A consistent higher order 
perturbation theory is presented with the only restrictions being 
that the Prandtl number is } and the viscosity and heat conductivity 
are monotonic functions of the temperature alone. ‘The theory is 
based on expansion of the disturbance fields in powers of an 
amplitude parameter «. ‘lhe non-linearity of the full Navier 
Stokes equations can be interpreted as interaction between the three 
basic modes; in order to help physical insight the interactions 
are classed as ‘ mass-like’, ‘ force-like ’, and ‘ heat-like’ effects. 

Besides the amplitude parameter « there is another subsidiary 
non-dimensional parameter « which indicates the relative 
importance of viscosity and heat conduction ettects as compared 
to the inertial effects, € is proportional to the ratio of the molecular 
mean free path and the characteristic length of the flow pattern 
(Knudsen number). ‘Che main contribution of the paper is thc 
outline of a consistent successive approximation for an arbitrary 
order in » and the presentation of explicit formulae for the 
second order (bilateral) interactions. 

A special case of rather general significance is treated in 
more detail. ‘This is when all three basic modes have intensities 
and length scales of the same orders of magnitude and in addition 
to % the parameter « is also small; the second-order interactions are 
then relatively few and easily identifiable and are shown in table 1. 

Ihe present analysis also sheds some light on the ‘ zero order ’ 
approximation which treats the vorticity and entropy disturbances 
as a ‘frozen pattern’ and the sound field as propagating non- 
dissipative waves. ‘lhe interpretation of hot-wire measurements 
relies heavily on these simplified models and the present paper 


lends som«e Support to these current hot-wire practices. 
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Ll. INTRODUCTION 

In fluid mechanics the non-linearity of the governing equations 
manifests itself in several ways. Some of the non-linear effects arc 
described by various special names, for instance, the continuous generation 
of vorticity in the phenomenon called turbulence, the gradual steepening 
of the slope of compression waves, the generation of aerodynamic sound, 
acoustic streaming, scattering of sound waves by turbulence, etc. The 
question arises of whether all types of non-linear phenomena governed by 
the compressible Navier-Stokes equations have been uncovered or whether 
there are still others that have escaped attention. ‘his can be answered 
only by a systematic study of the equations. 

There is another important question; by varying some overall parameter 
(e.g. the Mach number of the mean flow, the length scale of the fluctuation, 
etc.) the relative importance of the various non-linear effects may vary. 
A typical problem is the change in the nature of turbulence as the Mach 
number of the primary flow increases. In actual experimentation the 
measured quantities must be interpreted within a theoretical framework 


even if it serves only to guide our thinking. ‘The measurement of 


fluctuating quantities in a turbulent How has been extended to high speed 
flows (ISovasznay 1953) and the present approach was strongly motivated 
by the desire to generalize the crude but simple concepts that wer« 
obtained from the linearized problem. 

Both at iow speeds and at high speeds, the hot-wire measurements are 
carried out keeping the probe at a fixed location and the records are 
obtained as fluctuations in time and ail statistical averaging is performed 
in the time domain. ‘I’hese time records at low speeds are often regarded 
aS approximate instantaneous space traverses by virtue of ‘Taylor's 
hypothesis that the flow pattern is carried with the main flow and changes 
only relatively slowly. In high speed flow the simple linearized theory 
of the compressible case indicated that there are three fluctuation modes 
and for two of them (vorticity and entropy modes) this simple motion of 
‘frozen pattern’ can be immediately extended. We may call these 
‘parabolic’ modes as they obey parabolic differential equations. ‘lhe 
third mode, the pressure or sound mode, obeys a hyperbolic differential 
equation and can be interpreted as a system of acoustic waves travelling 
with respect to the main flow. 

The linearized theory does not indicate any interaction among the modes 
as long as the domain of interest is far from solid boundaries and it may h« 
valid as long as the fluctuations are weak. When the inten sity of fluctuation 
ceases to be small, the flow can still be conveniently interpreted as built up 
by superposition of the three basic modes of fluctuation. However, in 
addition to those originally accounted for, there are also fluctuations 
produced by the interaction of the basic modes. A systematic study of 


these interactions 1s expected to cover ull the possibl non-linear 


phenomena. 
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The analysis presented here does not attempt to ‘ solve’ compressible 
heat conducting viscous flow problems in the sense of an initial value or 
boundary value problem. It attempts rather to serve as a guide to assess 
all the non-linear interactions in a consistent framework in which 
experimental information (of the usual incomplete nature) can be 
understood. 


2. FUNDAMENTAL EQUATIONS 
The fundamental system of differential equations governing the motion 
of a viscous, heat-conducting, compressible medium consists of the 
continuity, momentum and energy equations. If p, p, 7, F, v, w and « 
denote respectively the pressure, density, temperature, entropy, velocity, 
coefficients of viscosity and heat conduction, the fundamental equations 
assume the following form : 


Op ,. ? 

ay +divpv = m, 2.4) 
p - = — grad p— } grad(u div v) + div[p def v] + pf, (2.2) 
ce = div[« grad 7] + p[3(def v) : (def v)— #(divv)*7]+Q0, (2.3) 


where D/Dt is Stokes material derivative and def v is the rate of deformation 
tensor. (If the components of the velocity v, referred to axes Ox, x, x3, 
are Vj, Us, V3, then def v is the second order tensor whose 7j-th component 
is (dv,/dx,; + 0v,/dx;) and (def v):(defv) denotes the double tensor product 
>> (euv,/0x; + dv;/ex;)?.) In the equations, m denotes the rate of mass 
production (or injection) per unit volume, f denotes the body force per 
unit mass and Q denotes the rate of heat addition per unit volume. 

‘The coefficients of viscosity and heat conduction « are assumed to be 
known monotonic functions of the temperature, i.e. uw = u(T) and « = K(T). 

For a perfect gas, the variables p, p, T and E are related by the equations 
of state 


i 


p = pRT (2.4) 


wa E-E, = Cylog(2)(22) (2.5) 
Pr}\P 


where y= C,,/C,, C, being the specific heat at constant pressure and 


mom 
C, that at constant volume. ‘The subscript r designates some standard 
reference state. 

If m, f, O are known functions of x, and t, the above five equations, 
together with the boundary conditions, are expected to be sufficient to 
determine the five unknowns p, p, 7, E and v. We shall be ultimately 
interested in the case m= 0, f= 0, and O = 0, but in the scheme of 
successive approximations the effect of the lower order terms on the higher 
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ones will be expressed most conveniently through these apparent ‘ external’ 
forcing functions. 

Let us consider a four-dimensional space-time domain G, the size of 
which is determined by several considerations discussed in the ensuing 


paragraphs. Let 
V = | | d&dr 
denote the four-dimensional volume of G where d& stands for the volume 


element dé, dé, d&, and dz for the element of time. We define the mean 
values of the pressure and entropy in G by 


1 < - 
Po =|) P(E. 7) dbdr, (2.6) 
; if « ” 
b= | | E(&,7) d&dr, (2.7) 
JJG 
and define py, Ty, uo, Ky by the equations* 

Po — iy att «, (2.8) 

1 ’ . Po Pr’ 
a ae slog(2)(2) (2.9 
P+} \Po | 
Ho = LT), Ky = K(14). (2.10) 


Finally, the coordinate system is chosen so that 
= T | v(&,7) d&dr = 0, (2.11) 
that is, the coordinate system moves with the mean velocity of the fluid. 

It should be noted that the mean values py, £,, etc., depend upon the size 
of the space-time domain G. Consequently the fluctuations p— py, E— Ky, 
etc., as well as the ‘ intensity ’ of the fluctuations, ( p — po)/Po, (E — £q)/ Eg, ete., 
also depend on the size of G. If we assume that p, FE and all other flow 


Vo 


variables are continuous functions ot the space and time coordinates, then the 
‘intensity’ of the fluctuations in G can be made arbitrarily small by limiting 
our consideration to a sufficiently small space-time domain G. Let « be 
a non-dimensional parameter characterizing the intensity of the disturbance. 
For example, we may take « as the maximum of |p—fp '/'po', |v/ay ete., 
where a, is the velocity of sound of the mean state. 

We now assume that in a sufficiently small space-time domain G, the 
pressure, density, temperature, entropy, p, p, 7 etc., can be expanded as 
power series in x. ‘This is the first condition which limits the extent of the 
domain G. ‘The actual size of G is determined by the values of the spatial 
and temporal gradients of the fluctuations relative to our coordinate system. 

Note that p,, Ty, etc. are not defined in the same manner as /, and /, mainly 


to simplify the algebra. ‘he ultimate physical conclusions are not affected by the 


manner in which p,, etc. are defined 
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For low speed turbulence Lin’s (1953) analysis indicates that the ratio 
of the temporal gradients to the spatial gradients may be small (compared 
with the mean velocity), in fact of the order of the velocity fluctuations 
themselves. ‘Che domain G would then be elongated along the time axis 
as can also be suspected from equations (6.5). For very small fluctuations, 
such as during final period of decay of grid turbulence, G may be quit 
sizeable. 
We write 
PD = pot p+ p™ +..., (2.12 a) 


P = pp t p™ + p® + ..., (2.12 b) 


where P I) p\! are of order . and p>, pare of order x" 
Substituting (2.12) into the fundamental system of equations (2.1) to 


(2.5) inclusive and collecting terms of the same order in x, we find 








op" ‘ i 
5 + pV. vl) = FO, (2.13 a) 
Avie) ; : A er 
Po yO Vp fe [yV?Vv LugV(V. vi" ) Py. (2.13 b) 
», Of" a : - 
a ra Kao I (2 13 c) 
me [" : 
i ae 
Pi Po ly 
kh y To p " pe . ; 
eee" s (2.18 6) 
Ry 0 Po 
The functions F'” for n | represent the non-linearity of the original 


equations. For n l, 


F\!) = m, = pl, Fi) = O, F') = 0), FY) = 0. 


For n > 1, the F® depend only on the p”, p™, etc., where k = n—1. 
For example, F'> V. [pv]. The principal advantage of the scheme 
presented here is that the differential equations are the same for all values 
of n, only the ‘inhomogeneous terms’ are different, and in the scheme of 
successive approximation these are computed from the solutions of the 
lower order equations. 

It is convenient now to eliminate two of the four thermodynamic variables 


and introduce non-dimensional quantities for the remaining two. We write 


(1) 
po =P (2.14) 

YP 

ki 
Ss) = ——., (2.15 
C pactey 

p 
‘The speed of sound of the undisturbed medium ts a5 = ypy py. We assume 
further that the Prandtl number (uC, «) = 3}, the sole purpose of this being 


to simplify the algebra. 





























Interactions in a viscous heat-conducting compressthle gas 499 


The governing equations now become 


a Pom as m” 


V vit) 4 = <n 2.l6a 
. ot ct Po en 
ove" ss = 
a ay VP — vy, V20™ — ly V(V. v) f, (2.16 b) 
aS Ow 
—— — fn V2S — (yp — 1)vyV?2P” = (2.16 c) 


3¥%o } T 7p > 
ot pol, | i 
where v) = /49/p, denotes the kinematic viscosity of the undisturbed medium 
and 
mi m, fo — f, OP =Q. 


The higher order m, £°, O”” are functions of the lower order solutions 


For the particular case n = 2, we find 
m? as bo. P She a 
— : ViCP® — SO?) + = = Uy 1)(P)? + (S®? + POS), (2.17 a) 
Py £40 
ovil 1 
f? (Pa SY) — (vv .V)vi ay] vv) | Tt 
ct " Po 
a 
+V. | det ve | (2.17 b) 
Po 
Or ago 
Ac se yp) ig Figen sv. | & vgs L(y Pe | 
pat p 1 ) : ct Po 
aa [}(def v”) : (def v) 3(V.v'))?] — 
2y I y 
“3 —— nhl P)2 — 2 PV S0 Sea (So |. (2.076) 
« 7 / 


In these equations, «’/p, is a function of P™ and S”, for, by ‘Paylor’s 


e€xpansion, 


du l d? L - P 
u( 7) = nw, +| += T—71,)+ sl +s (T—T,)? +... 
(T) fh (sr), ‘ (7 fips) xi a7), (7 o) 


se o(de) — Pya(dey — OOEW)7. 
~~ 77); , 77), i ee, 


so that pw) = T(dp/dT) pop, wo? = T(du/dT) pop, + (T)?(d2u/dT?) po, 


f 


etc. If we define uw, = (du/dT),_,., then 





Hy ray — Mo, Ty [S® +(y — 1) P}. 


A considerable simplification follows from the fact that the systems of 
equations (2.16) are identical in form for all values of m. ‘The ‘source 
functions’ m, f, QO correspond to the actual rate of fluid injection, 
body force and rate of heat addition, if any, while m®, f°, O® are calculated 
from the solution of the first-order equations and represent only apparent 
rates of mass, momentum, and heat injection. ‘This concept allows us to 
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study the behaviour of the fluid motion governed by (2.1) to (2.5) in two 
steps. Firstly, we treat the system (2.16) as if m”’, £™, OW were known 
and secondly, examine the structure of the quantities m”, f/, O™, 


"THE THREE BASIC MODES OF FLUCTUATIONS AND THEIR PRODUCTION 


Inasufhiciently small space-time domain G the intensity of the fluctuations 
in the flow field can be made arbitrarily small and terms of order 2? and 
higher powers of x in (2.12) may be neglected. We shall take » = 1 in 
equation (2.16), temporarily suppress the superscript (1), and obtain for 
the first-order fluctuations 

oP oS m 


V.v- = Ot = Os? (3.1 a) 

7 zs a‘, ve vyV°V 5 vyV(V. v) = f, (3. l b) 
as O , . 
es A\paveP = =. $c1:¢ 

Ot —3M%yVS— 3(7 1)vyV*P Pe a ZF (3.1 ¢) 


We can describe the kinematic field by two new independent variables, 
the vorticity Q and the specific dilatation rate q (the rate of volume increase 
per unit volume of the fluid) where 


B= ¥ xv, q= V.v. (3.2) 
‘aking the curl and divergence of (3.1b), we have 
dQ 
ry — v,V" Q — V 4 if (3.1 by’ 
aq 9 9 > ” 
or —sV)V? g= —& V2P+V.f. (3.1 b) 


‘The two equations (3.1b)’ and (3.1b)’ are equivalent to (3.1b), since 
once Q and gq are known, the velocity field can be calculated from (3.2) 
and the appropriate boundary conditions. 

By a slight manipulation as indicated below, the system of equations 
(3.1a), (3.1b)’, (3.1b)” and (3.1c) can be further combined into the 
following four equations : 





€ = — yV2?2=Vxf, (3.3 a) 
: V2P > = (V?P) Ls (;; a¥e¥" E V.£- = << ’ 
a® ot* 5 ascot aL \et Po dt\py C, Ty 

(3.3 b) 

os 1 v2 S 1 V2P O ake 
> My (y— 1)r% re sky (3.3 €) 

oP os 

ee ee) 3.34) 


Equation (3.3 a) is the same as (3.1 b); (3.3.c) is the same as (3.1.¢); (3.3 d) 
is obtained by substituting (3.2) into (3.1a); (3.3b) is obtained by 
eliminating S and g from (3.1 b)”, (3.3 d) and (3.3 c) 
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The system of equations (3.2) and (3.3) can be conveniently split into 


three sub-systems as follows. [Let us write 











RQ = Qo+ 2, + Q,, (3.4) 
P = Po+P, +P, (3.4b) 
etc.; then 
dQ. . en 
a —vV22, = Vx, (3.5.a) 
P, = 0, 5S, = 0, q (), (3.5 b) 
VXV9 = Qo, V.Vo = 0; (3.5) 
Q., = 0; (3.6a) 
oP eats oo ¢ , \m é O 
aa i — V2P = 37Vo >; Y"", “ag ghyV" im V+ = sar (3.6b) 
ot ’ ot ot Pry Ot pyC,, le 
os 
7 aa aVpV°S,, 3 ) l \yV7P.,; (3.6 C) 
oP os m 
q, —-—~"*+—-+-, (3.6 d) 
at a 
Vxv,=0, V.V, = Gp (3.6e) 
2. = 0, P 0) CS. /'2) 
os. O 
a = MeN, =a (3.7 b) 
ot - Po ( p F, 
aS, = 
q. = — ’ tice GC) 
V xv, = 0, V.V, = ,. (3.7 d) 


Of course, there are infinitely many ways in which we can split up the 
equations (3.2) and (3.4)*. However, the above choice is a convenient one 
since (3.5), (3.6) and (3.7) each represent individually familiar phenomena. 

Thus (3.5) is identical with the equations describing the production, 
convection and dissipation of weak vorticity fluctuations in a viscous 
incompressible medium and will be called the vorticity mode. Such a flow 
tield does not generate pressure fluctuations of the same order, since the 
latter is proportional to the square of the velocity fluctuation and hence is 
negligible when the vorticity fluctuation is weak. For a similar reason, 
there is no production of entropy fluctuation, since viscous dissipation 1s 
again a second-order quantity. Finally, the velocity field is solenoidal as 
it should be for an incompressible flow. 


It is possible to decompose the sound mode further into two types of physically 
distinct fluid motion: (1) that due to mass sources in an incompressible medium 
which is characterized by P os, = &y-.2 % = wipe GB, = Tx¥, = 9%; 
and (2) that associated with the propagation of pressure waves in a compressible medium 


which is characterized by (3.6 a), (3.6 b), (3.6 c), (3.6 e) and q,, eP,/ ot aS,/ ot, 
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;quation (3.6) is identical with the equations describing the production, 
propagation and absorption of pressure waves in a viscous compressible 
heat-conductive medium, and will be called the sound mode. As a result 
of absorption, there will be entropy changes accompanying the pressure 
Huctuations. ‘These entropy changes are denoted by S,, in (3.6). The 
velocity field is irrotational as is expected for a sound field. 

Finally, (3.7) can be recognized as the equations describing the production, 
convection and diffusion of hot spots in a heat-conducting fluid medium and 
will be called the entropy mode. Such a diffusion phenomenon does not 
generate pressure fluctuations and produces only a weak irrotational velocity 
field, necessary for the conservation of mass while the density varies. 

The vorticity mode is the type of fluid motion most frequently 
encountered in the mechanics of viscous incompressible fluid. ‘The sound 
mode is the type of fluid motion discussed in acoustics and in the theory of 
compressible fluids. ‘The entropy mode has been the main subject of 
investigation in the theory of heat transfer in fluids. We thus see that the 
solution of (3.1) can always be thought of as consisting of the superposition 
of three basic modes of fluctuations that are familiar in branches of fluid 
mechanics. Naturally, not all the three modes of fluctuations will enter 
significantly into a particular problem, for example, in turbulent heat 
transfer problems and high-speed shear flows, the fluctuations consist 
principally of the vorticity and entropy modes. 

We now assume that the space-time domain G does not contain any 
solid boundaries (this is the second condition which limits the size of G). 
In such a case, each of the three modes of fluctuations evolves as though the 
others were absent. Should G contain some solid boundaries, the last 
conclusion is no longer true, since the boundary conditions are imposed 
upon the resultant flow variables, such as the velocity vector vor 
temperature 7 and not separately upon the individual modes (see, for 
example, ‘Trilling (1955)). 

Decomposition of a small disturbance into different modes of fluctuations 
has been made previously by several authors, Carrier & Carlson (1946), 
Lagerstrom, Cole & ‘Trilling (1949), Cole & Wu (i952), Wu (1952), 
Kovasznay (1953). ‘Phe terminology and classification adopted here follows 
closely that of Kovasznay. ‘The modes of fluctuations are re-defined here 
to include the influence of body forces, mass and heat sources. Ihe fluid 
dynamic effects produced by these agents will be the subject of discussion in 
the next section. 

When the space-time domain G is not sufficiently small the intensity 
of fluctuation is not small and terms of the order of second and higher powers 
of x in (2.12) must be retained. Since equation (2.16) has the same form 
for n = 1 as for n 1, the flow field can still be thought of as consisting 
of the three basic modes of fluctuations. Since the apparent source terms 
m”, £, O"”, in (2.16) are functions of the fluctuations PH) S®,) vy where 
k <n, the terms P, S“, v’? may be interpreted as the pressure, entropy 
and velocity fluctuations generated by the nth-order interactions of the three 
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modes of fluctuations. Of all these interactions the most important are the 


bilateral interactions of the three modes which contribute to the values of 
P®, S® and v2, 
he entire analysis is partly motivated by the interest in compressible 
How turbulence and like phenomena. In those cases the fluctuation levels 
in the entire space-time domain are of the order of a few per cent of the 
corresponding mean values. ‘lhe theory described here is valid for 
interactions of any order provided the process converges in G. Its 
application to the bilateral interactions of the various modes of fluctuations 
will be given in §5 and § 6. 


4+, EFFECTS OF MASS ADDITION, BODY FORCES AND HEAT SOURCES 


(1) Effects of mass addition 

From equations (3.5), (3.6) and (3.7), we see that the effect of mass addition 
is to produce the sound mode. Closer examination of (3.6) reveals that the 
effects of mass addition are really two-fold. When the fluid is injected 
into the medium at a point or a region, it displaces the fluid that was 
originally there. Generally, the movement of the displaced fluid in turn 
generates pressure waves which propagate into the surrounding medium. 
However, mass addition with a spatio-temporal pattern obeying the diffusion 
equation does not cause pressure fluctuation, see (3.6b). It induces an 
incompressible source flow, see (3.6d), with a velocity potential ¢,, obeying 
a Poisson-type equation 





(4.1) 


(2) Effects of body force 


Referring to (3.5), (3.6) and (3.7), we conclude that the effects produced 
by body forces are two-fold. ‘They generate the vorticity mode and the 
sound mode. If the force field is decomposed into two component fields, 
an irrotational and solenoidal field, the irrotational force component produces 
the sound mode (cf. (3.6b)) and cannot produce the vorticity mode (cf. 
(3.5a)). On the other hand, the solenoidal component produces the vorticity 
mode and cannot produce the sound mode. Consequently, a harmonic 
force field (which is both irrotational and solenoidal) can produce neither 
vorticity fluctuations nor pressure fluctuations. Its only effect is to produce 
a harmonic velocity field given by (see (3.1 b)) 


—— = fy. (4.2) 


(3) Hffects of heat addition 

From (3.5), (3.6) and (3.7), we see that addition of heat to a medium 
generates the entropy mode and the sound mode. Physically this is clear 
tor heat addition to a gas produces two effects. Firstly, it increases the 
temperature of the gas (and hence the entropy) and, secondly, it causes an 
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expansion of the heated fluid element. The first action is a more localized 
one and is described by the entropy mode. On the other hand, the 
expansion of the volume occupied by the heated fluid produces in the 
surrounding medium pressure waves, much like those generated by the 
inflation of a solid body. ‘The action produced by the expansion of the 
heated fluid is thus described by the sound mode. 

At the end of the preceding section, we showed that the fluctuations 
Pw, S™, QO” can be interpreted as being generated by an apparent mass 
addition m”, body force f and heat sources O\, resulting from the 
nth-order interaction of the various modes of fluctuations. ‘The effects 
produced in the flow field by a particular type of interaction will then 
depend sac the manner in which the particular interaction affects the 
value of m”, f°, When the particular interaction gives rise pre- 
iulanele toa (> Sas to m, we shall describe the effects produced 
as ‘ mass-like etfects’, and we shall speak of the interactions as having a 
‘mass-like action’. Likewise, we can speak of ‘ force-like effects’, ‘ heat-like 
effects’ or ‘ force-like actions’ and ‘ heat-like actions’. ‘Thus, for example, 
the bilateral vorticity—vorticity interaction will be shown to produce mainly 
a force-like effect, and a slight heat-like effect. A systematic study of the 
effects produced by bilateral interactions of the various modes of 


Huctuations is presented in the next section. 


5. EPFECTS OF BILATERAL INTERACTIONS OF THE VARIOUS 
MODES OF FLUCTUATIONS 


General considerations 
In §3 we saw that the first-order system can be considered as the sum 
of three basic modes of fluctuations. ‘Their contributions to the thermo- 


dynamic and kinematic variables are, from (3.4), (3.5), (3.6) and (3.7), 
i = Be’, PR ae FS, SO = SY + SY, (5.1 a) 
(1 7 \ 5 
v v Vi +; (5.1b) 


Substitution of (5.1) into (2.17) leads to the following expressions for 
the second-order source functions in terms of the separate first-order mode 
contributions as indicated by the double subscripts. (For example, fo, 
represents the apparent body force produced by the vorticity-sound 


interaction. ) 


m®) = polMog + My» + Mgg t+ Mg, + Myg+ Mol, (5.2 a) 
a . : : ae 

f foo t+ fy +f + fon + fs tfso (5.2 b) 

Oo Po ( p TI Qoo os ' Q.. . Voy, ' Ons , Qa]. \ 5.2 c) 


(he 18 interaction terms can be computed from equations (2.17) and 
a brief examination leads us immediately to the conclusion that bilateral 
interactions of the three basic modes produce in general fluctuations in all 
three modes. In this case we have all 18 second-order interaction terms as 
indicated by Kovasznay (1953). 
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It is a natural question to ask whether all these 18 terms are equally 
important. Naturally such questions can be answered only if more is stated 
about initial conditions and boundary conditions. In a concrete flow field 
these are imposed by grids, coolers, fans, etc. One mode can often be 
neglected and this immediately reduces the number of interaction terms by 
one-half. 

In the next section an order of magnitude estimate of the non-linear 
interactions will be carried out for an important special example still having 
sufficient generality to illustrate the nature of these second-order interactions. 


6. BILATERAL INTERACTIONS OF THE THREE BASIC MODES OF FLUCTUATIONS 

Let us assume that the particular region of flow is far away from the 
boundaries and all three modes of fluctuations are present. Each mode of 
fluctuations can be pictured as a superposition of a large (or infinite) number 
of Fourier components. ‘Typical components for each of the three modes 
of fluctuations have been discussed by Kovasznay (1953), or can be found 
directly from (3.5), (3.6) and (3.7) with m, fand QO put equal to zero. These 
are summarized below. 


(1) The vorticity mode 
The typical component of the vorticity fluctuation is given by 
QY = Qi) explike.x— vy kot], (6.1 a) 
where Q(.k, = 0 and Q(!) is the complex amplitude fluctuation. kg is 
the wave-number vector, ¢ denotes time and x the position vector. ‘The 
corresponding velocity fluctuation is 


Ko x 926) 


vil) = io ae exp[ik, .x— v9 kot]. (6.1 b) 


(2) The sound mode 
The typical component of the pressure fluctuation is given by 


P® = PY ) exp[7k,, .x- ct], (6.2 a) 


where P%) is the complex amplitude; k,, the wave-number vector, and 


c is the complex number 


2y vyk ( 4y?/v,k,\2) v2 
ecg Ae a £0 es Cee ? 
cm ayk,( —. 14 | 9 ( : ) f ) (6.2 b) 


‘The real part of c gives the rate of damping while the imaginary part gives 
the frequency of oscillation. Note that the absolute magnitude of c 1s of 
the order of ayk,,, since in general (v9 k,/a)) < 1. (See also the paragraph 
following equation (6.4).) The entropy generated as a result of the 
dissipation of the pressure waves is given by 


a(y -] Vo ke 


p? 


Ss,” 


y 


P\) exp[ik, .x— ct]. (6.2 c) 





eaeay. Sp 
C— gk, 


2K 
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The velocity fluctuation associated with the sound mode is 


2 R 
_ « os Pi’ exp[tk,, .x — ct]. (6.2 


p —. 
C— uk, 


(3) The entropy mode 
The typical component of the entropy fluctuation associated with the 
entropy mode is 


SY = SO expfik,.x—4y,k?¢], (6.3 a) 


where S‘') is the complex amplitude and k, is the wave-number vector. 


8U 


The corresponding velocity fluctuation is 


vil) = ifvyk, S\) expltk, .x— vk? i]. (6.3 b) 


(4) Relative orders of magnitude 

We are now ready to estimate the relative importance of the various 
interaction terms. By inspecting (5.3), (5. 4) and (5.5) from the point of 
view of substitution of periodic solutions given by (6.1), (6.2) and (6.3), 
we can assert the following facts: 

(a) Fluctuations in each mode may be characterized by two quantities, 
namely, an amplitude and a length scale. 

(6) ‘The medium has an intrinsic length scale, namely vp/ap. 

(c) In the interaction formulae, differentiation with respect to space 
coordinates does occur but integration does not; consequently, the 
terms containing derivatives decrease monotonically with decreasing 
wave-numbers (with increasing length scales). 

We shall discuss the more complicated case where all three fluctuation 
modes have amplitude and scales, respectively, of the same order of 
magnitude. If one of the amplitudes or one of the wave-numbers is smaller 
than the other two, this can only simplify matters and never complicate 
them. 

Whenever a quantity is differentiated with respect to the space coordinates, 
factors of k and k?, etc., will appear. We shall assume that all three 
characteristic wave-numbers are of the same order of magnitude, i.e. Rg/k, 
k.,/k, k,/k are each of order one, where k is the largest among 9, k,, and ,. 

The amplitudes of the fluctuations have been characterized by the 
non-dimensional parameter «. ‘The length scale of the fluctuations will 
be characterized by another non-dimensional parameter e, defined as the 
ratio of the intrinsic scale of the medium v,/a) to the length scale of the 
disturbances k, that is, 

€ = vgk/ay. (6.4) 
Clearly « can be regarded as the Knudsen number of the disturbances, 
when it is small the inertial effects will dominate over heat conduction and 
viscous effects. ‘The magnitude of « can be estimated for atmospheric 
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conditions for which v, = 0-15cm*sec"! and a, = 30000cmsec"!. The 
intrinsic length (which is simply proportional to the mean free path of the 
gas) 1S v/ag=0-5x10-°cm. Even in the extreme case when the 
characteristic wavelength of the disturbance is as small as 1mm, 
e = 3-14x 10. 

The typical bilateral interaction term is of the order of x? and we shall 
consistently neglect terms of order « and higher. If, however, « < 1 
as is usually the case and especially if « < « we shall find bilateral interaction 
terms that are of the order «e. (‘These usually represent higher derivatives 
associated with viscous or heat conduction terms.) We propose to neglect 
these terms as they will be smaller than the third-order interaction terms. 
Furthermore, if « < « we may even neglect the viscous and heat conduction 
terms in the first-order problem and obtain the simplest set of equations 








oO Qo) 
a = Q, (6.5 a) 
a 
3S 
he 0, (6.5 b) 
A2 PO 2 
<a 8 V2PO = 0. (6.5) 


The superscript (0) emphasizes the difference between these equations 
and (3.5 a), (3.6b) and (3.7b). 

This set of equations was presented by Kovasznay (1953) and represents 
a consistent ‘zero order approximation’ for weak fluctuation fields when the 
time interval of observation is short. ‘The space domain of validity may 
or may not be also small depending on whether or not « increases with 
increasing space domain. (In homogeneous turbulence « may approach 
a limit for increasing domain.) 

Equation (6.5 a) represents a ‘frozen pattern’ of vorticity (or solenoidal 
velocity) and when it is applied to homogeneous turbulence in a wind tunnel 
it predicts an unchanged flow pattern carried by the mean flow velocity 
as suggested by ‘Taylor’s hypothesis. Experimental verification of the 
frozen pattern concept by space-time correlation has been carried out in 
low speed flow for grid turbulence (Favre et al. 1954) and also the outer 
portion of the turbulent boundary layer (Favre et al. 1957). 

Equation (6.5 b) indicates a ‘frozen pattern’ behaviour for temperature 
spots and was similarly used for evaluating hot-wire measurements in 
low-speed turbulent temperature fields. Equation (6.5c) represents the 
well-known equations of sound propagation of the pressure field. 

Equations (6.5) represent the most modest approximation for predicting 
the temporal behaviour of the three modes. ‘Their value in interpreting 
hot-wire measurements in a high-speed flow is equal to that of Taylor’s 
hypothesis as used for interpreting hot-wire measurements in low-speed 
turbulent flows. 

We shall consider the amplitudes of the three modes to have orders of 
magnitude when the velocity fluctuations in the sound mode are comparable 


2K2 
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to the velocity fluctuations in the vorticity mode and when the temperature 
fluctuations in the sound mode are comparable to the temperature 
fluctuations in the entropy mode, Le. |v}| = |vg)!, |T)”| = TY”). 

Space and time derivatives can be estimated from (6.1), (6.2) and (6.3) 
as follows. By choosing the reference space scale as k"! and the reference 
time scale as (a)k)~! all the non-dimensional space derivatives will have the 
same order of magnitude as the quantity which is differentiated (except for 
factors ko/k, k,,/k, k,/k that were assumed earlier to be of order unity). 

The orders of magnitude of the non-dimensional time derivatives differ 
depending upon the ‘hyperbolic’ (sound) or ‘parabolic’ (vorticity and 
entropy) character of the mode. For the sound mode the non-dimensional 
time derivative is of the same order as the space derivative (since the speed 
of sound was used as the reference velocity for forming non-dimensional 
quantities). For the vorticity and entropy modes, however, the non- 
dimensional time derivatives have order of magnitude (k9/k)e or (k,/k 
times the quantities that differentiated. Indeed, pac ie of on 
with respect to ¢ produces a factor vy k% and non-dimensionalization of t by 
(a)k)* results in the factor (v)k&)/(a)k) = (Rg/k)*e. Physically, this is 
clear for the rate of decay of the parabolic modes is an extremely slow process 
when their length scale is very large compared to the mean free path and 
when time is measured in units of the period of sound oscillation. 

The above estimate for the order of magnitudes* of the various quantities 
associated with the sound mode can be summarized as follows: 




















eo TA k ) 
i) ) CY: ) | 
po, ,—* = O(a), ss = Of ae }, 
a’ Ty k 
aP® gpa E avi) 1 dv) 1 aT 1 dT) _ k,, 5 (6.6) 
Om.” Gt "ty 0%: "Ay C2 ° Ty oe ’ T. @ ep 
4 O(1 9 O(1) 2 
asp ASP _ ofl. | 
ox; ’ ot “ae J 


According to our basic assumption in this section, the intensity of the 
vorticity mode must be such that v'g’/a) = O(v),?/a)) = O(a). The order of 
the vorticity properly non-dimensionalized is then Q®/a)ko = O(«). 

The estimates concerning the vorticity mode can also be summarized 
simply as follows: 


¢ Ji 1) er 7 v4 1 ) } 
aie, Of[ — ox |, us O(x 5 
ayko ( k x) Ay (x) 
1 dQ® k2, 1 ov) ko - 
= NO, —- ~~ = O| ~« ), > (6.7) 
ak Ox, : Ay OX; ia 
1 0 QY ke 1 ov 
en Ol — ae }, — — = Of ue }. 
ak ct 3 a ot k2 j 


The terms k,/k, etc., are retained to facilitate arguments in the sub-cases 
where one of the wave-numbers k,, kg and k, is small. 
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Finally, in accordance with the basic assumption, the temperature fluctua- 
tion associated with the entropy mode must be T}""/T, = O(T;”/T,) = O(). 
Now from (2.18) and (3.7a), TS?/T) = SY. It follows that S® = O(«), and 


from (6.3 b), we find _(1) h 
a o( jae). 
Ay k 


Consequently, the velocity associated with the entropy mode is usually a 
very small quantity. ‘This is also clear from (3.7) and (3.7d). Evidently 
this velocity is induced as a result of the diffusion of entropy spots, which 
is a rather slow process when time is measured in units of the period of 
sound oscillation of a comparable spatial scale. 

The estimates for the order of magnitude of the various quantities 
associated with the entropy mode are now summarized* : 











T?) uD k, ) 
S™, => — O a), —— o( = xe), 
é T\ ( ) ay k | 
aS? 1 38T® ' 1 dv k2 ; 
Ox; > Ty “OX; — o(: xe), ay “Ox; _ oF <x), (6.8) 














II 


as) 1 aT? fk 1 av) 3 | 
ol le ee ae 
at, @ k? a, ot k3 J 


By substituting the mode contributions in equations (2.17) the significant 
bilateral source terms can be readily obtained. Since only bilateral 
interactions of the first-order solutions are computed here we shall drop 
the superscript (1) in all subsequent equations. ‘The bilateral source 
terms are 





Mog = 9, (6.9 a) 
a 
yy = —V.(Pp¥,)+ —5- = P3 + Olen), (6.9b) 
m,. = O(ex’), (6.9 c¢) 
MQpy leds V.(P5 ¥o) + O(ex”), (6.9 d) 
oP 
m,, = V.(S,v,)+S, “ar + O(ex”), (6.9 e) 
Meg = V.(S, Vo) + O(ex?) ; (6.9 f) 
foo = —(Vg-V)VoQ, (6.10.a) 
ov, p , 
fp = —P, 7 - }V(v>) + O(ex?), (6.10 b) 
f.. = O(ea?), (6.10) 
fig = —[V(v,-Vo)—V, x Lg] + Ofex?), (6.10 d) 
... : ye ost t O(ex?), (6.10 e) 
fig O(ex?). (6.10 f) 


* Here we assume that the quantity «’ T, is of the same order of magnitude 
as tty. ‘This is indeed the case for air if Sutherland’s formula for viscosity is used, 
(0; /0) To 1-5— Ty/(T,+C) where p, (dp/dT)p,_ 7, and C 140° K. 
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Finally, (5.5) can be simplified to 


Qoo ais i Qo, = O(ex?), (6.1 l a) 
O,, = —(v,.V)S, + O(ca2), (6.11) 
Oxo = —(Vq.V)S, + O(ex?). (6.11 ¢) 


The various types of interactions can now be described in greater detail. 


Vorticity—worticity interactions 

According to (6.9 a), (6.10a) and (6.11 a) the bilateral vorticity—vorticity 
interaction produces mainly a force-like effect. In addition, there is a slight 
heat-like ettect (of order « smaller). ‘The apparent body force is —(v9.V)vo 


or, in tensor notation, 


: 0 0 
[ fool; “25 Oe, Voge dx, [Za:Vaj]- 


Since this force field is in general neither solenoidal nor irrotational, it 
generates both the vorticity mode and sound mode. ‘The production of the 
vorticity fluctuation is of fundamental importance in the mechanics of 
turbulent motion. ‘The production of sound by vorticity—vorticity inter- 
action has been analysed in a general fashion by Lighthill (1952). Our result 
with regard to sound generation could have been anticipated from Lighthill’s 
theory. However, speaking in terms of body force instead of stresses, this 
analysis brings out more clearly the fact that the solenoidal component of 
the force field plays no part in sound generation but is of basic importance 
in the energy transfer between different size eddies in the turbulent 


fluctuation field. 


Sound-sound interaction 
According to (6.9b), (6.10 b) and (6.11 a), bilateral interactions of sound 
modes produce primarily the mass-like and force-like effects. In addition 
there is a slight heat-like effect indicating that only weak entropy fluctuations 
are generated by sound-sound interactions. A typical example of 
sound-sound interaction is the ‘steepening’ of the finite amplitude 
compressive waves and it is well known that the deviation from the 
isentropic change through even a moderately weak shock wave is proportional 
only to the third power of the shock strength. ‘This particular conclusion 
also suggests that the sound-sound interaction may lead to only a slight 
production of the vorticity fluctuations. From (6.10b), we have 
f,, = — PS! — [Vor + O(c). 
Furthermore, the velocity induced by pressure waves propagating in a 
viscous compressible medium satisfies the equation 
= a2V P+ ,%V2v,,, 


so that f,,, can be rewritten as 


f., = —4V[az P? +07] + O(e2?). 
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Consequently, f,,, is principally an irrotational vector field, so that 
sound—sound interaction does produce only a slight vorticity fluctuation. 
There is, therefore, little possibility of altering the intensity of turbulent 
fluctuation in high speed flow through the random interaction of the 
pressure waves and, possibly, weak shock waves. We thus see that the 
main effect of sound—sound interaction is the production of the sound mode 
itself. 

The net strength of the sound source, according to (3.3b), (6.9b), 


(6.10 b) and (6.11), is 


a vs 00, P ‘ 
(5 _ $0}, —V.f, + og V. [v,(V. v,) +(v, -V)v, +45 V(P5)] + 
| dead ‘ 
Kk y) ap |e + O(ex?), 


where use has been made of (3.6d) and (3.6e) (with m= 0). In tensor 
notation, the last result can be written simply as 


72 
0 











0? re = i oe 1 2 29 OP, 2 2 2 
anda | eit + a a aP; 5, | t(y— | (SE) — aP,) | + O(ex?). 


Assuming that the flow fluctuations at infinity are zero, the distribution of 
source strength, when integrated over the entire flow field, in general yields 
a resultant strength different from zero indicating that there is a net generation 
besides self-scattering. (A discussion of self-scattering was given by 
Westervelt (1957).) 


Entropy-entropy interaction 

Our analysis clearly shows that the entropy—entropy interaction, though 
a theoretical possibility, is practically insignificant. From (6.9 c), (6.10 c) and 
(6.11 a), it is seen that it produces only a slight mass-like and a slight heat-like 
effect (both of order of ex”), and an even smaller force-like effect (of order 


€7a:”), 


Vorticity-sound interaction 

According to (6.9 d), (6.10 d) and (6.11 a), the vorticity-sound interaction 
produces a mass-like and force-like effect. It has only a slight heat-like 
effect and hence produces only very slight entropy fluctuations (i.e. of 
order of €x”). 

Since the force field in general is not an irrotational field, there will be 
vorticity production as a result of interaction. ‘The source strength in the 
production of vorticity fluctuations is V x (v,, x Q 9). This term includes 
the effect of both the stretching and expansion of vortex tubes and the 
convection of vorticity under the action of sound fluctuation. In the case 
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when a single shock wave is interacting with the turbulent flow, the 
theoretical calculation of Ribner (1954), which was explored experimentally 
by Kovasznay (1955), indicates that the effect is not very large, at least in 
that particular problem. 

The vorticity-sound interaction also produces.the sound mode. The 
resultant strength of the sound source as given by (3.7) can be simplified to 


d 2 Qo, ; 
(5; — 3% ‘mop ~WV fo, + > 2 = 2V. [(Vg.V)v,] + O(ex?). 





In the last step of the derivation, use has been made of the relation 


oP 
at = —V.v, + Olea). 





In tensor notation, the resultant source strength of the sound mode is simply 
2 0*(Vo;V,,;)/0x;0x;, a result which could have been anticipated from 
Lighthill’s quadrupole radiation theory. Assuming that the flow fluctuations 
die down fast enough at infinity, integration of the source strength over the 
entire flow field yields no net contributions. Consequently, the sound mode 
produced by vorticity-sound interaction can be interpreted as the scattering of 
the incident sound wave by the vorticity fluctuations. Scattering of a shock 
wave by turbulent fluctuations has also been calculated by Ribner (1954). 


Sound-entropy interaction 

By (6.9 e), (6.10 e) and (6.11 b), the interaction of the sound and entropy 
modes produce all three basic modes of fluctuations. 

Let us first examine the vorticity mode produced by the interaction. 
According to (3.3a) and (6.10e), the source term which produces vorticity 
fluctuation is 


V x | s. | + O(ea”) = VS, x = + O(ex?) = —a, VS, x VP,. 

‘This is a familiar term in the theory of rotational flow. Physically, the 
vorticity produced by it can be interpreted as the consequence of the 
fact that the resultant force acting on any lump of fluid fails to pass 
through the centre of mass, and thus produces a torque and an angular 
acceleration. ‘This result suggests that the action of pressure field on a 
stream carrying randomly distributed temperature spots in a turbulent 
stream may generate turbulent fluctuations and can modify the original 
turbulent fluctuations in the stream. In the only known experimental 
realization of this situation (spottily heated air accelerating through the 
nozzle of a supersonic wind-tunnel) the effect appeared unimportant (see 
Morkovin 1956, p. 8). 

Next consider the sound mode generated by the interaction. ‘The net 
sound source can be calculated from (3.6 b) and can readily be simplified to 


é oO 7 re 
= 3 v2), Verbs t — = v.56. v) | + O(ex). 
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If the disturbance dies down at infinity fast enough, the resultant source 
strength integrated over the entire flow field is zero, so that the phenomenon 
can again be interpreted as the scattering of the incident sound wave by the 
entropy spots. 

Finally, consider the generation of the entropy spottiness by the inter- 
action. It follows from (3.7b) and (6.11b) that the source of entropy 
fluctuation is given by the term —(v,,.V)S,. It can be regarded as a 
convection of heat by the motion of the sound field. 

Special cases of generation of all three fluctuation modes in the interaction 
between an entropy disturbance and an infinite or wedge-attached shock 
wave were discussed theoretically by Chang (1955). 


Entropy-—worticity interaction 

According to (6.9 f), (6.10) and (6.11 c), the entropy—vorticity interaction 
produces a mass-like and a heat-like effect as well as a slight force-like effect. 
Consequently, this type of interaction can only produce a very weak (of 
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Second-order non-linear interaction between modes of order «a. 


Table 1. 
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order ex”) vorticity fluctuation. ‘lhe sound field produced turns out to be 

weak also because the expansion generated by the apparent mass addition 

is exactly balanced by the contraction following heat subtraction. ‘his 

can best be seen if we calculate the net source strength for the sound mode 

produced. From (6.9f) and (6.11 ¢), one concludes that 

(; -- Ne hows Gro = - V. [vo S,] 
) 


ot 


| 


(vo .V)S, + O(ex2) 





) 


ct 

= O(ex?), 
since V.vo =0. Consequently, the main effects of entropy—vorticity 
interaction are in the production of the entropy mode. ‘This mode is 
generated as a result of the convection of the hot spots by velocity fluctuation 
associated with the vorticity mode. ‘This type of interaction has a leading 
role in turbulent heat transfer phenomena, see Corrsin (1952). 

In order to summarize all second-order effects table 1 was prepared. 
It gives the source strength for all three modes up to order «. ‘Terms 
of order x and #?e? are only indicated. ‘Table 1 includes also a more or less 
generation’ or ‘ scattering’. 


‘ 


intuitive label for each bilateral interaction, like 
‘These are not intended as accurate descriptions but are merely to suggest 


the corresponding physical processes. 
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ty 
IS SUMMARY 
- A simple approximate model is set forth for the flow field 
1g 


5 between the nose of a blunt body of revolution and its detached 
shock wave. ‘The model tends to explain the poor convergence 





a. of Chester’s solution, which is based on an improvement of the 
- Newtonian approximation. It suggests a modification of his series 
* for the body shape which appears to improve its convergence 
considerably. 
st ‘ 
1. INTRODUCTION 
S Chester (1956) and Freeman (1956) have recently advanced in this 
x journal an ingenious attack on the problem of supersonic flow past a blunt 
body. (See also Lighthill 1957, $3.5 et seg.) ‘Their point of departure 
t is the ‘ Newtonian-plus-centrifugal’ solution, which becomes exact as the 
free-stream Mach number M approaches infinity and the adiabatic index y 
approaches unity. After modifying this basic solution to be valid near the 
body, Chester and Freeman improve upon it by successive approximations, 
the key to their success being the adoption of von Mises’ transformation, 


in which the stream function is taken as one of the independent variables. 
The result of Chester (who has carried the approximation several steps 
further than Freeman for a perfect gas with constant specific heats) is a 
double series expansion in 6 = (y—1)/(y+1) and M™. 

Unfortunately, the series for the flow variables of interest appear to 
converge poorly and, surprisingly, less well for axisymmetric than plane 
flow. For example, for axisymmetric flow of a perfect gas with y = 7/5 
at M = o past whatever body it is that supports a parabolic detached 
shock wave, Chester’s series for the stand-off distance A of the shock wave, 











in terms of its nose radius R., is 
AJR, = 1(1 —0-6667 + 0-4333 — 0-3062 +...). (1) 


Although the expansion parameter 6 is here only 1/6, this series converges 
so slowly as to be of little practical value. (Chester has used his series only 
out to 6 = 1/10, and then attempted to extrapolate to such ‘large’ values 
as 1/6.) 

The question naturally arises whether this discouraging behaviour 1s 
symptomatic of some basic defect in the approximation, or simply indicates 
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that the series has an unduly small radius of convergence (if it in fact 
converges). ‘Che present note aims to suggest that the latter is the case, 
and that by appeal to a simple approximate model of the flow field the series 
for the shape of the body can be recast in a form that appears to converge 
much more rapidly. 


2. A SIMPLE MODEL OF THE FLOW FIELD 

Consider axisymmetric flow of a perfect gas with constant y behind 
a prescribed shock wave. For simplicity we take the shock wave to be a 
parabola (which is the only case considered in detail by Chester), but other 
shapes could be treated similarly. According to the Cauchy—Kowalewski 
theorem, the flow field is analytic in some region downstream of the shock 
wave, and there is no reason to doubt that this region extends to the surface 
of the body. One can therefore attempt to expand the flow field in ‘Taylor 
series starting from the shock wave. It is convenient to work with Stokes’s 
stream function % and parabolic coordinates (€,7) which are normalized 
such that the shock wave is described by y = 1. 

On the downstream face of the shock wave, the stream function has the 
free-stream value, its first derivatives are given by the Rankine—-Hugoniot 
shock relations, and its second derivatives can be found by substituting 
into the equations of motion. The present model consists of the resulting 
first three terms of the ‘Taylor expansion starting from the shock wave: 


l y+1)M? M?—-1\2 —1)M?-—2(2+ M? 
Bt A yg a a 














a 2 A = sf 2A 
M?+1.. A(y + 1) M4 

ra 24 : Wa, —1\)4e2 

2 + OT AE + 2(y — 1)Mt& x 


2yM? —(y—1)(1 +) . 
(i+2)48 lc —7), (2a) 


where A = 214+ &)+(y—1)M?. (2b) 





to 


| 
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y+ et Ae. , 
a a oe —7) | =p +31 hop ra | (3) 


In the Newtonian limit, y > 1 (and » —- 1), this reduces further to 


2 oe. 

~3(1 21) the (7) (4) 

¢ = Y Gey | 
On the axis (€ = 0) the point at which this last expression vanishes (corre- 
sponding to the nose of the body) lies at (1—7) = 3(y—1), which is just the 
result of the first approximation of Chester and Freeman. It is mainly 
this attribute that qualifies the model as a useful one. Off the axis (4) 
cannot vanish, which suggests that the model will prove less faithful away 











from the nose of the body. 
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It should be remarked that the idea of expanding the flow field in ‘Taylor 
series from an assumed shock wave is by no means new, but has been followed 
by a number of workers with uniformly poor results. ‘The innovations in 
the present model are, first, that the natural coordinates associated with 
the shock wave are employed, and, second, that the ‘Taylor series is terminated 
with the parabolic approximation to the stream function; and both of these 
are essential. Previous investigators, of whom the most assiduous has been 
Cabannes (1951, 1956), all carried out double ‘Taylor series expansions in 
Cartesian coordinates. Aside from the inferiority of a double series to the 
present single expansion, this yields series whose first few terms approach 
the true result much more slowly than the present series in parabolic 
coordinates. ‘This is illustrated in figure 1 for a parabolic shock wave at 
M = 2(with y = 7/5), for which Cabannes (1956) has computed seven terms 
of the double ‘Taylor series. 
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Figure 1. Taylor e¢ries expansions for stream function on axis of symmetry behind 
parabolic shock wave at M = 2,y =7/5. (a) Cartesian coordinates; (5) Parabolic 
coordinates. 


Figure 1(4) also illustrates the second point that, even in the natural 
coordinates, adding further terms to the Taylor series worsens rather than 
improves the approximation. The reason for this is that analytical continua- 
tion of the flow field upstream through the detached shock wave leads to 
a limiting line, or envelope of characteristics, as sketched in figure 2. The 
continued flow field has a square-root behaviour there. Since the ‘Taylor 
series represents this analytical continuation, the limiting line determines 
its radius of convergence. On the axis, the limiting line is about as far 
ahead of the shock wave as the body is behind it, so that the convergence is 
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marginal; and in the example of figure 1 it is actually closer, so that the 
series diverges at the body. (‘The divergence can be ameliorated by applying 
the transformation of Shanks (1955).) Off the axis the limiting line is 
even closer, relative to the distance of the body, so that the convergence 
deteriorates further. 


Limiting 


7 line 
Outgoing a 
P 


characteristic, 


Body 





7 
AX | 


Figure 2. Schematic diagram showing analytical continuation of flow upstream 
through shock wave. 


Poor convergence, or even divergence, is not a fatal flaw, however; 
though it is for this reason that the truncated series is regarded as a model 
rather than a systematic approximation. ‘The variation of the stream function 
is actually very nearly parabolic, so that (somewhat as with asymptotic series) 
the first three terms of the Taylor series provide a close approximation 


near the axis. 


3. APPLICATION TO CHESTER’S SOLUTION 
3.1. Body shape at infinite Mach number 
Consider first the case of infinite Mach number. Requiring the model 
stream function (3) to vanish determines the body according to 
| fo) 
( 7 )model ! + (38(1 —8) — 2(1 — 8)€?/(1 +4 &)]} 2? 





— 
on 
— 


where 6 = (y—1)/(y+1). For comparison with Chester’s solution, we 
transform to cylindrical polar coordinates in which the shock wave is 
described by x = 377, which means that x = $(1+€?—7?) and r= én, 
and then expand for small 6. ‘Then (with fractions unreduced to facilitate 


comparison) the body is described by 


2 yak Fd 
Ymodel of 1 V/ (98/3) 4 To $- 168 °% (Jo 3) + | + 
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whereas Chester’s result, of which the leading term was first given by 
Hayes (1955), is 


a 26. 463 gi 
x=0 1 \ (d0 Ss) t 10° 168 oO* (dod 3) are =) 








l fy 19. 47. an 11173 

— eel 4 — ( an = < ji j = \~ 
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- [36 874/(85 3). | (7) 
The agreement is remarkable for the terms independent of r, though (as 
anticipated) less satisfactory for those in r?. ‘The model series converges 
for 56 smaller than 3/9, which suggests that Chester’s series will converge 
for 5 less than 3/8, corresponding to y less than 11/5, a value not attained 
in reality. However, the convergence of both model and prototype is 
unsatisfactory at 6 = 1/6, corresponding to y= 7/5. In the model the 
poor convergence results from expansion of [1 +(36)!?]-' by the binomial 
theorem, which suggests that the convergence of Chester’s series might be 


] 


considerably accelerated by recasting it as a series for (1 —7) Doing so 





yields on the axis (€ = 0) 


o>) 


— eee l. 187 ~— ; 
(1-7) 0 1+ 4/ (80/3) 30° t 340 °% (00/5)+... |. (Sa) 


For 6 = 1/6 this gives 


(1—n)7= 6(1+0-6667 — 0-0722 + 0:0247 +...), (Sb) 





which appears to be a great improvement over the result of the original 
series for the related quantity A/R, given in equation (1). 


3.2. Body shape at high Mach number 

Consider now Chester’s general case of high Mach number. For 
simplicity, we restrict the model to the axis, where it is most faithful. “Then 
from (2a) it locates the nose of the body at 


\ 


(1—n)moaer = [2+ (y—1)M2][(y + 1)M2 4 {6(y — 1)M44 





In terms of Chester’s parameters 


6 = (y—1)/(y +1), d 6+ M-, 





this is 
(1—7)mode = (8+ M-2— M-%5){1 + [(1 —8)(35 + 4M-4 — 3M-48)}#23-1 (9b) 
or 


(1—17)moaer = d(1 — M2 + M-4/d){1 + [(1-d+ M~)(3d—30 
-4M-+4—3M-d+3) 


Now the last form can be expanded in a series containing only positive 


7) A (9c) 


powers of d'* and W-°, as in Chester’s solution, only under the assumption 
that M-* = O(d?), though Chester assumes merely that M-? < d. This 
minor discrepancy probably indicates a certain lack of fidelity in the present 




















520 M. D. Van Dyke 


model. ‘To play safe, however, we shall treat Chester’s solution only under 
the additional assumption that W~? = O(d?). For y = 7/5, this restricts M 
to values greater than perhaps 4 or 5. 

As before, the model suggests recasting Chester’s series as an approxima- 
tion for (1—7)'. The result, including now terms in &?, is 


i ee : | 
_— dj 1 IF) uae s0¢ | (3) | ; 
(7 1) Uf8d\ 299. M2 
: E a" al (5) i008 7+ 3a 
27269 8d 
- somo 4s/ (3) oe (10) 


(Note that on the axis the solution to this degree of approximation does not 
depend upon ™ explicitly, but only in so far as it appearsind.) ‘The apparent 
convergence of the terms in é*, though again not as rapid as that of the terms 
independent of €, is a great improvement over that of the original series, 
and probably adequate for practical purposes. 

For y = 7/5, this modified series (10) leads to ratios of stand-off distance A, 
initial shock-wave radius R,, and body nose radius R,, which are compared 
in table 1 with the results of the model. 
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Table 1. 


For M © an accurate numerical solution of the full equations has 
been found as the first part of a current programme of calculating the flow 
field behinda family ot detached shock waves. The results are A/R, = 0-0988 
R,/R, = 0-726, and A/R, = 0-1360. Hence the modified series predicts 
the ratio of stand-off distance to body radius to within a fraction of 1%, 
and even the model is only 2°, off at / = co. ‘The numerical solution also 
indicates that the axisymmetric body that supports a parabolic shock wave 
at J = oo is itself almost exactly parabolic (whereas in plane flow it is 
exceedingly close to circular). 

‘he tabulated values of stand-off distance are compared in figure 3 
with a number of measurements on spheres in air under conditions such that 
the adiabatic exponent is close to 7/5. Schlieren photographs show that 
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for spheres flying at Mach numbers between 2 and 7 the detached shock 
wave is nearly parabolic over the region that determines the stand-off 
distance, so that the comparison is justified. At lower Mach numbers the 
shock wave on a sphere departs significantly from a parabola, and the 
assumptions of Chester’s theory have also been greatly exceeded, so that the 
agreement shown in figure 3 is undoubtedly coincidental. (No values are 
shown for Chester’s original series, since it converges so poorly that he did 


not himself suggest using it at y = 7/5.) 
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Figure 3. Stand-off distance of shock wave from sphere. 
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3.3. Pressure distribution 
Chester’s series for the pressure at the surface of the body also sutters 


I 
from poor convergence (except at the stagnation point). ‘Thus at M L 
it gives, in units of p, V* 
# | l Fe a9 32 Pee 
p=({!1 zo + 5] ry? Z3 0 g OV (80 5) + 3360 5 | (11) 


The corresponding result for the model stream function is found, using the 


Bernoulli equation and the condition of constant entropy along streamlines, 


to be 
| co. a 747 

4 | ) A2 2} ( GY, (98/3 0* ) 

7 i ? Q ] } | D 4 \ IO}. ) 16 ) ‘ (12) 

The agreement is fair. In particular, the model confirms that the last 


coeficient is enormous, though Chester’s is still greater by a factor of four. 
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Unfortunately, the model does not in this case lead to a useful modification 
of Chester’s solution. It actually suggests that his series should be used 
to calculate the quantity »%p?(u? + v7), and that the pressure then be found 
from the Bernoulli and entropy relations; but working back from Chester’s 
solution for p gives a series for that quantity whose convergence is not 
appreciably improved. 

In the case of plane flow, the convergence of Chester’s series is nearly 
adequate. It would nevertheless be of interest to construct a model, if only 
to explain the source of the logarithmic terms. (Successively higher powers 
of logarithms in perturbation series may be the asymptotic representation 
of Bessel functions, inverse hyperbolic functions, small fractional powers, 
etc.). However, the author has been unable to devise a useful model, or 
otherwise to discover an appropriate transformation of the series. 
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The effect of a vertical lapse rate of temperature on the 
spiral flow of a fluid in a heated rotating cylinder 


By G. N. LANCE 


Computation Laboratory, University of Southampton 
(Received 8 October 1957) 


SUMMARY 

Experiments have been performed with a viscous fluid contained 
in a rotating cylinder heated from below. Previous theoretical 
explanations of the flow patterns obtained in these experiments 
assumed that the convective terms in the heat transfer equation 
were negligible. ‘The present paper gives a treatment which 
includes one of these convective terms. ‘lhe results confirm the 
physical reasoning that, when the lapse rate is positive, the stability 
is increased and that the motion is therefore decreased. 


1. INTRODUCTION 

Many experiments have been performed with a fluid in a rotating 
dish-pan which is heated from below. ‘lhe case of a shallow fluid has been 
covered very fully by Fultz (1951, 1956). Several authors have put forward 
theories which attempt to explain the results. Davies (1953) used a theory 
which neglected convective terms in the heat transfer equation, whereas 
Kuo (1954, 1955, 1956) included both horizontal and vertical temperature 
gradients, but both these theories needed, for their success, the assumption 
of a shallow fluid. 

Experiments with deep fluids have been described by Fultz (1956). 
Skeib (1953) performed experiments with a value of a, which is effectively 
a depth parameter (see $3), very nearly equal to unity. Unity is the value 
of a which is used, in this paper, to obtain the curves shown in figures 1-3. 
Davies’s theory was extended to cover larger depths by Lance & Deland 
(1955) but they neglected the convective terms in the heat transfer equation 
(in the present paper this reference will be referred to as 1). 

in the dish-pan experiments the fluid was heated near the circumference 
of the base of the dish-pan and cooled near the axis of rotation. All the 
theories so far put forward have used the method of small perturbations, 
the perturbations being superimposed on a steady state of rotation without 
heating and at constant basic pressure, density and temperature. Previous 
theories, with the exception of Kuo’s (1954), nave assumed that the vertical 
temperature lapse rate was zero but the experimental results have shown 
that in fact the temperature varies in a manner which is very nearly a linear 
function of height. ‘The purpose of the present paper is to determine the 
effect of this vertical lapse rate and to see whether the flow patterns previously 
given in I are greatly altered by the lapse rate. 

2L2 








524 G. N. Lance 


2. "THE EQUATIONS OF MOTION 
‘The notation used here is the same as that employed in I, where the 
equations of motion were presented in full. If it is assumed that the motion 
is axially symmetric and that a steady state has been reached, then the 
equations can be written, using cylindrical polar coordinates (r, 0,3), in 
the following form: 





2 OP 1a : 
a(pe-S)=- f+ V4 — +e, mn 
Dv, - U,V; _ ‘ v . 

alae + 37) ~ o(t- 3) (2) 
= = 2. 0 >] 
De =~ ag BOHR a (3) 
Dp 

Dr ee by ’ (4) 
Pr = po H(T, — To), (5) 

DT, Ory. : 
Pid eye =JSkV?7 1 +, + Ph. a 


In the above six equations, u,, v,, w, are the components of fluid velocity 
in the directions of 7, 6 and = increasing, respectively, relative to an axis 
system fixed in space with the origin at the centre of the base of the dish-pan. 
The pressure is denoted by p, and the density by p,; yu is the coefficient of 
viscosity, « is the reciprocal of the coefficient of cubical expansion, J is the 
mechanical equivalent of heat, c, is the specific heat at constant volume, 
kis the thermometric conductivity and Q, is the viscous dissipation function. 





Also, D r) 0 . 02 1. (oe 
— Uy ON Sag V2 == +- = mate 
Dt oz or? ror. og? 
and ou, u, CW, 
by = [Goce | ms 
OF r OR 


When the liquid is not heated, rotation of the cylinder causes solid 
rotation of the fluid and, if heat is supplied, the departures from the state 
of solid rotation will be small. ‘Thus, if Q is the angular velocity of the 
cylinder, we write 

uy, =U, v = OQr+y, em, Pi = Pot Ps 7) 

/ 

Pi=Pot?P, T, = T+ T. " 

The symbols without suffixes are the six perturbation quantities which 

determine the departure from the state of solid rotation. It was shown in I 

that, if attention is confined to slow rates of rotation (as in the experiments), 

the upper surface of the liquid is approximately the horizontal plane z = h, 
where h is the depth of the fluid. 

When the expressions (7) are substituted into (1) to (6), the following 
equations are obtained for the perturbation variables: 

—20p uv = —dp/or+ p(V2u — u/r?), (8) 
+ 20p)u = p(V20—v/r?), (9) 
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op/dzx— gp + pV2w = 0, (10) 
rps, = A(ru)/or + a(rw)/dz = 0, (11) 
p= —al, (12) 

V?T = (pc,/k)DT/Dt, (13) 


where all non-linear terms have been neglected. Now, for the purposes 
of the present paper, it is assumed that 07,/dz = © = const., so that (13) 
reduces to 


Ve (pc, k)Ow, (14) 


the work of I simply took © = 0. 

The term containing 07,/cr has been omitted and there is, in fact, no 
real justification for such an omission; however, it does enable the 
variables r and z to be separated and in this way the equations can be reduced 
to ordinary differential equations (see (17) to (20) below). Moreover, 
such a procedure has the advantage that the effect of only one advection 
term may be obtained and their relative importance may be determined. 
Finally, if the density is eliminated from (10) and (12) it is found that 


— dop/oz+gaT +pV?w = 0. (15) 
3. REDUCTION OF THE EQUATIONS TO A SYSTEM OF ORDINARY DIFFERENTIAL 


EQUATIONS 


It was shown in I that the variables y and z can be separated by the 


substitutions 
I(r, 2) = T,(2)J,(Br), 
u(r, 2) = U,(z)/,(Br), 
u(r, z) = Vi(s)4,(Pr), (16) 
w(r, z) = W,(z)J,(Br), 
p(r, 2) = P,(2)J9(Br), 
(r,s) = I(2)J,(Br), 


where fr, is a zero of J ,(Bry) = 0 and ry is the radius of the dish-pan. 

It is only possible to retain the vertical lapse rate in (14) because the 
r variation of both T and w is of the form./,(8r). If the form were not the 
same for both these variables the r coordinate would not cancel from this 
equation and the system of partial differential equations could not be 
reduced to a system of ordinary differential equations. As in I, all the 
quantities may be written in non-dimensional form and the pressure can 
be eliminated from (8) and (10). The final system of equations is 


U"~@U'+2RV' +aW" —a@W+aT = 0, (17) 
V" = 2RU+@V, (18) 
W' = —aU, (19) 


T" —a@T = (h?00/v)W = AW, (20) 
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where o is the Prandtl number, v is the kinematic viscosity, a = Bh and 
R is the rotational Reynolds number Qh?/v; the dashes denote differentiation 
with respect to € (= 2/h). 

The boundary conditions are the same as those used in I, since they 
do not depend on ©, namely 


U(0) = V(0)=WO0)=0, TO) =1, 


U"(1) = V'(1) = W(1) = T'(1) = 0. (21) 


4. METHOD OF SOLUTION 

Despite the fact that the system of differential equations (17) to (20) 
is a linear one, it is not practical to solve them without recourse to numerical 
methods. A numerical method was also employed in I but, since then, 
a better method for the solution of a system of linear differential equations 
subject to two-point boundary conditions has been described by Goodman 
& Lance (1956). The improved method was used for the present com- 
putations and found to be very satisfactory. ‘The calculations were 
performed on a Ferranti Pegasus digital computer. 

Numerical values are required for a, R and A © (h?o@/v). The value 
of a was taken to be unity. This implies a depth of liquid of about + cm 
in a dish-pan with a radius of 15cm. Since fry = 3:83 and a = Bh, 
h = r,/3-83 and so, under the conditions of the experiment, for which 
r,= 15cm, o =7 for water at a temperature of 15°C, v = 0-01 and 
© = 1, it follows that A = 10737. In the calculations the Reynolds 
number R was taken to be 0 (no rotation), 2, 4, 8, 16, 32 and 64. 


5. RESULTS AND CONCLUSIONS 
Radial velocity component 
Figure 1 shows the curves U(€) x 10 against €, fora = 1 and A = 10737. 
The curves shown are for R = 0, 16, 32 and 64. The results for R = 2, 4 
and 8 were computed, but they are omitted because they lie between those 
of R= 0 and R= 16. This figure is to be compared with figure 2a of I, 
which shows the curves when the lapse rate is zero, i.e. A = 0. 

The following differences may be noted. 

(i) In figure 1, U x 10* is plotted, whereas in the corresponding figure 
of I, Ux 10% is shown. Hence, there is a considerable drop in the 
magnitude of the radial velocity component. 

(ii) In the present case, for each value of R there is a zero of U when 
€ < 0-5; previously, the zero occurred when € > 0-5. 

(iii) For small R the greatest inflow (U < 0) is slightly be/ow the surface 
when A + 0, whereas when A = 0, the maximum inflow is at the 
surface. For large R there is a band, for which 0:25 < € < 0-75, 
of comparatively slow inflow, but when € > 0-75 the inflow increases 


and reaches a maximum at the surface. 
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(iv) In the case of R = 64 the radial motion tends to be concentrated 
in the upper and lower regions of the fluid. Such a behaviour 
was not apparent in the results of I. 


UxIo* 
+6 J 


+5 
+4 
+3 


+2 
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Figure 1. Non-dimensional radial velocity profile plotted against non-dimensional 
height. Cases for which 4 10 737 are shown as solid lines. ‘The broken 
curve is 4 0,R 16 and is the value of U » 103. 


Zonal velocity component 
Figure 2 shows the curves V(€) x 104 against €, fora = 1 and A = 10737. 
The curves shown, for R = 0, 2, 4, 8, 16, 32 and 64, may be compared 

with those of figure 3a in I. ‘The following points are notable. 
(i) As in the case of U(&) there is a reduction in magnitude by a factor 

of 10. 

(ii) The band of easterlies is confined to € 0-2 when A = 10 737. 
(iii) The maximum zonal velocity at the surface continues to increase 
as R: increases, whereas when A = 0 the zonal velocity reached 
a maximum value when R = 16. (This is only true for R 64, 
which was the limit to which the calculations were taken.) 
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Figure 2. Non-dimensional zonal velocity profile plotted against non-dimensional 


height. Cases for which A 10 737 are shown as solid lines. The broken 
curve is 4 0, R 16 and is the value of V’ « 10° 


Vertical velocity component 

Figure 3 shows the curves —W/(€)x 10° against € for a=1 and 
A = 10737. The curves for R = 0, 16, 32 and 64 are shown, the remainder 
are omitted for simplicity. The case A = 0 is figure 4a of I, and the 
following are the differences between the two cases. 

(i) ‘The scale factor is 10° in the present case as opposed to 10° in 
the case of A 0, 

(ii) The curve R= 64 shows a distinct flattening in the range 
0:25 <— € < 0-75, when A = 10737, which implies that w is 
tending to a constant value at intermediate depths. 

(iii) For each value of R the maximum vertical velocity is below € = 0-5 
when A 10737, but when A = 0 the maximum is in each case 


ibove & Ps 


lhe results, when looked at in toto, show that when a positive vertical 
lapse rate is present, the whole velocity field is reduced in magnitude. 
[his result agrees with physical reasoning because, if the upper layers of 
the fluid are the warmest, then greater stability would be expected and 


consequently less motion. 
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Since the results obtained here are an order of magnitude smaller than 
those obtained previously it should be noted that on all the diagrams the 
comparison curves, which are taken from Lance & Deland, are shown 
dotted and are plotted to a different scale. Enough has been said in the 
above remarks, about the results, to make the differences apparent and the 
interested reader is referred to Lance & Deland for complete diagrams of 


the case 4 = 0). 
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Figure 3. Non-dimensional vertical velocity profile plotted against non-dimensiona} 


IZul i i hhe-(€ll?) 
height. Cases for which 4 10 737 are shown as solid lines, The broken 
‘urve is 4 ) R 16 and is the value of WW 103 


A detailed comparison with the experimental results of Fultz and Skeib 
has not been made because, in all the experiments, the value of the Reynolds 
not correspond with any of those used in the present 

calculations. It is obviously very desirable to extend this work to higher 
f R and to other values of the depth parameter a. Such an 
at present in hand but it was considered worth while to 
] 


publish these results as being a valuabl 


number does 


investigation 1s 


extension of the work described 


thor wishes to express his thanks to Mr J. V. Davies for his 
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Properties of the Hugoniot function 


By ROBERT D. COWAN 
Los Alamos Scientific Laboratory of the University of California, Los Alamos, New 
Mexico 


(Received 22 August 1957) 


SUMMARY 


Differentiation of the Hugoniot function 
H(p,v) = E(p, v) — E( Po, %) + 3(P + Pow — %) 
and use of the first and second laws of thermodynamics leads to 
the relation dH = T dS+dA, where dA is the element of area 
in the (p,v) plane swept out (in a counter-clockwise direction) by 
the line segment (pp, v)) > (p, v) as the point (p, v) is moved from 
some point (p,,v,) to a neighbouring point (p, + dp, v, + dv). 
This relation, together with rather general assumptions regarding 
the shape of the isentropic curves dS = 0 for the material behind 
the shock, makes possible the geometrical derivation of a number 
of properties of the function H and of the Hugoniot curves dH = 0. 


1. INTRODUCTION 


‘The laws of conservation of mass, momentum, and energy lead to the 
following relation (first derived by Hugoniot (1889); see also Courant & 
Friedrichs (1948), pp. 116-146 and 204-218) which must be satisfied across 
any shock front: 


E— Ey, = 3(P +Po%o- 2); (1) 


where py, vp, and E, are the pressure, specific volume, and specific internal 
energy of the material ahead of the shock front, and p, v, and EF are the 
values of these quantities behind the shock front. (£ and Ey must, of course, 
be referred to the same energy zero—for example, the energy of the elements 
in the ideal-gas state at zero temperature.) If the equation of state of the 
material behind the shock is known so that FE and p can be found for any 
state, then (1) defines all states which can be reached from (Pp, v%, Ey) by 
means of a single shock. ‘The locus of all such states can be represented 
as a curve in the p-@ plane, known as the Hugoniot curve. 

Even without knowing the equation of state in detail, it is possible to 
derive several important qualitative properties of the Hugoniot curve on 
the basis of certain rather general equation-of-state characteristics. ‘This 
problem has been discussed at length by Bethe (1942) and by Weyl (1949) 
for the case'in which the materials behind and ahead of the shock are 
chemically identical, and by Courant & Friedrichs (1948) also in the case 
in which the shock induces a chemical reaction (as in the detonation of 
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explosives). However, though the assumed equation-of-state characteristics 
are essentially of a geometrical nature (in terms of the p-wv plane), the 
treatments of the above authors are primarily analytical. ‘The result is 
that one does not readily obtain much of a feeling for the geometrical 
relations between the original assumptions and the final deductions. ‘The 
treatment given below is basically geometrical, and the writer has found 
it quite useful in dealing with several problems which he has encountered, 
two of which are mentioned at the end of the paper. Many of the properties 
discussed below are well known, but are included for the sake of 
completeness; the proofs given serve to illustrate the simplicity of the 
geometrical approach. 

‘l'wo of the basic assumptions made by Courant & Friedrichs in treating 
the reaction case are rather obscure; these points are cleared up in the 
discussion following equation (2) and Property 15. 


2. ASSUMPTIONS 

It is assumed that there is no external electromagnetic field present, 
which might affect the problem in the case of shock-induced ionization. 
The state of the material behind the shock (which may be a mixture of 
several compounds) is then described completely by the internal energy 
function E = E(S,v,n;), where S is the specific entropy and n; the number 
of moles of compound 7 per gram of mixture. The difference in internal 
energy between two neighbouring states is 


dE = T dS — pdv+ > pL; dn,, (2) 


where py, is the chemical potential of the ith component. It is assumed 
that any shock-induced reaction takes place rather rapidly, and we shall 
consider only those states which occur to the rear of the reaction zone, 
where the material is essentially in chemical equilibrium. ‘The term 
> pu; dn, in (2) is then zero, the change in internal energy becomes simply 
dE = T dS-—p dv, (3) 
and the internal energy can be considered as a function of S and wv only. 
Assuming further that 
(op/eS), > 0 (4) 
so that no two isentropic curves ever cross, the energy can be considered 
as a single-valued function of pressure and volume 
E = E(p,2), (5) 
which will be supposed to be everywhere continuous with continuous 
derivatives. (Assumption (4) is invalid for water between 0° and 4 C, 
and for phase changes generally.) Finally, it is assumed that all isentropes 


have negative slope and positive curvature: 
(cp/ev)s < 9, (6) 
(0?p/dv?) 0. (7) 
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These two conditions are satisfied rather generally; they have been 
investigated at some length by Bethe (1942) for substances not reacting 
chemically and shown to be valid usually except at some phase trans- 
formationst. However, (7) is occasionally violated in the case of chemically 
reacting substances. 

From (3) and (4) it follows that 


Since (4) and (6) imply 


(aS/év), > 0, (9) 
it follows from (3) that 
(2E|20), > —p. (10) 
In some cases for p > py and v < vy (which will be mentioned specifically 
below), the stronger condition 
(0 , dv), Pass (p + po) (11 a) 
or (from (3)) 
(0S Ov), > (p—Ppo)/T (11b) 
will be assumed to hold, and in other cases, the still stronger condition 
(cE/dv),, > 0 (12a) 
or 
(0S/ov), > p/T (12b) 


will be applied. Even this last inequality (which, with (8), is equivalent 
to (dp/dv), <0) has been shown by Bethe to be usually valid except at 
some phase transformations. 

It is assumed throughout that py) > 0. 


3. DEFINITIONS 
In order to investigate the properties of the Hugoniot curve defined 
by (1), it is convenient to define the so-called Hugoniot function 


H(p,v) E(p, v) — E(Po, %) — 3(P+Po)(Yo— 2); (13) 
where E(p ,@ ) is the internal energy which the material behind the shock 
would have if it were returned (always under conditions of thermodynamic 
and chemical equilibrium) to the pressure and specific volume (fo, Up) of 
the unshocked material. By assumption (5), H isa single-valued, continuous, 
and differentiable function of p and vz. In terms of the Hugoniot function, 


the Hugoniot relation (1) becomes 
H = E,( Po. %) — E(Po, Uo) = constant. (14) 


+ The inequality examined by Bethe is actually 
( 0p/év) gz (ép/ ev) g+(p/T)( ep/ 2S), < 0, 
which is less general than (6) if (4) holds. (However, Bethe’s third assumption 
(@p/aS), >—2T/v is more general than (4).) 
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‘Thus the Hugoniot curve is a curve of constant H, the constant value being 
zero if the shock induces no chemical reaction (or a reaction such as 
dissociation or ionization for which the shocked material can be returned 
reversibly and spontaneously to the original state of the unshocked material), 
and the constant being greater or less than zero if a reaction is induced 
which is exothermic or endothermic, respectively. (Note that if the 
unshocked material is, for example, a complicated organic high explosive 
in a condensed phase, it may be impossible to reduce to py the pressure of 
the reaction products at volume vy even by reducing the temperature to 
absolute zero, so that E(py,v)) is undefined. In this case, the zero- 
temperature isotherm passes above the point (Pp, 7%), and H(p,v) must be 
defined in terms of say E(T=0,v,) or E(T=0, p=0) instead of E(Po, vy). 
However, this only changes the value of the constant in (14), and requires 
no essential revision in the discussion below for Hugoniot curves 
H = constant > 0. It is convenient to suppose that E(p»,v») does exist 
so as to permit simultaneous treatment of the cases H = constant 2 0.) 
The differential of the Hugoniot function (13) is, on using (3), 


dH = dE+3\(ptpo) dv+(v—vp9) dp} 
T dS +3{—(p—pp) dv+(v—w) dp}. (15) 


‘This is the expression used by Weyl (1949) in his analytical discussion of 
the problem. However, it is much more convenient for a geometrical 
treatment to write (15) in the form (Cowan 1953) 


dH = T dS+daA, (16) 


where dA is the area swept out in the p-—wv plane by the radius vector 
(Po. Up) > (Pp, z) in moving from the point (p, v) to the point (p + dp, v+ dv), 
dA being positive if the radius vector moves in a counter-clockwise direction. 
‘That the final term in (15) is the element of area just defined can be seen 
by noting that this term is just half the alternating (two-dimensional vector) 
product of the vectors (U—vt»,p—p») and (dv,dp), or alternatively by 
introducing polar coordinates 


U—Uy = rcos6, p-—pPo =rsing 


which show that this term is 


1,2 do = dA. (17) 


It may be pointed out that although dH is a perfect differential, neither 
T dS nor dA is one. Equation (16) is in fact quite similar to (3), both p dz 
and dA being elements of area in the p-v plane—the one in Cartesian, the 
other in polar coordinates. 

It will be convenient in the discussion below to speak in terms of a 
coordinate system with origin at (Pp, Uv») and so we introduce the coordinates 
p* = p—Ppo, v* =v—wy. The isentropic curve which passes through the 
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origin will be denoted by Sy, and the tangent line to S, at the origin will 
be called 7... For each isentrope S lying above Sy, then, as a consequence 
of assumptions (6) and (7), there exist exactly one point on S in the second 
quadrant and one point on S in the fourth quadrant at which a straight line 
through the origin is tangent to S. ‘The locus of all such points of tangency 
will be denoted by L7, and can be seen to be a continuous curve which is 
tangent to S, at the origin, but otherwise lies entirely above S, (figure 1). 
Because of assumptions (4) and (7), any straight line through the origin 
cuts LT in at most one point, and it is evident from the definition that 
entropy increases monotonically in travelling along L7' outward from the 


origin in either direction. 


S o 


the tangent to S, at the origin 7’; , and the locus LT 


Figure 1. he isentrope S,, 
of all points at which a ray through the origin is tangent to an isentrope. 


The curves L7, Sy, and 7, each pass through the origin but are otherwise 
confined to the second and fourth quadrants. ‘The portions of these curves 
lying in the second and fourth quadrants will be denoted by the superscripts 
(2) and (4), respectively. 

It should, of course, be kept in mind that the region of the (p*, v*)-plane 
in which H is defined is bounded from below by the zero-temperature 
isotherm. ‘This (along with any regions in which the assumptions (3) to (7) 
are invalid) may introduce certain limitations on the validity of some of the 


statements made below. 
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4. PROPERTIES OF THE HUGONIOT FUNCTION 
General properties of H(p*, v*) 
Property 1. H(0,0) = 0, by definition. 


) 
Travelling to the right along an isobar p* < 0 or upward along an 
isochor 7* > 0, both dS and dA are positive, and hence dH > 0) 
from (16). 


PROPERTY 2. For p* < 0, (¢H/dv), > 0; for v* > 0, (CH/op), > 0. 


Property 2’. Wherever condition (11) is satisfied, (¢H/ov), > VU even 
though p* > 0. 
This follows immediately from (15). 


LT ip*™ 








Figure 2. On any straight line through the origin, H increases away from the origin 
on the first-quadrant side of LT’ and increases toward the origin on the third- 


quadrant side of LT (Property 3) 


PROPERTY 3. On any ray through the origin travelling away from the 
origin, ZH > 0 on the first-quadrant side of L7’ and dH < UO on the third- 
quadrant side of L7 (figure 2). Thus H passes through a maximum 

Gd. 

Onsucharay, dA 0, and hence from (16), dH rds. he qualitative 


behaviour of His thus identical to that of S and follows immediately trom 


the definition of L7 














2°77 
Ry) 


ul 


Properties of the Hugoniot function 


Property 4. ‘Travelling along an isentrope in the direction of increasing 
pressure, dH > (0) on the first-quadrant side of LT and dH <0 on the 
third-quadrant side of LT (figure 3). Thus H passes through a minimum 
at LT and through a maximum at LT, 

Along an isentrope, dS = 0, and hence dH = dA. 











a aiid 
ia "ee “Se saa - (4) 


— 


— 


ag ° 


Figure 3. On an isentrope, H increases in the direction of increasing pressure on the 
first-quadrant side of LT and increases in the direction of decreasing pressure 
on the third-quadrant side of LT (Property 4). 


PROPERTY 5. Except for the origin, H > 0 on LT and on S§” and 
H <Qon Sand T,. ‘Travelling away from the origin, dH > 0 on LT@ 
and S\” and dH <0 on SM“, TY and TY’. 

The statements regarding LT and 7, follow from Properties 1 and 3, 
and those regarding S, from Properties 1 and 4. ‘Travelling away from 
the origin on LT, both dS and dA 
Note, however, that a similar statement cannot be made on LT since dS 


and dA have opposite signs ; H 0 on LT. but does not necessarily 


are positive and hence dH (). 


il 


increase monotonically in a direction away from the origin—see the 
5 


discussion following Property 1 


Properties of the curve H = 0 


In the case Of a shock which includes no chemical r action (in the sense 
; ; : ; “4 . 
described in §3), the possible thermodynamic states behind the shock are 


those for which H U: such States are thus of conside¢ rable physical 


interest. 


F.M. 
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PROPERTY 6. There exists a continuous curve H = 0 (denoted by H, 
for brevity) which passes through the origin O, but is otherwise confined 
to the second quadrant between LT and S\ and to the fourth quadrant 
between S\\? and 7'~?. For any point P on S§, the ray O > P cuts H, in 
exactly one point; for P on S\”, O > P cuts Hy in at most one point and 
a vertical or a horizontal line through P cuts H, in exactly one point. 

It is evident from Properties 3 and 5 that H > 0 on the first-quadrant 
side of LT® and Sj”, and H < 0 on the third-quadrant side of S" and 
Ts. ‘Thus any points for which H = 0 must lie in the portions of the 
second and fourth quadrants stated above. For any point P on S\’, the 
ray O-> P cuts LT™ in exactly one point O; since the values of H at P 
and O are less than and greater than zero, respectively (Property 5), and 
since H increases monotonically from P to O (Property 3), the line segment 
P->O contains exactly one point at which H=0. Similarly, for any 
point P on S$” (where H > 0), a vertical or a horizontal line through P 
intersects 7', at a point O (where H < 0)f, and by Property 2 there exists 
exactly one point between P and O at which H =0. Continuity of the 
curve H, follows from that of the curve S, and that of the function (13). 
For P on Sj}, O > P cuts H, in at most one point by Property 3. 

Note that there is nothing which precludes the possibility of maxima 
and minima in volume along H'>); Bethe has, in fact, discussed phenomena 
(molecular rotation, vibration, dissociation, ionization) which tend to 
produce just such effects. ‘The question of pressure maxima and minima 
on H{ is dealt with in Property 10 below]. 


Property 7. H > 0 on the first-quadrant side of the curve Hy, and 
H <0 on the third-quadrant side of Hp. 
This is evident from the statements made in proving Property 6. 


PROPERTY 8. Asa point P moves along H, in the direction from fourth 
to second quadrant, the straight line through O and P turns continuously 
clockwise with no stationary values of its slope. 

In the vicinity of the origin, Hy has positive curvature. On the second- 
quadrant side of the origin, this follows from Property 6 and assumption (7) ; 
on the fourth-quadrant side, H, cannot have negative curvature because 
of Property 6, and it cannot have zero curvature because to the second 
order dS (and hence also dH) is negative going down Ty”. ‘Thus Property 8 
holds in the vicinity of the origin. Since the line O- P cannot cut H, 
in two different points (Property 6), the direction of rotation of O- P 


+ If O lies below the zero-temperature isotherm (isentrope), then the line through 
P must be terminated at the S = 0 curve, where H may or may not be negative 
(Property 4). his indicates the likelihood of an upper limit to volume on Hp. 
Such a limit indeed exists (as pointed out by Bethe) since the final term in (13) tends 
to plus infinity with v whereas E(p, v) certainly has a finite lower bound. 

[ Bethe states that pressure maxima and minima may occur on H“), This 
is a result of his assumption (¢p/éS),, > —2T/v, which is less restrictive than (4). 
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Properties of the Hugoniot function 
can never reverse itself. A point of inflection P at which H, is tangent 


to O + P is also impossible, from Properties 3 and 6. 


Property 9. ‘Travelling along H, in the direction from fourth to 
second quadrants, S increases monotonically (dS > 0) except for a 
stationary value at the origin, where AS is of the third order in shock strength. 

Since dH = 0 along Hy, then, from (16) and (17), 

T dS -dA ty do. 


In the vicinity of the origin 


* _ (6 1/2 
tan’ = Fe ={— (eae ll o* 
ye dv) 9 2.\ dv? a 
and 72 1/22 
dtan@ = —— dé = 5 (sa) dv*, 


so that 





“_ 1) —: 1 (&p «)3 1 (*?\ 
S(v*)— S(O) = — iif. (54) (ey nara 127, (4) (ey (13) 


which is of the third order in v*f. It may be noted that the change of 
entropy along H, near the origin is directly related to the non-zero 
curvature of the isentrope S,. (In fact, if S) were a straight line, then 
along S,, dS = dA = 0 and hence dH = 0; thus Hy would be coincident 


with S, along the entire length of the latter.) 


Property 10. In travelling along H away from the origin, the pressure 
and energy vary as follows: (a) if dv <0, then dp > 0 and dE > 0; 
(b) if dv 0, then dp 0 if and only if the inequality (11) is satisfied 
on H,, and dE > 0 if the inequality (12) is satisfied on Hp. 

If dv < 0 in travelling up H®, then Property 8 implies dp > 0, and 
from (3) 

dE = TdS—pdv 0 
since dS > 0 from Property 9. If dv > 0 in travelling up Hf, then 
TdS > 0 and —p dv 0 and nothing definite can be said about the 
sign of dE. However, if the inequality (11) is satisfied, then (CH/dv),, > 0 
from (15), and this implies dp > 0 from Property 7. If the even stronger 
condition (12) is satishied, then dz 0 and dp 0 along H, imply from 
(8) and (12) that also dE > 0 along H,. (Note that assumption (12) is 


a sufficient but not a necessary condition for dk > 0.) 


Properties of a curve H c (c = constant > OV) 
Constant-H curves for H < 0 are not of great interest since such a 
curve (which will be denoted by H_) represents the locus of possible states 
+ This expression can also be obtained by direct differentiation of (15) at constant 
H to give the coefficients (4,S/év\),, in a ‘Taylor series expansion for S, 
2M 2 
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behind a shock which has induced an endothermic chemical reaction, and 
such cases are not of much physical importance. ‘The properties of such 
curves are similar to those of H, and will only be summarized briefly. 


Property 11. H_ lies everywhere on the third-quadrant side of H,; 
H > —c everywhere on the first-quadrant side of H_ and H < ~c 
everywhere on the other side (Properties 7 and 3). For p* <0, a 
horizontal line cuts H_ at most once; for v* > 0, a vertical line cuts H 
at most once (Property 2). Any ray from the origin cuts H_ at most once 
(Property 3). In travelling along H_ in the direction from fourth to second 
quadrant, dA <0 and dS > 0 (see (16)). On H®, the behaviour of 
pressure and energy is as described in Property 10. 


Properties of a curve H = +c (c = constant > 0) 


Constant-H curves for H > 0 are of great physical interest since one 
such curve (which will be denoted by H_,) includes the locus of possible 
states behind a shock-induced exothermic chemical reaction (as in the 
detonation of an explosive, either gaseous, liquid, or solid) or behind such 
a reaction which is not shock-induced (slow combustions or deflagrations) ; 
see, for example, Courant & Friedrichs (1948). 


Property 12. H., lies everywhere on the first-quadrant side of H,, 
or, if H is not defined at the origin, then everywhere above the zero- 
temperature isentrope. In any case, H/, is confined to the second, first, 
and fourth quadrants. 

If H is defined at the origin, this follows from Property 7; if not, then 
it follows from the physical requirement 7 > 0. In addition, it can be 
shown (Courant & Friedrichs 1948) that the only physically significant 
parts of H, are those in the second and fourth quadrants, and so we shall 
not consider H‘” in any detail. 


PROPERTY 13. Given any point P on LT, then, according as the value 
of H at P is less than, equal to, or greater than c, the ray O > P cuts H, 
not at all, at the one point P, or at (at most) one point on either side of P. 
Similarly, the isentrope through P cuts H. zero, one, or (at most) two 
times according as H(P)c if P is on LT® or according as H(P)  « 
if Pis on LT®, 

‘These statements follow immediately from Properties 3 and4. Note the 
essential difference between these properties of H., and the corresponding 
properties of H, and H_, directly connected with the fact that only H 
can intersect 7 at a point other than the origin (Property 5). 


PROPERTY 14. For any point P on LT), there exists a continuous curve 
H_,, through P, which at P is tangent from below to S{? and tangent from 


~ 9 


above to O- P; (0®p/dv®), > (@p/ev?), > 0. On H, the entropy 
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has a maximum value where it crosses L7 and decreases monotonically 
on either side of this point. 

The statements regarding tangency follow directly from Properties 4, 
3, and 2, and the continuity of the Hugoniot function (13). ‘The entropy 
maximum follows from the tangency relations, and the absence of any 
other stationary points in entropy is obvious since only at points on LT 
can dS and dA be simultaneously zero (see (16)). 

The intersection of H“ and LT represents the state of the reaction 
products for a Chapman-Jouguet deflagration (Courant & Friedrichs 1948). 


* 





|P 





Figure 4. Diagram showing the mutual relation of isentrope, L7\2), and curves H., 
defined by H = constant > 0. Under the assumptions (4) to (10), the relations 
shown at (a), (6), (c), and (d) are possible (Property 15), but those at (e) are 
not (Property 3). Under assumption (11), only (a) is possible. 


Property 15. For any point P on LT®, the constant-H curve, Eee 
which passes through P may behave in any of several ways: (a) H. may 
be tangent from above to the isentrope through P; (6) H, may consist 
of the single point P; (c) H, may be tangent from below to the ray O > P; 
(d) H. may have two branches, one of the form (a) and one of the form (c). 
H. cannot have a point of inflection at P, behaving as in (a) on one side 
of LT® and as in (c) on the other side. 

At any point P on LT, then along O > P (or along the isentrope 
through P) dS = dA = 0 and hence dH = 0) from (16); thus H. must be 
tangent to O-> P (and to the isentrope S,). However, H, cannot lie in 
between O -> P and S, because, travelling away from P on such a curve, 
dS <OanddA < Oand hence dH cannot be zero. ‘Vhus H, must behave 








542 Robert D. Cowan 


as shown at (a), (d), (c), or (d) in figure 4. An inflection point at P is not 
possible; see figure + and Property 

It is evident from figure 4+ (which may be viewed as a contour-map of 
a three-dimensional (v,p,H) surface) that H, has the form (a), (4), (c), 
or (d) according as in travelling up LT® H increases, passes through a 
maximum, decreases, or passes through a minimum, respectively. ‘Thus 
the necessary and sufficient condition that only the relation (a) prevail 
is that H increase monotonically in travelling up LT (see the discussion 
following Property 5). Since on LT (dp/dv), = (ep/ev)y <0, an 
equivalent condition is that (CH/dp), > 0 or (CH/dv), > 0. ‘This last 
inequality, from (15), is equivalent to (11). Jn all the following discussion 
of the properties of H.,, it will be assumed that on LT 


(=) — 3 (p z p ) or (5) iol (v —v) 
Ov 3 z\ . 0 op * 2 0 Ud. 


Proofs of most of the following statements have already been given or are 
obvious. 


Property 16. In travelling along LT in the direction away from the 
origin, dH > 0. 


PROPERTY 17. For any point P on L7™, there exists a continuous 
curve H, passing through P which is tangent from above to the isentrope 
through P ((0%p/ev*),, ( 0*p/cv*), > O) and which meets LT only at P. 
For any point P’ on LT between O and P and for any point P” on LT“ 
beyond P, the ray O-> P’ and the isentrope S». intersect H‘ not at all, 
O-> P and Sp meet H‘ only at P, and O— P” and S,. each intersect 
H® at most twice (once on either side of LT). 


Property 18. In travelling along H® away from LT in either 
direction, S increases monotonically with no stationary values except at 
the point of intersection P of H“’ and LT, where S is a minimum. For 
a point O near P, the difference S(Q)—S(P) is of the second order 
Av = Ug — Up. 


Similarly to the derivation of Property 9, one hast 


tan 0 op ote. 1 wl (Av)* 
; ov H 2 ov" H U* ; 





2 | r? . 0? Av 
_ dtand an av = Use!) age d(Av), 
i pail Cv H re 
S(O) — S(P) = So 5), (Av)?. (19) 
\ov 


PROPERTY 19. At the point of intersection of H‘? and LT, the internal 
energy increases in the direction of increasing pressure. On the part of H‘: 


t See the footnote to Property 9. 
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which lies to the left of LT, the behaviour of pressure and energy are as 
given in Property 10 for H®. 
The first statement follows from (3) and (6); 
(CE/0p)q = (2E/op)s = —p(dv/ep)s > 9. 
Proof of the remainder is exactly as that for Property 10. 
The intersection of H® and LT® represents the state of the reaction 
products for a Chapman—Jouguet detonation (Courant & Friedrichs 1948). 


5. SUMMARY AND APPLICATIONS 
In figure 5 is sketched a series of constant-H curves illustrating the 
most important properties developed above, assuming that condition (11) 
holds for p* > 0. ‘These properties—subject to limitations which may 





Figure 5. Diagram showing isentropes, the curve L7, and the general form of the 
constant-H curves assuming conditions (4) through (11) to be valid. (In order 
to illustrate volume minima on H@) on the scale of this figure, the upper part 
of H‘2) has been drawn in such a way that the requirements of (21) are 
violated.) 

arise from the fact that the zero-temperature isotherm does not pass below 

the origin, or from phase transformations or other causes of violation of 

assumptions (4) to (10)—are briefly as follows. 
The (p*, v*)-plane can be thought of as divided into three parts by the 
half-line (v* = 0, p* <0), and by the curves LT and LT, the loci of 
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all possible Chapman—Jouguet states for exothermic detonations and 
deflagrations, respectively. No constant-H curve has pressure maxima 
or minima anywhere, nor volume maxima or minima for v* > 0. No 
constant-H curve is cut by a straight line through the origin nor by an 
isentropic curve more than once in each of the three regions of the (p*, v*)- 
plane. The curve H = 0) passes through the origin, in the second quadrant 
lies between LT® and S®, and in the fourth quadrant lies between S\ 
and 7. Above and to the right of any constant-H curve H = c, the value 
of H is greater than c; on the other side of this curve, H is less than c. 
‘Travelling up H = 0 in the direction from fourth to second quadrant, 
the entropy increases monotonically except for a point of inflection at the 
origin. ‘Travelling up a curve H = constant <0, dS > 0 everywhere. 
A curve H = constant > 0 in general crosses both LT and LT®; in 
travelling up this curve, starting below LT, entropy increases mono- 
tonically until it reaches a maximum value at L7™, then decreases 
monotonically to a minimum value at LT, after which it again increases 
monotonically. 

Two applications of the properties developed above will be given. 
Both pertain to the curve H}* (the case in which the shock produces no 
chemical reaction), and serve to elaborate on certain characteristics of this 
curve. 

1. The equation of state of a gas may be written in the form 


E = po|(y—1). (20) 


In the case of a monatomic ideal gas, y= 5/3 for all p and v; the effect of 
excitation and ionization in a real gas is to increase EF for given p and v 
and thus produce an effective value of y which is less than 5/3. ‘Thus for 
a point (p,v) which lies on the Hugoniot curve for the ideal gas, it is seen 
from (13) that H > 0 for the real gas. From Property 7, it follows that 
the Hugoniot curve for the real gas lies to the left of (or below) the ideal-gas 
curve. At very high pressures, the real gas is completely ionized and the 
potential energy becomes negligible relative to the kinetic energy, so that 
the equation of state of the real gas becomes identical with that of the ideal 
one and, consequently, the two Hugoniot curves also become identical; 
the possibility of a volume minimum on the Hugionot curve of the real 
gas (due to ionization, etc., as mentioned earlier) is thus evident. 

2. For a certain problem which the author once encountered, it was 
important to know for an arbitrary point (p,, 7,) on H' whether the isentrope 
through this point intersected the chord (pp, v») — (pj, @;) or lay everywhere 
above it. It is immediately evident from Property 6 that the former is the 
case. In fact, the point of intersection is closer to (pp, vp») than to (p,, %) 
as can be seen in the following way. In figure 6, S, is the isentrope 
through (p,,7v,) and (ps, v2) is the point on S, for which py = py. Then 
from the Hugoniot relation (1) 


E, — Ey = 3(P — Po)(Yo — U1) + Poo — 21), 
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and from (3) 


so that 


E,— Ey = A, —Az—po(v2—%); (21) 


where 4, and A, are the shaded areas in figure 6. From assumption (10), 
it follows that 4, > Ay. Thus it may be seen that not only does Hj” hie 
between LT® and S, but tends to lie closer to the latter than to the former 


(see figure 5). 






Pp 
te 
“2 (paw) 
PY, 
ltr. 2 2 
S, 
0 
a Vv 
Figure 6. Relation of the Hugoniot chord (po, vo) — (pi, v;) to the isentrope through 
(pi, ¥y); the area A, is greater than the area Ap. 
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Viscous Flow Theory, II—Turbulent Flow, by Suin-I Pai. Princeton: 
D. Van Nostrand Co., 1957. 277 pp. $6.75 or 48s. 

In a time which most people would find barely adequate for the writing 
of one book, Professor Pai has published no fewer than three on different 
branches of fluid mechanics. Beginning with Fluid Dynamics of Fets (1954), 
he followed up with Viscous Flow Theory, [—Laminar Flow (1956), and 
now we have this new work on turbulent flow. ‘The two parts of the 
treatise on viscous flow theory are really separate books, since there is not 
much in common, pedagogically speaking, between laminar and turbulent 
flow; indeed, there is not even much of a case for the use of a common 
title in view of the relatively unimportant part played by viscosity in most 
cases of turbulent flow. It is no mean achievement for a man to find the 
determination and knowledge needed for three books in rapid succession, 
and it might be supposed that they have a purpose which is obvious. I 
have not in fact been able to discern clearly the purpose served by the 
present book or to classify it as one of the common types. I do not think 
that it can be called a research monograph on turbulent flow; the author 
does not reveal important new advances, and his book does not have the 
physical insight and critical appraisal that would be needed for it to say 
anything new to those in touch with the subject. Nor can it be said to be 
an exposition for students; the material is too raw and incomplete, and 
includes too many unrelated items of current research which are of unknown 
significance, for the book to be suitable for teaching purposes. ‘There is 
undoubtedly a need for a textbook on turbulent flow, but the present work 
does not have the regard for fundamentals and for logical coherence that 
such a text should have. Despite the claim made on the dust-jacket (with 
characteristic punctuation and spelling) that ‘‘ Engineers, research scientists 
as well as advanced students will find this complete treatment of turbulent 
flow an indispensible source of fresh information”’, it seems to me to lack 
both the indispensability of a good textbook and the freshness of a good 
research monograph. 

So tar as technique of writing is concerned, the author has acted 
simply as a reporter, absorbing large masses of paper and regurgitating it 
in a uniform style in the manner of a journalist. ‘The present book, like 
a good newspaper, can fairly be said to be comprehensive in scope, up to 
the minute, impartial, and discursive. It also has the defects of some 
journalistic reporting—a limited respect for the truth, haste in composition 
and printing, some irresponsibility, and an apparent reluctance to think 
about the material for longer than it takes to read it and write it 
down in different (not always!) words. Professor Pai does have a 
genuine flair for reporting scientific developments; as many visitors 
to the University of Maryland have found, he makes a creditable 
job of writing up their lectures for later reproduction and distribution by 
the University. ‘This kind of writing requires mental agility and elasticity 
to an unusual degree. But though the rapid reporting of many different 
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scientific subjects demands considerable skill, the job encourages super 
ficiality and this is unfortunately the hall-mark of the book under review. 
Passing lightly from one snippet of news to another, the text always stops 
short of the point of illumination and fails completely to provide a clear 
picture of the basic mechanics of turbulent flow. 

‘The scope of the book is wide, and embraces all aspects of turbulent 
motion; the classical cases of turbulent shear flow, mixing-length theories, 
boundary layers with pressure gradients, compressible turbulent flow, 
turbulent diffusion, decaying isotropic turbulence, spectrum theory, 
Kolmogoroff’s similarity theory, and even turbulence in an electrically 
conducting fluid, all find some kind of place in the book. ‘he treatments 
of these topics are not thorough (nor could they be in a book of this size), 
and a more cautious author would have hesitated to try to compress a whole 
subject like turbulent diffusion into 18 pages. But granted that it is useful 
to say something about everything within the space of one book, there is an 
acute problem of selection of material to be faced. ‘lhe pieces of work 
which Professor Pai chooses for inclusion show him to be attracted strongly 
by those having a mathematical form and hardly at all by physical thought ; 
for instance, the book includes fairly full accounts of Reichardt’s *‘ inductive 
theory of turbulence” (pp. 37, 129), the method of determining profile 
drag from wake measurements (p. 135), Goldstein’s hyperbolic diffusion 
equation (p. 186), correlation and spectrum tensor relations (pp. 214-230), 
and Chandrasekhar’s ‘ propagation’ of density fluctuations (p. 256), but 
on the other hand makes no mention whatever of what is perhaps the 
most fundamental (and practically useful) formula in the theory of 
turbulence, namely, that the rate of dissipation of energy is of the 
same order of magnitude as u//, where u and / are representative 
velocity and length scales of the turbulence. ‘lhe more straightforwardly 
mathematical investigations in any subject are undoubtedly the easiest to 
comprehend, and are often seized upon with relief by students; for that 
very reason, it is desirable aiways to explain their significance, the assump- 
tions on which they are based, their implications, and their agreement with 
observation. I am sure Professor Pai would agree with that precept, but 
I do not think he has done much about it; he says (p. 209), “‘... we have 
no general methods for solving such nonlinear equations. At the present 
time, we try to obtain an understanding of the mechanism of turbulence 
by studying experimental results ...”’, but in fact he does not provide a 
single diagram or table of observations relating to homogeneous turbulence 
even though about half the book is concerned with that case. His bias to- 
ward highly theoretical enquiries sometimes leads him into rather pretentious 
and irrelevant generalities, as for instance when (ch. 8) he drags in phase 
space and Liouville’s theorem and other concepts of statistical mechanics. 

The book falls into two non-overlapping and non-interacting parts, 
the first being concerned with ‘semi-empirical theories of turbulence” 
and the standard cases of turbulent shear flow, and the second with the 
‘‘statistical theory of turbulence”’ and isotropic turbulence in particular. 
‘This dichotomous approach to the subject of turbulent flow was remarked 














Reviews 


548 


in a recent review of Momentum Transfer in Fluids by Corcoran, Opfell 
and Sage (7. Fluid Mech. 2, 1957, 204); indeed, had Professor Pai’s book 
appeared a little earlier, it would have been convenient to review them 
together since they have a good deal in common quite apart from their 
concern with turbulent flow. ‘The division is made quite sharply here, the 
mean velocity distribution in jets, wakes, etc. being regarded as belonging 
to the “‘semi-empirical theories” and velocity fluctuations in those same 
flows as belonging to the “statistical theories”. Unless someone soon 
writes a book which integrates the old and the new, the means and the 
fluctuations, the “‘ semi-empirical” and the “‘ statistical”’, we shall have a 
generation of students who believe that there must be some basis in nature 
for this dual approach and that jets and wakes really do have two faces. 
Mixing-length theories and the statistical terminology represent the contri- 
butions of two separate periods of research (both of them now past), and 
we ought by now to be able to make up our minds about what is good in 
each and to allow what is not good to be quietly forgotten. 

The first part of the book on the “semi-empirical theories” follows a 
familiar and appropriate pattern. ‘There is first a chapter introducing the 
momentum and vorticity transfer theories and hypotheses about the mixing 
length, and containing a rather perfunctory account of the facts that they 
are intended to explain. ‘The two key empirical laws concerning the mean 
velocity are the universality of the distributions in the central region of a 
channeland near arigid planeboundary. ‘Ihesetwolaws deserve prominence, 
as any number of writers have pointed out, but Professor Pai passes over 
them lightly and does not even distinguish them clearly (p. 28). The next 
five chapters consider in turn the cases of channel flow, a flat plate without 
and with a pressure gradient, a compressible boundary layer, and what the 
book calls ‘“‘turbulent jet mixing regions and wakes”’. | cannot see here 
any improvement on the corresponding parts of the books by Schlichting 
(1950) and Goldstein (1938), any advantage conferred by the date of 
the present book being more than balanced by the greater perspicacity 
of these earlier writers. ‘The author perpetuates the misleading half-truth 
that near a wall ‘‘there is a laminar sub-layer in which the flow behaves 
as a laminar flow” (p. 40) and that the relative velocity fluctuations fall to 
zero at the wall (p. 242). He makes the calculation of turbulent boundary 
layer parameters seem even more like a set of arbitrary rules than it usually 
does, and even though the chapter on the compressible boundary layer is 
a brave attempt to be up-to-date, he has not succeeded in giving the various 
contributions to that problem any order and coherence. In his description 
of free turbulent flows, the wonderful phenomenon of the sharp transition 
from rotational to irrotational flow at an irregular fluctuating boundary 
is mentioned only once (p. 75) and it is there referred to as a postulate. 

The second part of the book, on the “statistical theory”’, is evidently 
more congenial to the author, and he makes full use of the greater 
opportunities for presenting mathematical analysis. He has had the useful 
idea of including an account of elementary probability theory, and he has 
taken particular trouble over the account of harmonic analysis of random 
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functions. He uses the box device of first imagining the random function 
to be zero outside some distant boundary in order that the Fourier coefficients 
should exist, and of allowing the boundary to recede to infinity at an appro- 


priate stage of the analysis. ‘lhe approach has the advantage of directness, 
although in Professor Pai’s hands it is rather lengthy (and he surely unnerves 
his readers when he says (p. 217)—wrongly—that the shape of the box 
in three dimensions could affect the results after the limit is taken). He 
makes a number of portentous remarks about ergodicity in general and 
about the relation between theoretical and experimental averages in 
turbulence, but he fails to say clearly (indeed it is not evident that he knows 
it) that stationarity in the variable concerned in a necessary, and usually 
in practice a sufficient, condition for equality of the ensemble average and 
an integral average. ‘lhe whole treatment of stationary random functions is 
poor, and there is confusion (pp. 161-164) between the various values taken 
by a random variable and the various values taken by a random function of ¢ 
at different values of ¢ in one realization of the function. ‘There are several 
gaffes, as when he states that an auto-correlation function is positive definite 
(p. 213) and that non-normality of the probability distributions of velocity 
derivatives indicates that “they are not entirely random in nature”’. 

After these mathematical preliminaries, most of the standard analysis 
of homogeneous turbulence is presented. Provided no physical arguments 
are involved, Professor Pai gets by without much error, although the 
selection of material is curious and undiscriminating. Where physical 
justification of the analysis is needed, as in the account of Kolmogoroff’s 
similarity theory, and in the brief excursion into turbulence in a conducting 
fluid, the text is quite inadequate. ‘The introduction to chapter 13 reflects 
confusion between the various wave-number ranges of the spectrum and 
the various stages of decay; it may have been this confusion which led 
Professor Pai into the persistent error that inertia forces have no effect on 
the smallest eddies (p. 234) and that as a consequence the energy spectrum 
has a Gaussian form at large wave-numbers (p. 238). Nowhere does one 
hear about the law of decay of turbulent energy and the spectral distribution 
that are actually observed in real wind tunnels. As for the chapter on 
turbulent diffusion, I can say only that evidently Professor Pai and I have 
quite different views about what are the real problems involved. 

I am sure that my general opinion of this book has already revealed itself 
sufficiently. Like the book by Corcoran, Opfell and Sage, I do not think it 
can be recommended, either as a text for postgraduate students or as a 
source of information or enlightenment to more senior people. Again like 
this earlier book, it is written in rudimentary English, and with little regard 
for the rules of grammar and punctuation. 

As a less sober postscript to this review, 
of the present work (Ff. Fluid Mech. 2, 1957, 
Pai’s diverting habit of incorporating slices of other people’s papers in his 


I, like the reviewer of Part I 
515), was struck by Professor 
book. This is perhaps carrying the reporting technique a little too far. 
I was able to recognize three whole paragraphs and some odd sentences 
from two of my own papers in the chapter on diffusion, and more of my 
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own paragraphs in other chapters (with no hint that a direct quotation 
was being given and often without reference to the source). He even uses 
this strange technique on his own writing, and on page 146 reproduces 
what he has given on page 10. How many other writers are represented 
in this book? As a method of constructing a book it is so absurd that one 
cannot be indignant at being copied in this way. Scissors and paste may be 
useful accessories at a lower level of writing, but at this level it is impossible 
to integrate existing pieces of prose into a continuous reasoned narrative. 
An author cannot avoid mastering a subject and rewriting it himself. 


G. K. BATCHELOR 


Magnetohydrodynamics, by J. G. CowLinc. New York: Inter- 
science Publishers, Inc., 1957. 115 pp. $3.50. 

The motion of a fluid or a gaseous medium possessing electrical 
conductivity in the presence of a magnetic field has become of great 
importance in the last decade. ‘The theoretical treatment is variously 
designated as ‘magneto-hydrodynamics’ or ‘magneto-gasdynamics’, or 
more generally ‘magneto-fluid dynamics’, and sometimes, léss character- 
istically, as ‘hydromagnetics’. Although a few experimental investigations 
on the influence of magnetic forces on the motion of a conducting liquid 
date from the previous century and the discovery of the magnetic field of 
sunspots was made in the first decade of the present century, it has taken 
some time for a wider interest to develop. One feature leading to this 
wider interest seems to have been the discovery of interstellar polarization 
of light, which for its explanation required the assumption of interstellar 
magnetic fields of wide extension. Problems connected with cosmic 
radiation led to the same conclusion. Once the presence of such fields was 
recognized, many other problems called for investigation. Almost at the 
same time the wide application of shock waves for the investigation of high 
temperature gas motion with its consequent ionization led to the problem 
whether such flows could be influenced by means of magnetic fields; 
research on the behaviour of magnetically driven plasmas arose out of this. 
Perhaps two landmarks in the development of the subject were the first 
Symposium on Cosmical Gas Dynamics (Paris, 1949, later followed by 
meetings at Cambridge, England, 1953, and Cambridge, Massachusetts, 
1957); and H. Alfvén’s book Cosmical Electrodynamics, published in 1950. 
From 1950 onward papers on astrophysical problems and on results of 
experimental work have come out in rapid succession; the number of 
conferences is multiplying, and now technical applications are inview. ‘Thus 
there is a need for a treatise which may help interested scientists and engineers 
to become acquainted with the many aspects of magneto-fluid dynamics. 

Professor Cowling’s book is not intended as a treatise, but as a condensed 
survey. ‘The topics treated are briefly as follows. Chapter 1 gives the 
basic equations and general ideas, stressing the distinction between cases 
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where the electrical resistance of the medium is a determining influence, 
and those cases where the conductivity of the medium can be considered 
as being practically infinite. ‘lhe meanings of the magnetic Reynolds 
number and of other non-dimensional quantities are pointed out. Since 
one of the factors occurring in the magnetic Reynolds number is the length 
scale of the field, the case of ‘infinite conductivity’ is of great importance 
in astrophysical and geophysical systems with their large dimensions, the 
decay times for electric current systems being here of cosmical duration. 
As an illustration of a different case a simple treatment is given of laminar 
flow under the combined influence of ordinary viscosity and what is 
sometimes termed ‘magnetic viscosity’. Chapter 2 is about ‘magneto- 
statics’ and considers equilibrium states in which the mechanical forces 
are balanced by magnetic stresses. Applications are made to some 
astronomical problems (sunspots, streamers and filaments, and spiral arms 
of the galaxy). Chapter 3 treats wave motion, mainly for an incompressible 
fluid where the velocity of propagation is wholly determined by magnetic 
effects (the so-called ‘ Alfvén waves’, connected with the effective tension 
along the lines of magnetic flux). Again some astrophysical applications 
are mentioned. Chapter + considers the influence of magnetic forces on 
the stability or instability of fluid motion, with examples taken from cases 
of laminar (incompressible) flow and from problems of convection. 
Chapter 5 is entitled ‘ Dynamo theories’ and considers the question whether 
fluid or gas motions in the presence of an originally weak magnetic field 
can lead to an increase of the field strength, in the way in which this occurs 
in a self-exciting dynamo. Elsasser’s, Bullard’s, and Parker’s theories of 
the origin of the magnetic field of the earth are considered and criticized. 
The chapter also includes a brief discussion of turbulent fields of motion 
and of the distribution of energy between turbulent motion and magnetic 
field. Finally chapter 6 gives attention to the effects of molecular structure, 
as they appear in ionized gases. Each chapter is concluded with a short 
list of the most important and relevant publications. 

From this summary it will be seen that Cowling’s book is of value to all 
who wish to become acquainted with various topics of magneto-fluid 
dynamics. But on the other hand, considered from the point of view of 
the needs of this moment, the book leaves many things to be desired. 
Although it is evident that in the short space set for the Interscience ‘Tracts 
on Physics and Astronomy (of which this book is one) there is no room for 
the writing of a textbook, one has the feeling that in preparing it the author 
has limited himself mostly to a critical discussion of a series of selected 
topics, in order to point out wheretheoriesare incomplete. Professor Cowling, 
who is an authority in the field of gas motion under the influence of magnetic 
fields, of course has full right to do so, but as a result of this the reader is not 
greatly stimulated to enter into the subject and to attempt to work out 
problems for himself. It is clear moreover that the emphasis in the selection 
is on the geophysical and astronomical applications, as Professor Cowling 
himself observes in his preface; although he adds the hope that the 
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publication of the book may expedite developments also in the engineering 
aspects, one cannot say that he is giving much help in this respect. 
Chapter 6, for instance, on ionized gases, indeed gives the relevant 
equations and puts them into various forms, and these forms are discussed 
from a number of points of view, but at the end the reader still may ask: 
where do I stand? How can I apply this to a particular problem? Mean 
collision times are introduced, but there is no indication concerning the 
methods in use for calculating or estimating them (although some of these 
methods are considered in The Physics of Fully Ionized Gases, by L. Spitzer, 
Jr., which appeared earlier in the same series of tracts and which may in 
a sense be regarded as a companion volume). In some of the other chapters 
of Professor Cowling’s book, in the treatment of problems concerning the 
sun, the magnetic field of the galaxy, the magnetic field of the earth, 
and turbulence theory, the ultimate impression obtained is rather a negative 
one; it looks as if almost nothing fruitful can be done in these matters. 
Problems concerning the behaviour of acoustic waves, and still less 
problems referring to shock waves in a magnetic field, have not been 
considered at all. Here a great many cases are awaiting adequate treatment 
and it would have been useful if at least some indications had been given. 
There are, for instance, three-dimensional problems concerning the 
propagation of shock waves approaching a magnetic field (in experimental 
situations one always has to do with non-uniform fields); there is also the 
problem of a shock wave in a cold gas subjected to a magnetic field, the 
electrical conductivity appearing only behind the wave front as a consequence 
of the high temperatures generated there. ‘There is the problem of 
stationary supersonic flow of an ionizable gas around a magnetized body. 
Boundary layer flow under conditions leading to ionization and with 
magnetic field influencing the motion likewise calls for attention. Still 
another group of problems refers to plasmas driven by magnetic fields ; 
this is under investigation at present on account of the very high temperatures 
which can be producedinthem. ‘The magnetic field due tothe plasma current 
itself, as well as fields derived from outside sources, are both of importance 
in guiding the motion of the plasma. ‘lhe question of the stability of the 
plasma current involves the so-called ‘ pinch effect’ (the lateral contraction 
of a plasma current as a consequence of its own magnetic field); Professor 
Cowling does make two short references to this effect (pp. 27, 29) without, 
however, making any reference to the intriguing questions connected with 
it. (Some discussion of the ‘ pinch effect’ is to be found in Spitzer’s book.) 
Certainly we are grateful for the material Professor Cowling has put 
so lucidly in 115 pages. ‘The book seems well adapted to the needs of 
astrophysicists and geophysicists who wish to obtain a rapid view of the 
present state of the application of magneto-fluid dynamics to their field. 
But one who comes from the aerodynamical side still feels the need for 
a work of more comprehensive character, which can serve as a guide for 
those who wish to develop their ability to work in this interesting domain. 


J. M. BurGERs 
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